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This work discusses charmless inclusive B decays and its application to the ex- 
traction of \Vub\- Starting from first principles we relate the differential decay rate 
to the hadronic tensor in terms of optimal choice of kinematical variables. We 
review the traditional methods of calculating the hadronic tensor, expansion in as 
and HQET, and discuss their shortcomings. In the kinematical region relevant for 
experiment ("shape function" region), the hadronic tensor can be factorized, at 
each order in l/rrib, as a product of calculable hard functions and a convolution of 
calculable jet functions with non perturbative shape functions. 

Using SCET, we calculate the leading order hard and jet function to first order 
in as- Large logarithms are resummed in RGE improved perturbation theory. 
Local OPE is used to relate moments of the renormalized shape function to HQET 
parameters, defined in the "shape function scheme" . Beyond leading order in 1/m;,, 
several subleading shape functions arise. We derive them at tree level, where they 
can be expressed as forward matrix elements of bi-local light-cone operators. 

Based on these theoretical calculations we present two applications. In the first, 
we present the "state-of-the-art " expressions for the triple differential B — > X^li^ 
decay rate and the B — > Xs^ photon spectrum. These expressions include all 
known contributions and smoothly interpolate between the "shape-function" and 
" OPE region" . Based on these an event generator can be constructed, from which 
the theoretical prediction for any experimental cut can be extracted. In the second, 
a weight function is constructed that relates the P+ spectrum in 5 — > W to the 



normalized B — > Xg'y photon spectrum. This weight function is independent of 
the leading order shape function. At leading power the weight function contains 
two loop corrections at the scale ^/mbAqcB- Power corrections from phase space 
are included exactly, while the remaining subleading contributions are included at 
first order in Aqcd/^6- 
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Chapter 1 
Introduction 

1.1 Preface 

The standard model is a theory that describes our knowledge of three of the four 
fundamental forces in nature: the strong, the electromagnetic and the weak. The 
strong and the weak forces especially pose interesting challenges to theoretical 
physics. 

The strong interaction can be described at high energies by the QCD La- 
grangian, where the quarks couple to each other via an exchange of gluons. The 
coupling strength is universal, i.e. independent of the flavor of the quarks. At en- 
ergies above roughly 1 GeV we can use perturbation theory to calculate the effects 
of the strong force. At lower energies the perturbative description, based on the 
use of the interaction Lagrangian, breaks down. As a result, at low energies one 
has to use either non perturbative methods such as lattice calculations and chi- 
ral perturbation theory, or extract information from experiment. The separation 
of scales and the perturbative from the non perturbative physics will be a major 
theme in this work. 

The weak interaction is one of the richest parts of the standard model. As 
opposed to the strong interaction, the weak interaction is not universal and is 
different for different flavors. More specifically the strength of the interaction 
between quarks of different flavors is described by a unitary matrix known as 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix. The CKM matrix elements are 
fundamental parameters of the standard model and therefore it is important to 
measure them as accurately as possible. Indeed, Measuring the CKM matrix 
elements has been one of the main goals of research in particle physics, both for 
theory and experiment. 

This work aims at precision measurements of Vub, one of the smallest CKM 
matrix elements . Currently the highest precision on \Vub\ can be achieved through 
inclusive measurement of charmless decays of B mesons. By charmless we mean 
both semileptonic (B — > Iv) and radiative [B — > -^^57) decays. The former 
is used to actually measure |Ktf,|, while the latter is used to extract, directly or 
indirectly, the main non perturbative object which is the shape function. Of course, 
better understanding of radiative decays is important by itself, as it has important 
implications for constraining physics beyond the standard model. 

The outline of this work is as follows. We begin by an introduction, which 
takes the reader from "first principles" to the starting point of chapters 2 and 3. 
These two chapters are more theoretical, where in 2 and 3 the perturbative and 
non perturbative calculations are presented, respectively. The last two chapters 
are concerned with applying these calculations. Chapter 4 gives the necessary 
ingredients for a construction of an event generator, which allows experimenters 
to calculate any semileptonic spectrum they would like to measure. In chapter 5 a 
weight function is constructed, that directly relates the P+ spectrum in i? — > li' 
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to the photon spectrum in S — > Xg^. We present our conclusions in chapter 6. 

1.2 Chapter outline 

In this chapter we begin our study of inclusive charmless B decays. We start with 
the description of the kinematics of these decays and discuss the optimal choice 
of kinematical variables. We then turn to the dynamics. The basic quantity we 
are trying to calculate is the hadronic tensor, which cannot be calculated exactly. 
We review the "traditional" approximation methods, namely, expansion in and 
Heavy Quark Effective Theory (HQET) and discuss their shortcomings. In order 
to overcome these problems we introduce the appropriate effective field theory suit- 
able for inclusive charmless B decays, the Soft CoUinear Effective Theory (SCET). 
This effective field theory will be used in the following chapters to analyze both 
the perturbative (chapter 2) and non perturbativc (chapter 3) corrections to the 
hadronic tensor. This introduction is meant to be pedagogical, so we present most 
of the derivations explicitly. 

1.3 Kinematics made easy 

In this work we focus on B ^ Xu Iv decays, but also discusses B — > decays, 
both by itself and in relation to B ^ X„ Iv. These decays are inclusive in a sense 
that one of the final states is something that contains a specific quark {u or s) but 
is otherwise unrestricted and can even be a group of particles. 

We begin analyzing these decays by looking at the simple question, "How many 
kinematical variables are needed to describe an inclusive event?" . We are looking 
at a decay oi&B meson into n particles, where n — 1 of them have known mass. The 
"particle" X has unspecified invariant mass. In general for a decay into n particles 
there are 3n — 4 kinematical variables, but not all of them are determined from the 
dynamics. Since we do not have information about the spin of the hadronic state 
X or the spin/polarization of the other decay products, the only four vectors we 
have at our disposal are Pg, the four momentum of the B meson, pi = Px the four 
momentum of X, and p2, . . . , Pn- Lorentz invariance implies then that the possible 
variables are the (n + l)(n + 2)/2 scalar products of these n + 1 vectors. Since the 
masses of all the particles apart from pi are known, we have n constraints of the 
form pf = mf, (i = B,2,. . . ,n). Conservation of momentum would allow us to 
eliminate all the pairs that contain Pb (since Pb = XliLiP*)' ^-^ eliminating n + 1 
variables. That leaves us with n{n — l)/2 independent variables. The rest of the 
variables are undetermined from the dynamics and can be integrated over. This 
argument holds for n < 3, since at 4 dimensions at most 4 vectors (corresponding 
to n = 3) can be linearly independent. For n > 3 the argument needs to be 
modified. 

The decay mode B — > Xg'y has n = 2 and therefore only one relevant variable. 
This variable is usually taken to be E^, the energy of the photon. As we shall see. 



a related variable, P+ = Mb — 2E^ will turn out be useful. For B — > IP we 
have n = 3 and three relevant variables (see table 1.1). There are various choices 
in the literature for these variables, and in order to understand them we have to 
say a few words about the dynamics of these decays. 

Table 1.1: Number of kinematic variables for inclusive B decays. Using the con- 
straints on the right column we can eliminate some of the variables on the left 
column. 
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First we need to distinguish between the hadronic level, which looks at the 
decay as a decay of hadrons, and the partonic level, which looks at the decay as 
a decay of quarks. At the hadronic level we have a B meson carrying momentum 
Pb = Mbv decaying into a leptonic pair (the lepton and the anti-neutrino) carrying 
momentum total q, and a hadronic jet carrying momentum Px- Conservation of 
momentum implies Mbv — q — Px- At the partonic level we look at this decay as a 
decay of a 6 quark, carrying momentum rUbV, into a u quark, carrying momentum 
p (at tree level), and a virtual W carrying momentum q. The W in turn decays 
into a lepton I and an anti-neutrino z/;. (More accurately we write the momentum 
of the b quark as rribV + k where k is O(Aqcd)- and expand in powers of k). 
Conservation of momentum implies rubV — q = P- If we define A = Mb — mij we 
find that Px = p + Av. Beyond tree level p would be the momentum of the jet of 
light partons created in the decay. We will elaborate on the exact definition of A 
in chapter 2. 

There are generically two common choices of variables: 

• Leptonic: The energy of the lepton Ei, the energy of the neutrino E^, and 
the invariant mass of virtual W boson q^, which is also the invariant mass of 
the lepton pair. This choice focuses on the leptons created in the decay and 
was used for example in [1] 

• Partonic: the energy of the lepton Ei, the energy of the partonic jet v ■ p, 
and the invariant mass of the partonic jet p^. This choice focuses on the 
partons created in the decay and was used for example in [2]. 



We will use a different choice of variables, which in some sense is the optimal 
one. We start by noting the fact that any measurement of B ^ X„ W suffers 
from a large background from the decay B — > Xf-lv. The reason for the large 
background is that |Vc6| 3> \Vuh\i meaning that the h quark "likes" to decay to a 
c quark much more than to a m quark. In order to eliminate the background we 
have to look at regions of phase space where particles containing charm cannot 
be produced. For example we can only look at events for which Pj- < M|,. This 
fact implies that the "typical" state will have large energy of order 7715, since it 
originates from a decay of a heavy quark and intermediate invariant mass, because 
of the experimental cut. Since ~ ^foAgcD we can write P|- ~ /WftAqcD- This 
imphes that some of the components of Px are larger than others. We would hke 
the choice of our variables to reflect that. In other words wc look for two variables, 
one of which scales like and one of which scales like Aqcd- A pair of variables 
which satisfy this condition is: 

P-^Ex+ \Px\ , P+^Ex- \Px\ , (1.1) 

where Ex and Px are the energy and momentum of the hadronic jet, respectively. 
Note that P+P_ = P+ + P_ = 2Ex and P+ < P_. The scahng of these 
variables will therefore be P_ ~ nib and P+ ~ Aqcd- In order to complete our 
choice we should add another variable. Specifying P_|_ and P_ would determine the 
energy and the invariant mass of the lepton pair, but not the individual energies 
of the lepton and the neutrino. We therefore have to add another variable that 
would distinguish between the two. We will choose: 

Pi^Mb- 2Ei , (1.2) 

where Ei is the energy of the lepton. These three variables allow us to specify any 
B — > Xu Iv event. 

Having chosen the variables, the next task is to determine to phase space in 
terms of P+, P/, and P_. In the rest frame of the B meson, conservation of energy 
and momentum gives us the following equations (the notation is self explanatory): 

Ei + E, + Ex = Mb (1.3) 
K + K+Px = 0. (1.4) 

We take the leptons to be massless, implying Ei = \Pi\ and Ey = \Pi,\. The limits 
of phase space are determined by the extremal values of the angles between the mo- 
menta. Because of conservation of momentum (1.4) we have only two independent 
angels. 

Let 6 be the angle between P; and P^. Equation (1.4) then implies P^ = 
Ef + El + 2EiEyCos9. Since -1 < cos^ < 1, we have {Ei - E^f < \P x? < 
{El + E^y. Using P+ < P_, the upper limit gives us P_ < Mb, and the lower 
hmit P+<Pi<P-. 
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Let a be the angle between Pi and P x- Equation (1.4) then imphes El — 

[Me-Ei-Exf = Ef + \'Px\^ + 2Ei\'Px\cosa. Since -1 < cosa < 1, the upper 
hmit gives {Mb - - Pi) > 0, and the lower hmit (Mb - P-){Pi - P+) > 0. 

We have also an additional constraint from the QCD spectrum. The lightest 
state containing a u quark is the pion, so we must have < Pjr = P+P-. 
(Without this constraint the lower hmit of P+ would be 0) . Combining all of these 
constraints we finally have: 

-p^<P+<Pi<P-< Mb. (1.5) 

One of the benefits of this choice of variables (apart from the easy derivation) is 
that phase space has an extremely simple form, probably the simplest form possible 
(compare for example the phase space limits in [2]). This has practical implications: 
since we ultimately integrate over phase space, some of the integrations can be done 
analytically. We note that a similar choice of variables would be useful also for 
W. 

What are the phase space limits for B Xg'-jl Conservation of momentum 
and energy imply that Mb — Ex = E^ and E^ = \ P x\-i oy P_ — Mb and P+ = 
Mb — 2E^. The limit of the phase space for P_,_ are determined from E^ > 0, and 
the QCD spectrum constraint Af|.. < P| = P+P- = MbP+, where X*(892) is the 
lightest strange meson that can be produced in this decay. All together we have: 

§^<P.<Ms (1.6) 

Notice the similarities between equation (1.5) and (1.6), especially if we take the 
heavy quark limit, which sets the lower limit of P+ to be 0. In the following we 
will explore these similarities extensively. 

Having discussed the kinematics let us say a few words on experimental cuts and 
the regions of phase space they probe. For measurement of the photon spectrum 
B — > Xsj, experiments must impose a cut on the photon energy of about 1.8 GeV 
(measured in the rest frame of the B meson) [3, 4, 5]. This means that the typical 
values of P+ = Mb — 2Ej are of order Aqcd- 

For B — > XuW, as we mentioned earlier, experimental cuts are imposed to 
eliminate the charm background (figure 1.1). Some of the cuts are: 

• Cut on the charged Icpton energy Ei. If Ei > (M| -M|,)/2Mb ^ 2.31 GeV, 
the final hadronic state will have invariant mass smaller than Md- For this 
cut Pi < 0.66 GeV, which implies that P+ is of order Aqcd 

• Cut on the hadronic invariant mass Mx ■ To eliminate the charm background 
we need M^ < M|, . The cut M^ = M|, is depicted as a solid black hyperbola 
in figure 1.1. 




P_ [GeV] 

Figure 1.1: The (-P+, P-) phase space for B Iv. The charm free region is the 

dark grey region below the black hyperbola, which correspond to M'j^ = P^P_ = 
Ml. The sohd blue line is = {Mb - P+){Mb - P-) = (Mb - Mof. The red 
dashed hue is P+ = M^/Mb- 



• Cut on the leptonic invariant mass g^. Any cut of the form < (Mb — Mbi)'^ 
would not contain charm events. The cut = {Mb — P+){Mb — P-) = 
{Mb — MoY is depicted as a solid blue line in figure 1.1 

• Cut on P+ will be discussed in chapter 2. The cut P+ = M^/Mb is depicted 

red dashed line in figure 1.1 



1.4 From the optical theorem to the hadronic tensor 

The basis for the calculation of an inclusive rate is the optical theorem, which 
relates the decay rate to the imaginary part of the forward scattering amplitude. 
In terms of the T matrix [6] 

{f\m = {2t:Y5\P, - Pf)M{i ^ /), (1.7) 

the optical theorem states that 

/^3 p 

\M{i - f)\' (27r)V(P. - P;). (1.8) 



Using the optical theorem we obtain the following expression for the total rate 

^ = 2]ki:/2|lFi^('-^'iw(/'.-/'.) 

^ lmM{i^i) (1.9) 



This formula holds both for B — > Xu W and B — > Xg'^, but }A is different in each 
case. 

In order to proceed we need to find the matrix element of the T matrix, which 
is determined by the Hamiltonian. The relevant general Hamiltonian contains 
both weak interactions and strong interactions (unless otherwise noted, we neglect 
contributions from electromagnetic interactions) . For the weak interactions we can 
use perturbation theory. At the lowest order in the weak interaction coupling 
(or alternatively Gp) we have ior B ^ X^W decays: 

Hw = ^VubU^'^il - 75)&^^Tm(1 - + h-c. (1.10) 

This Hamiltonian can be thought of as an effective Hamiltonian, obtained (at tree 
level) after integrating out the heavy W bosons [7]. For B — > Xg'f we have (see [7] 
for derivation) 

Hw = ^C,^^Vt,V*m-, (^^^ sa^'^il + ^,)F,,h + h.c, (1.11) 

There are other contributions from the weak Hamiltonian to i? — > ^sT, but their 
contribution to the rate is suppressed by gl and we will ignore them. 

Using S = Tie-'^'^"'"'^} (T is the time ordering symbol) and = 1 + if, the 
non vanishing lowest- order contribution to the T matrix element is given by: 

{B\f\B) = {B\'- j d^xid^X2T{Hw{xi), Hw{x2)}\B). (1.12) 

If we write = + 0^, the only non vanishing contribution to the matrix 
element would come from T{0\ O}. Using the delta function in equation (1.7) 
we can eliminate one of the Xi integrals to obtain the result: 

M^iJ d^x(B|T{O^(0), 0{x)}\B). (1.13) 

This can simplified further, by noting that up to a constant O can be written 
as a product of two "currents" : 

O = CqVbJNH = C J Jnh- (1.14) 

For B — i> XuW transitions J = u-^^{l — 75)6 and Jnh = ^7i/(l — 75)^';, and for 
B — > Xgj we have J = sa^"{l + 75)6 and Jnh = 2d()Ay. In both cases the fields 
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in Jjvif do not have any strong interactions, which would allow us to factorize M. 
to a product of hadronic and non hadronic tensors. We insert a complete set of 
states in the form of 

(The notation is such that the sum over the hadronic states includes the integration 
over the hadronic momenta) Thus we get: 



M = i d^xlC 



NH 



{B\j\0)\X,){0\4M\NH){X,\J{x)\B){NH\JMm 

d^PNH 



V AT LJ *^ 



{B\J{x)\Xhbg){0\JNH{x)\NH){X,,g\j\Q)\B){NH\4Hm^^ 

(1.16) 

where Xg (Xhi^) is a complete set of hadronic states containing a q quark (two 
b quarks and a q). NH is a non hadronic state, i.e., a photon for B — > Xg^ 
and a lepton pair for B — > Xulv- The second term in equation (1.16) would not 
contribute to B decays, since after taking the imaginary part one gets a factor of 
5^{Pb — Px^iq — q) (see appendix A). For the first term we distinguish between 
B Xgj and B Xy, IP- For B Xgj we have: 

J2{Q\4^imH){NH\JMm = ^J2^oKMQ)MQ)){^M\dfiMm 

NH A 

= -iqe^q^g^.e^''', (1.17) 
where q is the photon momentum. For B — > Xu W we have: 

J2(^\4HmNH){NH\JMm = 

NH 
spins 

8 (P^P;' + P^P; - g,.P' ■ P^^ - te,.afsPrP.^,) e^'\ (1.18) 

where q = P^ + P'^^ and we take eoi23 = 1- 
Using this we write A4 as: 

^ = E/^|^^V(S|r-|S>. (1.19) 



Q 



The current correlator Tf^" is defined as 



T""" J d^xe^«^r{jt'^(0), rix)}, (1.20) 

where the currents are: 

JLs = «-^[fl,7l(l + 75)&, (1.21) 

and we have written the photon momentum as — n ■ (see below). We have 
also defined: 



^2 



ub\ 



2 



and 



' 6— >u 



8 (p^p;' + Pip;^ - g,.p' ■ p'^ - ie.^pPrp^;) 



LIT = -Ag,^{n-qf. (1.23) 

We should note that in writing M. in this form we have used the fact that the second 
term in equation (1.16) does not contribute to B decays, so we have changed its 
phase factor to e**^. 

It is customary to define the hadronic tensor W^^ as: 

TT 2Mb 

which allows us to write the differential rate as: 

= E J^IL. 27rA Viy-. (1.25) 

In order to continue we separate the discussion of S — > and B — > IV. 

For B decays, we have a real photon (g^ = 0) recoiling against the 

hadronic jet. In the rest frame of the B meson we can take the jet to be moving in 
the +z direction and the photon to be moving in the —z direction. The momentum 
of the B meson in its rest frame is then Pb = Mbv, where v'^ = (1, 0, 0, 0). We also 
introduce two conjugate light-like vectors n'^ — (1,0,0,1) and n'* = (1,0,0,-1), 
which satisfy n-n — 2, n-v — n-v — 1, n-\-n — 2v. Any 4-vector can be 
decomposed as: 

= n ■ a h n ■ a \- a'^ = a't + a'^ + a^, (1-26) 



n 



where a'^ are the x,y components of (for our choice of n and n). In terms of 

op 

2 ~ ^-^T 2 



these vectors wc can write the photon 4-momentum as q'^ — n ■ q^ — 2E^^. We 



can also express P+ and P„ as P+ = n ■ Px, P- = n ■ Px- 

For P — > the Patj^ momentum in (1.25) is q. Integrating over the angles, 
we find that the photon spectrum in S — > Xs^ is: 

^ 27rAViy^'^ = I^IK^V;:! Waemty;, (1-27) 

We now turn to P — > X„ /i/ decays. The P^h momenta in (1.25) are Pi and Pj.^. 
As we saw before, of the 6 variables only 3 are relevant, which can be taken to be 
El = {Mb-Pi)/2, = {Mb+Pi-P+-P-)/2 and cos^ = (Pj,-Ef-El)/{2EiE,) 

{9 is the angle between P; and P^). Integrating over the three irrelevant angles 
and changing variables to P^,Pi,P-. We find that 



dV = dP+dPi ^p_ ^^ l^ "'-i'(^- — ^L.,W^"'. (1.28) 

The hadronic tensor has the property (proved in appendix A) that W^^ = {W'^'^)*. 
Using this identity we can decompose the tensor as a sum of 5 possible Lorentz 
structures. Given an independent pair of 4-vectors a and b, the 5 structures are 
^n^u^ jjujju^ (^^i^jju ^ a''h^'),g^"' and -^e^"'''PaJ}fj. The coefficients of these structures 
are called structure functions, denoted as Wi. 

Typical choices for a and h that appear in the literature are v and q [1], and v 
and p [2] . We will choose v and n, which has the benefit that the coefficients of the 
structure functions in the triple decay rate are independent of mf,. Decomposing 
the hadronic tensor in this basis we write: 

-g'"'W2 + v''v''W3 + {n'^v'' + 7fv^)Wi + n^rfWr, (1.29) 

The reason for defining the coefficient of Wi in this way is that then at tree level 
only Wi is non zero. The structure functions are (generalized) functions of 
P+ and P_ (or n • p and n ■ p). Having decomposed the hadronic tensor, we can 
write the triple rate in terms of the structure functions Wi. In the expression 
LfiuW^^ we have several scalar products: v ■ Pi,v ■ Pui,n ■ Pi,n ■ Pi,^, Pi ■ Pi,^. Using 
Pui = MbV — Px — Pi we can eliminate all the scalar products that contain P,y^ . Of 
the remaining scalar products v ■ Pi = Ei and n ■ Pi can be obtained from Px ■ Pi, 
via Px ■ Pi = [P-(n ■ Pi) + P+(2P; — n ■ Pi)]/2 (since Px has no _L components). 
After a httle bit of algebra we find the following compact result: 



d'r Gl\Vuk\' 



dP+ dP^ dPi 167r3 



{Mb - P+) 



(P_ - Pi){Mb -P. + Pi- P+) m + {Mb - P-){P- - P+) ^ 
+ {P_-Pi){Pi-P+)(y^W; + W^ + ^W^], (1.30) 



1 n 



where we have defined 



y = ^. (1.31) 



Notice that the triple rate only depends on one linear combination of W?,, and 
VF5. The reason is that in the limit M; — > (justified for electrons and muons) we 
have the two constraints q^L^^, — q^L^^ — 0. 

Equations (1.27) and (1.30) are the main results of this section. Of course they 
are not useful without the knowledge of W^^ . The (approximate) calculation of 
the hadronic tensor will be one of the main topics of this work. In the rest of 
the chapter we discuss mostly B — > X„ W and mention in passing what are the 
equivalent formulas for B — > ^^7. 

1.5 Free quark approximation - tree level and as correc- 
tions 

In calculating the hadronic tensor, the simplest approximation we can make is to 
look at the h quark as a free quark, ignoring the spectator quark in the B meson 
and the non perturbative interactions that binds them together. As we shall see in 
the next section, this approximation would be the zeroth order term in an HQET 
based local Operator Product Expansion (OPE). 

In this approximation the B meson states can be written as \B) = \h{mi,v)). 
Writing the current as J'^ = qVb, the matrix element of T'^" is: 

(K™..)|r'.>M> = ./.'.e-(6|r{.no).rw}|6> 

^^^e^(™^>^^26(mbv)r'^r{g(0),g(x)}r'^ii6(mbv) 

b quark spin 

= -m, J d'p5\m,v - q-p)Tr {r^^r'^i±i} , 



(1.32) 



where in the second line we have averaged over the b quark spins and F'* = ryO-piJ.\^o_ 
Using the identity [8] 

lm-{x±ier' = ^^^<)>-i)(a;), (1.33) 

we find that ^ 

^^^TiS.r'':^^^^] . (1.34) 



(Notice that for the free quark approximation rrib ~ Mg). For B X^lv decays 
we can write 5{p^) as 5{n ■ p)/{n ■ p). Since 

5{n ■ p) = 6{n ■ p) In- p— + n ■ p— j = 5{n ■ p) in- p— \ , (1.35) 



1 1 



and r'* = 7''(1 - 7^) we find 

W^"" = S{n-p){n''v'' + n''v'' 
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Wi = S{n-p). (1.36) 



This simple result is the reason that Wi is defined as it is. 

As a quick check let us plug this result into (1.30). Since P+ = n • p + A, we 
can integrate over the total phase space to find the well known result, that total 
rate at tree level is Gp\Vub\^ ml/ (1927:^). 

For B X,7 decays we have f = i[7^, ^](1 - 75) and = 7''](1 + 75). 
Calculating the trace in equation (1.34) we find that 



-2 6{n-p). 



(1.37) 



Using n-p = n- {mi,v — q) = mb — 2E^, we can plug this result in equation (1.27) and 
find that the total B X,-f rate at tree level is Gl\VtbV*^\'^\C^^\'^a^^ml/{32n'^) 

A much more complicated calculation is to find the 0{as) corrections to the 
structure functions Wi, in the free quark approximation. Such a calculation was 
preformed by De Fazio and Neubert in [2]. The hadronic tensor in [2] was decom- 
posed using the vectors p and v: 



{p^v" + p'^v^' - g^'^v ■ p - ie^^°Vt'/3)W"i 



;i.38) 



The relation to our basis is: 
^^^n^p-n^ 



W2^W2 + n-pWi 



W3 = W3 + 2n-p{Wi + W4 + n ■ pW^) 
n ■ p — n ■ p 



;i.39) 



n ■ p — n ■ p 



The kinematical variables z,t,p^ and x in [2] are related to ours via: 

z = ifi • p + n • p)/mi„ t = {n ■ p — n ■ p) / {n ■ p + n ■ p), fP' = {n ■ p){n ■ p)/ml, and 

X = {Mb — Pi) /nib, where again = n ■ p + A, P_ = n ■ p + A. 

We will not transform the complete 0{as) result from [2] to our basis. Instead, 
we would like to illustrate the well known problem with the QCD calculation, 
namely the appearance of double and single logarithms that might spoil the con- 
vergence of the perturbative expansion. Let us look on one of the terms of Wi: 



6{p') + 



An 



-4 



Inp^ 
p2 



[1] 



+ . . . 



where the "star" distribution is defined in chapter 2 and the ". 
0{as) terms. In terms of our basis we have: 



;i.4o) 



denote other 



Wi = S{n ■ p) + 



An 



-4 



n ■ p — n ■ p 
n ■ p 



ln(n ■ p/rub) 



n ■ p 



[rrib] 



+ 



(1.41) 
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The star distribution has two parts. The second part, which is proportional to a 
delta function, is the source of the double log. Inserting equation (1.41) into (1.30) 
we have: 



ub\ 



|2 



dP+ dP_ dPi 16n^ 



{Mb - P+) (P- - Pi){Mb -P- + Pi-P^ 



47r 



S{n-p) ^ 
2 rrib 



(1.42) 



In a somewhat implicit notation, u is the upper limit of n • p in the P+ integration 
and depends on the specific experimental cut. As an example, let us find the P/ 
spectrum. The integral over P_ is automatic. For the P+ integral wc note that 
the upper limit of n • p is u — Pi — A. Changing variables to x — {Mb — Pi) /nib we 
find that: 



dT ^ Gj,\Vub\'ml 
dx 1927r3 



2x^(3 - 2x) |l - ^ [ln^(l - x) + . . . ] I (1.43) 



which was first derived in [9]. It means that for high values of i?/, i.e. a; — 1, 
the so called "endpoint region", the double log gets large. Notice that equation 
(1.42) implies that this is a generic problem. Anytime the value of u is not of 
order m;,, the double log would not be small. But, because of the experimental 
cuts the values of u are never of order rribl A similar situation arises for the single 
logarithms. Since these logarithms are multiplied with ag, we might question the 
convergence of the perturbative series. 

The same double logarithms appear for B — > Xg^. Prom [10] and earher 
references cited there, we find that the decay rate for B — > with a lower cut 
on the photon spectrum > (1 — 5)m.i,/2 is: 

r [B - x.t],,,,.^.,„./, = ^\v.K,m-fml {l - [Wi + ...]} . 

(1.44) 

Because of the experimental cut 6 ~ Aqcd/''^6; so we encounter large logarithms 
again. (The universality of the double logarithms is not a coincidence and has its 
origin in the properties of soft Wilson loops, see below). 

Having large logarithms is of course a well known problem, encountered, for 
example, when one takes into account the QCD corrections to weak decays [7]. The 
solution to the problem is to "resum" these logarithms. This is done by matching 
to an appropriate effective field theory and then solving the RG equations. As we 
shall see, the appropriate effective field theory is SCET. 

One more comment before we conclude this section. Usually when one talks 
about "large logarithms" the question arises "how large are the logarithms? do 
we need resummation?" In order to answer this question we should note that 
the scale dependence of as{^) is unspecified at this order. Some prefer to set 
II ^ rrih, but as we shall see, the more appropriate scale for these logarithms 
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is ~ ^/rn}JVqp£) ~ 1.5 GeV. Recall that because of QCD running as{Mz) — 
0.112, ^^(mt) = 0.22, «s(1.5GeV) ~ 0.375. So we have the effect of logarithms 
multiplied by large value of a^. This further motivates the need for resummation. 

1.6 Effective field theories I: HQET 

Effective field theory (EFT) is a very useful tool in general and in particular for 
B physics. It allows us to extract as much perturbative physics as possible from 
the matrix elements, and to break down the problem of "solving QCD" into man- 
ageable pieces, that we can try to model or extract from experiment. 

For our specific goal, the calculation of the hadronic tensor, given an effective 
field theory the procedure is as follows: 

• A QCD current J is matched onto a series of the form Jqcd = CqJo + 
Ci Ji + . . ., where the Wilson coefficients Cj are calculable as an expansion 
in as and each Jj is suppressed by i-th powers of a large scale. 

• The QCD Lagrangian is matched onto a series of the form C = Cq-\-Ci-\- 

• The \B) meson states are matched onto eigenstates of Cq, called |-B)o . 

• T^^^ can then be written as a series 

t{ E ^^Jl E ^^-^^ e^-^'^'^^i+-} = t[c,JI Co Jo} + 

i i 

+t[c,jICoJo} +T^^CoJlC,J,} +t[Co4,CoJo,i J rfVi} + ... 

• W'^'^ is then written as a series of non perturbative matrix elements with 
calculable coefficients. 

The effective field theories that we will use HQET and SCET, are based on an ex- 
pansion in inverse powers of the b quark mass. Since the hadronic matrix elements 
scale like Aqcd and Aqcd << fT^-b, our expansion parameter will be Aqcd/^^^ ~ 0.1. 

We start by looking at HQET and review its predictions for the hadronic tensor. 
In order to do so we go over some of the elements HQET that are important for 
our purpose. For an extensive review of HQET see [11, 12]. We start by writing 
the b field as: 

b{x) = e-*'""''-" [K{x) + H,{x)] , (1.45) 

where 

K{x) = e^"*''^-^^^6(x), H^(x) = e^"^'''"^^^6(x). (1.46) 
In terms of these fields the QCD Lagrangian for the b quark is 

jC = b{ip - mb)b 

= hyiv ■ Dhy -Hy{iv D + 2mi,)Hy + hJpH^ + HyipK, (1.47) 
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where — D^^ — v^v ■ D. From equation (1.47) it is clear that the field Hy is to 
be integrated out. In zeroth order we find that 

Co = Kiv ■ Dhy. (1.48) 

For the light quarks and gluons there is no change, so we can write the leading 
order current as: 

qVbix) = e-'"''''-''qT''K{x). (1.49) 

Armed with this information we can calculate W'^'^ at tree level: 

= - / d*pd\mbV -q- p)hS^ o ^ r'^fe^ (1.50) 
J + ie 

In order to find the matrix element of T^'^ we need [11]: 

{B\ h^FK \B) ^ Mb Tr\r^-^\ . (1.51) 



Proceeding as before we find that the zeroth order HQET result is identical to the 
free quark result, namely: 

Wi^5{n-p), W^ = -2 5(n-p) (1.52) 

The as corrections to this zeroth order result would give us the De Fazio-Neubert 
hadronic tensor. 

The reason for using HQET is to calculate the 1 / corrections to the hadronic 
tensor. Such a calculation was performed in [1, 13] for B — > XclP. We will present 
the results as they appear in [12] after taking the limit rric — > 0. These corrections 
are known only at 0(q;°). 

First let us look at the 1/mj, term in the HQET Lagrangian: 

Ci = -hy- — hy - ghy-^ hy (1.53) 

zm^ 4m6 

The matrix elements of the operators in Ci define two HQET parameters: 

{B\hyD^hy\B) = -2MbAi 
{B\hya^,G^^hy\B) = -12MbX2 (1-54) 

Intuitively we can think of Ai as measuring the kinetic energy of heavy quark 
inside the B meson and of A2 measuring its chromomagnetic interaction. Both 
parameters have mass dimension 2 and a priori we expect Aj ~ ^qcd- 
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Instead of matching the current onto HQET current, the current correlator 
T^^'^ itself can be matched onto HQET. Doing so one finds that the zcroth order 
term is identical to the free quark result. The l/m^ correction to T^*^ vanishes. 
The reason is that the only operator that can appear at this order is v'^hyiD'^hy. 
At order 0(q;°) the matrix element of this operator would vanish because of HQET 
equations of motion [14]. This implies that the first non zero correction appears 
at order l/ml. 

The hadronic tensor in [12] is defined slightly different from ours: 

^'^^ = "^27^^"^ / ^{■^i(^)' (1-55) 

where = ?i7^(l — 75)6/2. This hadronic tensor is decomposed using the vectors 
q and v: 

= -g^^Wi + v^v^W2 - ie^yocfiv'^q^W^ + q^^q^W^ + {y^q^ + v^q^) (1.56) 
The relation to our basis is: 

^W^ + Wi+-W^ = {n-p-n- p) ( ^ 21^3 ) (1.57) 



4 y ym^ — n ■ p 

After the change of basis we find, including the zeroth order result Wi = 5{n- p), 

2Ai-3A2\ ,/2Ai-3A2 5Ai + I5A2 



3n • p Grrih 



- S"{n-p)^, 



W2 = 5{n-p) 



4Ai - (3A 



3{n ■ p) 



2 



y~ , ~ l~ 5(n-p) /2A1 + I2A2 4A1 + 9A2 



4 y n ■ p \ 3n ■ p Sm^ 

^ M(H^^4A.). (1.58) 
n ■ p \ 6 J 

As we now show this calculation poses a problem when applied for charmless 
semileptonic decays. Define A ~ Aqcd/'"^6; using this small parameter we can 
define a systematic power counting. For example would be an order 1 quantity, 
while the HQET parameters Ai, A2 are order A^. 

For B — > X(.lP both n • p and n • p scale like m?,. This means that 5{n ■ p) 
and its derivatives are order 1 quantities. Considering the scaling of the HQET 
parameters, we have an expansion of the form ~ 1 + A^ + . . . . We can imag- 
ine extending this calculation to higher order in A, by introducing more HQET 
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parameters. Each term in the expansion would be suppressed by more powers of 
A. 

The situation is fundamentally different for B — * Iv- We recall that because 
of the kinematical cuts enforced on us due to the large charm background, n ■ p 
never really scales like mi,. In fact, it scales like Aqcd- The appropriate power 
counting for B — > X„ lu is therefore n-p ~ A and n-p 1. Since the n-th derivative 
of S{x) scales like we have three types of terms in (1.58): terms which scale 

like A~^, e.g. 6{n ■ p),6"{n ■ p)Xi; terms which scale like A*^, e.g. 5'{n ■ p)\i\ and 
terms which scale like A\ e.g. 5{n ■ p)\i. It follows that the expansion in (1.58) 
must be rearranged as: 



Our notation is such that in each square bracket we group terms which scale the 
same way. (Notice that in W2, for example, the lowest order terms are of order A). 
Even after this rearranging we have a problem. By matching to higher orders in the 
HQET expansion and by introducing more HQET parameters, some of the terms 
might not be suppressed at all and will contribute to lower orders. In fact, in each 
order in A we have to sum infinite number of terms that contribute at that order. 
(We have included ... in each square bracket to reflect this fact). This problem 
was first observed in [15, 16, 17] and is often referred to as "breaking down" of 
the OPE. A more accurate description would be to say that we have tried to use 
an HQET based OPE in a kinematical region where it is not valid. This OPE 
assumes that all the components of p are large compared to Aqcd, which is true 
for B ^ XcW but not for B Iv- 

Let us look at the lowest order in A, i.e. the order A^^ in equation (1.59). The 
infinite number of terms that are summed give us a new non perturbative object 
which is no longer just a matrix element of an HQET operator, but a function 
of the variable n ■ p. This function is known as the leading order shape function, 
denoted as S. Beyond leading order there are more than one shape function and 
one would talk about several subleading shape functions. This analysis is the 
subject of chapter 3. 



5{n ■ p) 



S"{n-p)^ + ... 
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The fact that we had to do this resummation imphes that HQET is not the most 
appropriate effective field theory to describe charmless inclusive B decays. We will 
introduce the more appropriate theory (SCET) in the next section. As chapters 
2 and 3 will show, in this effective field theory the various shape functions arise 
naturally as non local matrix elements and there is no need to resum the HQET 
OPE. For the rest of this section let us talk a httle more about the leading order 
calculation. 

At 0{a^) the leading order shape function can be defined as the sum of a series 
of distributions: 

S{u;) = 6{uj)-^6"iuj) + ... (1.60) 

D 

(When no confusion can arise we will refer to the leading order shape function 
simply as the "shape function"). Moments of this function would then be related 
to HQET parameters. In order to find the relation, let us write a function f{x) 
as a sum of infinite series of the form f{x) = '^^aiS^^\x). Using "integration by 
parts" [8] we find that the n-th moment is J dxx"'f{x) = (— l)"'n! a^. Equation 
(1.60) implies that the first three moments of the leading order shape function are 
1, 0, and — Ai/3. 

In [16] it was shown that we can write the leading order shape function as: 

{B\KS[.-in.D)K\B)^^^^^ (1.61) 

2A/g 

Notice that the shape function is defined as a matrix element of a non local op- 
erator. Another important property derived in [16] is the support of the shape 
function. In the absence of radiative corrections (see chapter 2) we can think of 
the shape function as a probability distribution for the b quark to have residual 
light cone momentum co = n ■ k. For the total b quark light cone momentum 
fraction we have 

n-pB~ Mb ~ Mb ' ^ ' ' 

Since < < 1, (the limits corresponding to the b quark carrying none or 

all of the B meson light cone momentum), we have — < A;+ < Mb — nib- 
In the 777.6 —> oo limit this correspond to — oo < n ■ k < A, where at tree level 
A = limm^^ooiMB — rrib)- In chapter 2 we will define A beyond tree level. 

So far we have ignored the as corrections. How does the picture change when 
we include them? The answer was given by Korchemsky and Sterman in [18] who 
argued that at leading order a factorization formula holds. The factorization for- 
mula states that the differential rate can be written as a product of "hard" function 
and a convolution of a "jet" function with the leading order shape function, where 
both the hard and the jet function are calculable in perturbation theory. We will 
prove this factorization formula using SCET in chapter 2. At this stage we can 
illustrate this formula by using, yet again, the free quark approximation. Let us 
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write Wi as: 



(0) 



d{n ■ p) = 1 ■ J du) n ■ p d[n ■ p {n ■ p + cu)] d{u)) (1.63) 

We now identify the term 5{ou) with the first term in the moment expansion of 
S{u). We also define a hard function as Hui{n ■ p) = 1 + 0{as) and a jet function 
as J{u) = S{u) + 0{as). The factorization formula is then: 



Wi — Hui ■ J dun ■ p J[n ■ p {n • p + uj)] S{uj), 



(1.64) 



or symbolically Wi ■ J ^ S. 

We can motivate this procedure in the following way. At tree level we can think 
of the jet function as related to the imaginary part of the u quark propagator. Not 
neglecting the residual momentum k of the b quark, the light quark propagator is: 



in ■ p— ■ -. (1.65) 

2 n ■ p [n ■ p + n ■ k) + le 

Expanding this propagator as a power series in n ■ k would generate the leading 
order shape function S(n ■ k). Following this procedure we can "implement" the 
shape function in any parton level calculation as: 



Wi = f{n ■ p, n ■ p) ^ du f{n ■ p,n ■ p + u)S{u) 



;i.66) 



Such procedure was used for example in [2]. We should stress though, that this 
procedure is justified at 0{a^) only. A rigorous implementation of the shape 
function at 0{as) appears in chapter 2. 

For B — > Xg^ we find similar features. The HQET calculation preformed in 
[19] gives: 



wi: = -2 



5{n-p)- y5"(n -p) 



+ 



Ai — 3A2 
2mh 



S'{n ■ p) 



:i.67) 



Because of the experimental cuts n -p scales hke A (see section 1.3). We therefore 
have to group the terms according to their scaling in an analogous way to S — > 

Xu Iv. We also see that at leading order the hadronic tensor is proportional to the 
leading order shape function S. Including radiative effects the W^^ (and therefore 
the photon spectrum) can be factorized at leading order as 



Wj^^H^- j du mi, J[mb {n-p + u)] S{u). 



(1.68) 



Comparing equation (1.64) to (1.68), we see that the factorization formulas are 
very similar. The shape and jet function are the identical, but the argument of the 
jet function is different, since for B — > Xg^, n ■ p — rrib (this is a result of the fact 
the kinematics is simpler for radiative decays). The hard functions are different, 
with turning out to be more complicated. We will utilize these similarities in 
chapters 4 and 5. 
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Table 1.2: Different modes in SCET-I. We define A = Aqcd/''^6 and write the light 
cone components of vector as (n • a, n • a, a^). 



Mode name 


Momentum scaling 


Examples 


Hard 


(1,1,1) 


hard gluons 


Hard-CoUinear 


(A,1,a/A) 


hard-coUinear quarks (^) 






hard-collinear gluons {A^^ 


Soft 


(A, A, A) 


heavy quarks {h.^), spectator quarks (g^) 






soft gluons (^45) 



1.7 Effective field theories II: SCET 

In the previous sections we have encountered two problems in trying to calculate 
the hadronic tensor. On the perturbative side we had large double and single 
logarithms that threaten the convergence of the perturbative expansion. On the 
non perturbative side we had to sum infinite numbers of HQET matrix elements 
into non perturbative functions, called shape functions, which arc matrix elements 
of non local operators. This summation procedure is defined at tree level and it is 
unclear how to include radiative correction in a systematic fashion. Both of these 
problems can be solved by using SCET, which we now introduce. 

We will present a brief review of SCET, focusing again on its apphcation for 
charmless B decays. For a more complete discussion see [20, 21, 22]. The relevant 
effective field theory is known in the literature as SCET-1. Our presentation of 
SCET-1 is based on the "position space" formalism first presented in [23, 24]. In 
SCET-1 we distinguish between three type of modes according to the scaling of 
their momenta: hard, hard-collinear, and soft. Two of these modes, namely hard 
and soft, are familiar from HQET. The new mode, hard-collinear, is designed to 
describe the light quark produced in the h quark decay. The various modes are 
listed in tabic 1.2. 

In constructing SCET, we first integrate out the hard modes, which leaves us 
with hard-collinear and soft modes only. For the soft quark modes the Lagrangian 
is the regular QCD Lagrangian for the soft quarks and the HQET Lagrangian 
for the heavy quarks. (The Yang-Mills Lagrangian for the soft and hard-collinear 
modes can be found in [23, 24]). We will therefore focus on the Lagrangian for 
hard-collinear quarks. In order to construct it we decompose the hard-collinear 
fields as: 

V^(x) = ^^+^V'^e + r/, (1.69) 

where {i/),Tfi) / A and are projection operators and n, n are the light like vectors 

introduced before. In terms of ^ and r\ the QCD Lagrangian for massless hard- 



on 



coUinear quarks is 



Chc^i>ipi^^i^in-Di + fi'!^in-Dr] + lip^r] + nip±i- (1-70) 

(The covariant derivative in (1.70) contains both hard- coUinear and soft gluons). 

In order to decide which field is to be integrated out, let us look at the two 
point function for ^: 

(0|r{«x).f(0)}|0> = /|^e--|L||. (1.71) 

Since p is hard-colhnear momentum wc have d^p ~ \\ ~ A, and n ■ p ~ 1 , which 
imphes that ^ scales hkc A^/^. In a similar way we can show that i] scales like A 
and it should be "integrated out". Once this is done equation (1.70) becomes: 

Ac = e I Mri • + ip^^-L-zpA e (1-72) 
2 \ m ■ U J 

This Lagrangian contains terms of different orders in A and we will expand it to 
the lowest order, shortly. 

What about the currents? Our naive expectation is that the current q{x)Th{x) 
would be matched, at tree level and leading order in A, onto e~'^™'^'"'^^{x)T hy{x). 
The problem is that when a soft field interacts with a hard coUinear field only one 
component of the momentum is conserved, namely n • p. As a result, in position 
space, the soft fields need to be "multipole expanded" [23]. At leading order this 
amounts to placing ,^ at a; while is at a;_ = |n ■ xn. As soon as we put the two 
fields at different points the current is not gauge invariant. In order to make it 
gauge invariant we introduce a Wilson line: 



W 



P exp J dt n ■ Ahc{x + tn)^ (1-73) 



The Wilson line has the right transformation properties such that the current 
^-tmtv-x ^^W^ {x)T hy{x^) is gauge invariant. We can think of this Wilson line 
as summing infinite number of hard-colhnear gluons that can interact with the 
hard-colUnear quark at leading order. (Of course each interaction "costs" us one 
power of g which is perturbative at the hard coUinear scale). Apart from this 
hard-coUinear Wilson fine it is also convenient to introduce a soft Wilson fine: 

S^Pex.p(ig[ dt n ■ As{x + tn)) (1.74) 



Using this Wilson line we redefine the hard-colhnear fields according to 

^x) = S{x^)e'\x), Al = S{x.) (x) S\x^). (1.75) 



The fields ^(°) and are "sterile" in the sense that they do not interact with 
soft gluons. In term of the A'^f'\ the hard-coUinear Wilson line can be written as 
W{x) = S{x-)W^'^\x)S^{x-), where 

= P exp (^ig j dtn- Af^{x + tn)^ . (1.76) 

Thus the QCD current q{x)Th{x) is matched onto g"*™''^'^ where 

(1.77) 

and X = ^^'^^W^^^ = ^WS. Beyond tree level and at leading order in A the current 
would be matched onto 



q{x)Tb{x) e-^"*"^-^ [ dsQ{s)X{x + sn)T,:K{x_). (1.78) 

Notice that non local currents arise naturally in SCET. 

We can now go back to the hard-coUinear Lagrangian, equation (1.72), which 
we want to expand to the lowest order in A. By looking at the two point function 
for the hard-collinear and soft gluons, one can show that they scale like their 
momentum. That is, A^^ ~ (A, 1, ^A) and ~ (A, A, A). Since iDf" = id^" + 
d^hc + 9^s- ^^^^ lowest order in A 

in ■ D — > in ■ = in ■ d + gn ■ A^f. 

in ■ D — > in ■ L>hc + gn ■ As — in ■ d + gn ■ A^^ + gn ■ Ag. (1-79) 
The leading order coUinear Lagrangian is therefore: 

^hl = e I ( m • L'hc + • ^(x_) + if±hc ._ ipi.hc] e (1-80) 
z \ m ■ L>hc J 

Note that we have multipole expanded the soft gluon since it interacts with collinear 
quarks. We now introduce the "sterile" fields into the leading order Lagrangian. 
Using the identity 

S\in-d + gn- A^S^in-d (1.81) 

the Lagrangian becomes: 

= l''' If-- ^i? + ^Mt-^^^it) i'''. (1.82) 

Thus at leading order there are no interactions between soft gluons and hard- 
collinear quarks. 



oo 



This concludes the leading order "derivation" of the SCET-I Lagrangian and 
currents. Equation (1.78) will be used in chapter 2 to derive the factorization 
formula Wi ^ H ■ J ^ S. Renormalization group equations will be used to resum 
the large logarithms, solving our first problem. In chapter 3 we will use higher 
order corrections to the currents and Lagrangians to calculate the subleading shape 
function contributions. The shape functions will be described as matrix elements 
of non local operators that arise naturally in the process of integrating out hard- 
collinear modes, so we do not need to sum infinite number of HQET parameters. 
Since SCET is an effective field theory it is clear how to carry on this procedure 
to any order and include radiative corrections, leading to a well defined double 
expansion in as and nib. 
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Chapter 2 

Perturbative Corrections 

2.1 Chapter outline 

In this chapter we analyze the perturbative corrections to -B ^ Iv hadronic 
tensor at leading order in l/rrih. We begin by presenting a derivation of the factor- 
ization formula for inclusive rates in the shape-function region using the position- 
space formulation of SCET. We then perform a two step matching calculation 
QCD — > SCET — > HQET at next-to-leading order in perturbation theory, extract- 
ing the hard (H) and jet (J) functions. We then derive renormalization-group 
equations governing the dependence of the functions H, J, and the shape function 
S on the renormalization scale, and solve these equations analytically in momen- 
tum space. Next, we derive several model-independent properties of the shape 
function S. In particular, we present the precise form of the relations between 
renormalized shape- function moments and HQET parameters such as A and Ai, 
and we give an analytical formula for the asymptotic behavior of the shape func- 
tion. An unexpected outcome of this analysis is the finding that the shape function 
is not positive definite, but acquires a negative radiative tail at large values of |a;|. 



2.2 Factorization theorem for inclusive decay rates 

Recall that 

^ -Im (-^l^n-^) ^ T^'' f d^xe'i'' T{J^''{0),r{x)}, (2.1) 
71 2Mb J 

where for B l^v decays = -07^(1 — 75)6. 

Below a matching scale fi^ ~ the semileptonic current can be expanded as 

^ r ~ - 

u{x)-f>'{l - -t^)h{x) = e-^"^"^-" Y I ds Ci(s) X{x + sn) Tf J{(a;_) + . . . , (2.2) 

i=i 

where the dots denote higher-order terms in the SCET expansion, which can be 
neglected at leading power in Aqcd/'"^6- The position-space Wilson coefficient 
functions Cj(s) depend on the variable s defining the position of the hard-coUinear 
field. Here and below we denote functions in position space with a tilde, which 
is omitted from the corresponding Fourier-transformed functions in momentum 
space. A convenient basis of Dirac structures in (2.2) is 

ri'^ 7^^(1-75), r^ = ^;''(l + 75), r^ = n'^(l + 75). (2.3) 
The current correlator in (2.1) then becomes 




xT{j{(0)r^^X(tn),X(a; + sn)rj'J{(a;_)} + ... . (2.4) 
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In a second step, the hard-coUinear fluctuations associated with the hght-quark 
jet can be integrated out by matching SCET onto HQET at an intermediate scale 
/^j ~ ^mfeAqcD- At leading order the SCET Lagrangian (when written in terms 
of the gauge-invariant fields such as X) does not contain interactions between 
hard-coUinear and soft fields. Since the external S-meson states only contain soft 
constituents, we can take the vacuum matrix element over the hard-coUinear fields, 
defining a jet function 

{Xk{tn),Xi{x + sn)} \n) = 6ki J{x + {s - t)n) + . . . , (2.5) 

where k, I are color indices, and we have used translational invariance to determine 
the dependence on the coordinate vectors. At higher orders in SCET power count- 
ing, additional jet functions would arise, but their contributions can be neglected 
at leading power. Shifting the integration variable from x to z — x+{s — t)n, with 
— X-, and introducing the Fourier-transformed Wilson coefficient functions 

Ci{n-p) = J dse^'^'-PCiis), (2.6) 

we then obtain 

^ r ^ ~ 

T^^-^iJ2 q{n ■ p) Ci{n ■ p) / e"^^'^ ^C(O) f ^ J{z) T\ J{(z_) + . . . . (2.7) 

In the next step, we rewrite the bifocal heavy-quark operator as [15, 16] 
5i{Q)VK{zJ) = {hSs){Q)Te'--^+ {Slh){Q)^h{Q)Te^--^+h{Q) 

= J duje-^'^''-^h{0)r5{uj-in- D)h{0), (2.8) 

where T may be an arbitrary (even z-dependent) Dirac structure, and we have used 
the property in- D Sg = S g in -d of the soft Wilson line Sg, where iD^ = id^ + Qs^^ 
is the covariant derivative with respect to soft gauge transformations. When this 
expression is used in (2.7), the resulting formula for the correlator involves the 
Fourier transform of the jet function, 

j d^z e-*f J{z) = Jip") , (2.9) 

however with p^ replaced by the combination p'^ = p^ + ^utn^. In defining the 
momentum-space jet function J{p^) we have taken into account that the matrix 
element in (2.5) vanishes when multiplied by i/i from either side, since yiX = 0. 
Using that p^;- — P-i we now obtain 

3 r 

T'"'^iJ2 Hij{n -p) du J{pI) h r'^j^-F^i S{u; - in ■ D) h + . . . , (2.10) 

OK 



where Hij{n • p) — C*{n ■ p) Ciin ■ p) are called the hard functions. 

In order to compute the hadronic tensor we take the discontinuity of the jet 
function, 

J{p')^-l^[iJ{p')], (2.11) 

TT 

and evaluate the S-meson matrix element of the soft operator using the HQET 
trace formalism, which allows us to write [11] 

(.|.r.(.-.„..).|.> ,,,^,i,(,i^^,,„ ,,,, 

at leading power in the heavy-quark expansion. The soft function S{uj) coincides 
with the shape function /(/c+) introduced in [15, 16]. This gives the factorization 
formula 

= X; ^^^(^ • P) Y ^) / Jipl) + ■■■■ (2-i3) 

At leading power one could replace |(_ — > and n • p — 2v ■ p^ ^ 2v ■ p in this 
result. 

In the final expressions (2.10) and (2.13) the dependence on the three scales 
n ■ p ^ rUb, pIj ~ mfeAqcD and uj ~ Aqcd has been factorized into the hard, 
jet, and shape functions, respectively. Large logarithms associated with ratios 
of these scales can be resummed by solving renormalization-group equations for 
the scale dependence of these component functions. The factorization formula 
(2.13) was derived at tree level in [15, 16, 17], and was generalized to all orders in 
perturbation theory in [18, 25]. The derivation presented above is equivalent to a 
proof of this formula presented in [22] (see also [20]). The limits of integration in 
the convolution integral are determined by the facts that the jet function defined 
in (2.11) has support for > 0, and the shape function defined in (2.12) has 
support for — oo < a; < A with A = Mb — rrib, where nib is the heavy-quark pole 
mass. The argument of the jet function can be rewritten as 

pI ^ p^ + n ■ pcu ^ n ■ p{P+ - {A - cj)) = n ■ p{P+ - cu) , (2.14) 

where the variable a) = A — a; > 0. We shall see below that expressing the convo- 
lution integral in terms of the new variable cu eliminates any spurious dependence 
of the decay spectra on the 6-quark pole mass. 

Using the fact that the Wilson coefficients Cj are real and hence Hij is sym- 
metric in its indices, we find 



E 

2j = l 



H^Jtr(T^^T^,^-±^] = n.p 



+H22 v^v" + (//12 + H23) (n^v" + n^'v^) 



(2.15) 



This result may be compared with the general Lorentz decomposition of the 
hadronic tensor: 



-g^'^'m + v^'v'^Ws + (n^'v'' + rfvi')Wi + ni'rfW^ (2.16) 

We see that the structure function W2 is not generated at leading order in the 
SCET expansion. Since only the Wilson coefficient Ci is non-zero at tree-level, the 
structure function Wi receives leading-power contributions at tree level, whereas 
W4 and W5 receive leading-power contributions at 0{as{mij)). The function W3 
receives leading-power contributions only at 0{a1{mi,)), which is beyond the ac- 
curacy of a next-to-leading order calculation. 

In the next section we will calculate the Hij from [2]. This paper uses the 
decomposition 

W>''' = (pV+pV-^''^t;-p-ie'^^"V^^^)W"i 

-g^''W2 + v^v^'W^ + {t^v'' + p''v^)Wi + p^p'^Wr,. (2.17) 

Equation (2.13) then implies that: 

= 2//n j dwJ{pl)S{u) + ... 

= 2H^2 j (kjJ{pl)S{Lj) + ... 

W, = ^Jhl [ cLjJ{pI)S{uj) + ... (2.18) 
n-p J 

2.3 Matching calculations 

In this section we derive perturbative expressions for the hard functions Hij{n • p) 
and the jet function J(p^) in (2.13) at next-to-lcading order in a^. To this end, 
we match expressions for the hadronic tensor obtained in full QCD, SCET, and 
HQET, using for simplicity on-shell external 6-quark states. We also present results 
for the renormalized shape function in the parton model, which are needed in the 
matching calculation. Throughout this chapter we use the MS subtraction scheme 
and work in d = 4 — 2e space-time dimensions. 

2.3.1 Hard functions 

Perturbative expressions at 0{as) for the structure functions Wj in the decom- 
position (2.17) have been obtained in [2] by evaluating one-loop Fcynman graphs 
for the current correlator T^^" using on-shell external quark states with residual 
momentum k (satisfying v ■ k = 0) in full QCD. The leading terms in the region of 



hard-coUinear jet momenta are 



2 = ^(pD 



' ' S\n^ y - lOlny + ^ + 4L2(1 - y) + — + 5 



47r 



+ 



1 

2 



W4 



47r 



2s CpOi-s 
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Att 1-y \ l-y 



+ (81ny-7)U 



1 + 



3 

+ ... , 



lny-1 +... , 



(2.19) 

whereas W2 and do not receive leading-power contributions at this order, in 
accordance with our general observations made above. Here as = as{lJ>), pi = 
p^+n-pn-k, and in the shape function region, where P+ <^ P_, we have y ~ n-p/mi,. 
We have used that 2v • p/nib — n ■ p/mb + 0{X) in the hard-coUinear region. The 
star distributions are generalized plus distributions defined as [2] 



ln(x/M) 



dx ■ 



Fix 



J dxF{x) 



X 



^ + F(0)lni 

X u 



F{x)-F{0) , X F(0) , 2 z 

dx — In - H — In^ - 

X u 2 u 



(2.20) 



where F{x) is a smooth test function. For later purposes, we note the useful 
identities 



A 



hi{\x/u) 
\x 



[u/\] 



1 

ln(Ax/^) ^'"/^' 
x 



+ 5{x) In A 
\q.{x/u) 



X 



(2.21) 

In order to find the hard functions Hij, we calculate the discontinuity of the 
current correlator (2.4) between on-shell heavy-quark states in momentum space. 
We work to one-loop order in SCET, keeping i,j fixed and omitting the Wilson 
coefficient functions. The corresponding tree diagram yields 



D^°^ = KS{pI), with K = Ub{v)T^^_Tl 



Ub[V) 



(2.22) 



where ui){v) are on-shell HQET spinors normalized to unity, and the quantity 
K corresponds to the Dirac trace in (2.13). The interpretation of this result in 
terms of hard, jet, and soft functions is that, at tree level, J^^^pl,) = 5(p^) and 
'S'parton('^) = ^{'^ — 71 ■ k). (The secoud result is specific to the free-quark decay 
picture.) Then the convolution integral / du;J{p1)S{u;) in (2.13) produces 5(p|), 
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Figure 2.1: One-loop diagrams contributing to the current correlator in SCET. 
The effective current operators are denoted by crossed circles, and hard- 
coUinear propagators are drawn as dashed lines. Mirror graphs obtained by 
exchanging the two currents are not shown. 



and comparison with (2.19) shows that = 1, while all other hard functions 
vanish at tree level. 

The diagrams contributing at one- loop order are shown in Figure 2.1. They are 
evaluated using the Feynman rules of SCET. The first three graphs contain hard- 
collinear gluon exchanges, while the last two diagrams contain soft exchanges. 
The wave-function renormalization factors of the external heavy quarks equal 1 
on-shell. For the sum of all hard-coUinear exchange graphs, we find 
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(2.23) 
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The sum of the soft contributions is given by 



Di'^ = K 



An 



~--L+^--AL^ + AL-^]5{pi) 
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- + 8L 
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(2.24) 



where L = \n{n -p/fi). The 1/e poles in the sum of the hard-collinear and soft 
contributions are subtracted by a multiplicative renormalization factor Zj applied 
to the bare current correlator in (2.4), where 



1 + 
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(2.25) 



is the (momentum-space) current renormalization constant in SCET [21]. Taking 
the sum of (2.23) and (2.24) after subtraction of the pole terms, and matching it 



on 



with the results in (2.19), we find that at one-loop order 
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(2.26) 



In deriving these results we have used the identities (2.21) to rearrange the var- 
ious star distributions. The remaining hard functions start at 0{al). Using the 
relation Hij — CiCj, one can derive from these results expressions for the Wilson 
coefficients in the expansion of the semileptonic current in (2.2). We confirm the 
expressions for these coefficients given in [21]. 

2.3.2 Jet function 

After the hadronic tensor is matched onto HQET as shown in (2.13), the SCET 
loop graphs in Figure 2.1 determine the one- loop contributions to the product of 
the jet function and the shape function in (2.13). We may write this product 
symbolically as J'^^^ ® S^^^ -\- J^^^ S'^^\ where the (8) symbol means a convolution 
in CO. Whereas the jet function is a short-distance object that can be calculated in 
perturbation theory, the shape function is defined in terms of a hadronic matrix 
element and cannot be properly described using Feynman diagrams with on-shell 
external quark states. In a second step, we must therefore extract from the results 
(2.23) and (2.24) the one-loop contribution J^^^ to the jet function. To this end, 
we must compute the renormalized shape function at one-loop order in the parton 
model. This will be done in the next subsection. We may, however, already 
anticipate the result for the jet function at this point, because at one-loop order 
the graphs in Figure 2.1 can be separated into diagrams with hard-coUinear (first 
line) or soft (second line) gluon exchange. (This separation would be non-trivial 
beyond one-loop order.) The hard-coUinear contribution in (2.23) thus determines 
the convolution J^^^ S'^'^^ = J^^\pl), whereas the soft contribution in (2.24) 
corresponds to J*^°) <^ S^^\ It follows that the renormalized jet function is given by 
the distribution 



This result disagrees with a corresponding expression obtained in [35]. It will often 
be useful to separate the dependence on n ■ p and P+ in this result by means of 
the substitution p"^ — yp^, where = rnb{P+ — u) according to (2.14). Using the 



ApI) = 5{pI) + 
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Figure 2.2: Radiative corrections to the shape function. The bilocal HQET 
operator is denoted by the black square. A mirror copy of the first graph is 
not shown. 



identities (2.21), we find 



+ 4 
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+ (41n|/-3) - 



(2.28) 



The jet function is non-zero only if y > and P+ > cD, which ensures that pI > 0. 



2.3.3 Renormalized shape function 

Having calculated the short-distance objects Hij and J in the factorization formula 
(2.13) at one-loop order, we now turn to a study of radiative corrections to the 
shape function S{uj). There is considerable confusion in the literature about the 
renormalization properties of the shape function, and several incorrect results for 
its anomalous dimension have been published. We therefore present our calculation 
in some detail in this subsection and the following section. 

According to (2.12), the shape function is defined in terms of a hadronic ma- 
trix element in HQET and thus cannot be calculated perturbatively. However, 
the renormalization properties of this function can be studied using perturbation 
theory. To this end, we evaluate the matrix element (2.12) in HQET using exter- 
nal heavy-quark states with residual momentum k. For the time being, we keep 
V ■ k non-zero to regularize infra-red singularities. The relevant one-loop graphs 
are depicted in Figure 2.2. Adding the tree contribution, we obtain for the matrix 
element of the bare shape-function operator 0{u) = JiT S{u! — in-D)h (expressed 
in terms of renormalized fields) 

= Z,(5(u;-n-A;)-i^r(l + e) 



X 



(47r) 

5{uj -n-k) - 5{uj -n-k + Z)] ^1 + 
+ e{n-k-u){n-k- uj)-'{n ■ k - u + S)-^-' \ , (2.29) 




where 5 — —2v ■ k, and 



Zh-i + r(2e) r(i - e) r^^ 



(47r) 



(2.30) 



is the off-shell wave-function renormalization constant of a heavy quark in HQET. 
The next step is to extract the ultra-violet poles from this result, which determine 
the anomalous dimension of the shape function. We define a renormalization factor 
through 



'S'ren('^) 



duj' Zs{lJ,Uj') 5'bare(t^') , 
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(2.31) 



zs{u;,u;') + ... . 



The result for Zs following from (2.29) must be interpreted as a distribution on 
test functions F{u') with support on the interval — oo < u' < A. We obtain 
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Note the peculiar dependence on the parameter A setting the upper limit on the 
integration over oo' in (2.31), which combines with the plus distribution to form a 
star distribution in the variable {uo' — uj). 

We can now determine the renormalized shape function from (2.31). The result 
must once again be interpreted as a distribution, this time on test functions F{uj) 
integrated over a finite interval — Ahad < OJ <A. In practice, the value of Ahad is 
set by kinematics or by virtue of some experimental cut. The result is 
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While it was useful to keep the heavy quark off-shell in the calculation of the ultra- 
violet renormalization factor, the limit 6 = —2v • A; — > can be taken in the result 
for the renormalized shape functions without leading to infra-red singularities. 
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This gives 
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where the star distributions must now be understood as distributions in the variable 
{n ■ k — u). 

We stress that these results for the renormalized shape function are obtained in 
the parton model and can in no way provide a realistic prediction for the functional 
form of S{u;). This should be obvious from the fact that our results depend on a 
single "hadronic parameter" n ■ k, corresponding to a fixed residual momentum of 
the heavy quark. Only the dependence on the ultra-violet renormalization scale 
11 can be trusted. However, the one-loop result in (2.34) is needed to complete 
the matching calculation of the jet function described in the previous subsection, 
which can legitimately be performed with on-shell external 6-quark states. Given 
the expression for the renormalized shape function, we obtain 
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With the help of the identities (2.21) this can be shown to be equal to the finite 
part of (2.24), as we claimed above. 



2.4 Renormalization-group resummation 

Equations (2.26) and (2.27) determine the short-distance objects Hij and J in the 
factorization formula (2.13) at one-loop order in perturbation theory. However, 
there is no common choice of the renormalization scale that would eliminate 
all large logarithms from these results. Likewise, the shape function, being a 
hadronic matrix element, is naturally renormalized at some low scale, whereas 
the short-distance objects contain physics at higher scales. The problem of large 
logarithms arising from the presence of disparate mass scales can be dealt with 
using renormalization-group equations. Proceeding in three steps, our strategy 
will be as follows: 

i) At a high scale jih ~ f^h we match QCD onto SCET and extract matching 
conditions for the hard functions Hij. The corresponding one-loop expressions 
have been given in (2.26). At that scale, they are free of large logarithms and 
so can be reliably computed using perturbation theory. We then evolve the hard 
functions down to an intermediate hard-coUinear scale /Xj ~ a/^«6Aqcd by solving 
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the renormalization-group equation 



Hij{n • p, //) = 27j(n • p, //) Hij{n • p, //) , (2.36) 

where 7j is the anomalous dimension of the semileptonic current in SCET. 

ii) Next, we start from a model for the shape function S{u!, iiq) at some low scale 
//o = few X Aqcd large enough to trust perturbation theory. Such a model could 
be provided by a QCD-inspired approach such as QCD sum rules or lattice QCD, 
or it could be tuned to experimental data. We then solve the integro-differential 
evolution equation 
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dluji 



S{ijj, 1^^ ~ ~ J 'ysi'-^, ■, i^) S{ijj\ 11) (2.37) 



to obtain the shape function at the intermediate scale /Xj. 

iii) Finally, at the scale /ij we combine the results for the hard functions and 
for the shape function with the jet function J in (2.27), which at that scale is free 
of large logarithms and so has a reliable perturbative expansion. The dependence 
on the matching scales //ft and //j cancels in the final result (to the order at which 
we are working). 

We now discuss these three steps in detail. 



2.4.1 Evolution of the hard functions 

At one-loop order, the anomalous dimension 7j for the SCET current is twice 
the coefficient of the 1/e pole in the renormalization factor Zj in (2.25). More 
generally [21, 26], 

'^j{n-p,\i) = -rcusp(as) InJ^ + 'f^as) = — ^ ( -In^^^^ 7 1 + • • • > (2-38) 

n • p \ n ■ p 4 J 

where Fcusp = Cfois/tt + ■ ■ ■ is the universal cusp anomalous dimension governing 
the ultra-violet singularities of Wilson lines with hght-like segments [27, 28]. The 
exact solution to the evolution equation (2.36) can be written as 

Hij{n ■ p, Hi) = Hij{n ■ p, Hh) Uy{iJ,h, l^i) , (2.39) 

where 

Uyifih, f^i) = UifXh, f^i) y-^^rC/^..^^) . (2.40) 
The exact expression for the evolution factor reads 

lnU{iih,l^i) = 25'r(/Xft,Aii) - 2ar(Aift,Aii) In 2ay (///,,//,) , (2.41) 
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where the functions of the right-hand side are solutions to the renormahzation- 
group equations 

^j^5'r(i^,Ai) = -rcusp(Q;s(A*)) In^ , 

^^ar(z^,/i) = -rcusp(a.(/i)) , ^^ay(z/, /x) = -7'(a,(/i)) , (2.42) 

with boundary conditions /i) — etc. at /i — u. These equations can be 
integrated using that d/dln/i — (5{as) d/dag- The solutions are 

ois(lJ) a OisilA 

Tcusp (Q?) 

(2.43) 

and similarly for a^i. Explicit results for the Sudakov exponent Sy and the func- 
tions Or and a,^, in (2.41) at next-to-leading order in renormalization-group im- 
proved perturbation theory are given in Appendix B. 

2.4.2 Evolution of the shape function 

At one-loop order, the anomalous dimension for the shape function is twice the 
coefficient of the 1/e pole in the renormalization factor Zs- From (2.32), we obtain 
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The evolution equation (2.37) can be solved analytically using a general method 
developed in [29]. It is convenient to change variables from u tou = A — u G [0, cxd[ 
and denote S{cu) = S{A — cD). The renormalization-group equation then reads 
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where the anomalous dimension can be written in the general form 
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(2.46) 

The logarithmic term containing the cusp anomalous dimension has a geometric 
origin. Since the heavy-quark field h{x) in HQET can be represented as a Wilson 
line along the v direction, the field J{(a;) entering the SCET formalism contains 
the product of a light-like Wilson line (along n) and a time-like Wilson line (along 
v), which form a cusp at point x. The shape function contains two such cusps. 
According to the renormalization theory of Wilson lines with light-like segments, 



each cusp produces a contribution to the anomalous dimension proportional to 
Feusp In ^[27, 28]. 

The general solution of (2.45) can be obtained using the fact that 

- J^^dw' {cu'r reusp(«.) (%^f|7^) = reusp(«.) [pil + a) + 7e] (2.47) 

where ijj{z) is the logarithmic derivative of the Euler F function and the integral 
on the left-hand side is convergent as long as Re a > —1. 
Relation (2.47) imphes that the ansatz [29] 

f{u, /io, r) = I — ) U2in, l^oWsir, //q) (2.48) 



with 



In 1/2(11, iiq) = 2S{iio, IJ) + 2a^(//o, jj) , 



/dcx 
rcusp(Q;s) |V' [1 + T - 2ar(//o, A«a)] + 7fi} , 



(2.49) 
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provides a solution to the evolution equation (2.45) with initial condition, 
/(a), yUo, /^o, t) = {Oj/ ixqY at some scale //q- Here is defined such that as{Ha) = 
a-y is defined as in equation (2.43) with Fcusp replaced by 7, and r can be an 
arbitrary complex parameter. Note that ar(/i, /^o) < if > /xq. 

We now assume that the shape function S{Co, /io) is given at the low scale jiQ 
and define its Fourier transform with respect to ln(a;//xo) through 

S{uj,i,o) = ^ r dtSo{t)(-\ . (2.50) 

The exact result for the shape function at a different scale ji is then given by 

1 f°° 

S{Lb, //) = —/ dt So{t) fiu, 1^0, it) . (2.51) 

With the help of this formula, it is straightforward to derive explicit expres- 
sions for the evolution of the shape function from the hadronic scale up to the 
intermediate scale /Xj at any order in renormalization-group improved perturbation 
theory: 

/^ \it-2ar{no,Hi) r(l -\-if'\ 

fico, /.o, ^^) = t/2(/x. /xo)e^-('^-'^^)^- - ' + 
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(2.52) 



The leading-order result presented above can be simplified further. When (2.52) 
is inserted into (2.51), the integration over t can be performed analytically. Setting 
r] = —2ar{iJ>o,iJ>i), the relevant integral is 
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where 

(2.54) 



Mt)^ / —S{u;',fXo)i — ) 
Jo ^ \/^o/ 



is the Fourier transform of the shape function as defined in (2.50). The integrand 
of the t-integral has poles on the positive imaginary axis located at t = in with 
n > 1 an integer. For uj < uj' the integration contour can be closed in the lower 
half-plane avoiding all poles, hence yielding zero. For cu > uj' we use the theorem 
of residues to obtain 



7 = / dw'R{u,u')S{u;',iio), (2.55) 
Jo 



where 



l°°/cj'\-' 1 1 1 

Note that Rioj^uj') —>■ S{u; — uj') in the limit — > 0, corresponding to /iq, as it 
should be. Our final result for the shape function at the intermediate hard-coUinear 
scale, valid at leading order in rcnormalization-group improved perturbation the- 
ory, can now be written in the simple form (valid for //j > /xq, so that r] > 0) 

Siu, /..) = U,{^^, ^o) I du ■ (2.57) 

From the above equation one can derive scaling relations for the asymptotic 
behavior of the shape function for a) — > and a) — > oo (corresponding to a; — > A 
and u —oo). If the function S'((u),/io) at the low scale /Xq vanishes proportional 
to uj^ near the endpoint, the shape function at a higher scale > n vanishes faster, 
proportional to (2;^+''. Similarly, if 5 (a), /xq) falls off like (2;~^ for (2; — > oo, the shape 
function renormalized at a higher scale vanishes like u;~^™('^'^)+'n. Irrespective of 
the initial behavior of the shape function, evolution effects generate a radiative tail 
that falls off slower than 1/a). This fact implies that the normalization integral 
of S{uj,fi) as well as all positive moments are ultra-violet divergent. The field- 
theoretic reason is that the bilocal shape-function operator in (2.8) contains ultra- 
violet singularities as — > 0, which are not subtracted in the renormalization of 
the shape function. The situation is analogous to the case of the S-meson hght- 
cone distribution amplitude discussed in [29, 30]. These divergences are never an 
obstacle in practice. Convolution integrals with the shape function are always cut 
off at some finite value of a) by virtue of phase-space or some experimental cut. 
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2.5 Properties of the shape function 



In this section we discuss how moments of the shape function are related with 
HQET parameters. This will lead us to propose a new, physical scheme for defining 
a running heavy-quark mass, which is most appropriate for the study of inclusive 
spectra in the shape-function region. We will also present a model-independent 
result for the asymptotic behavior of the renormalized shape function (defined in 
the MS scheme), finding that it is not positive definite. 

Most of our discussion in this section is phrased in terms of the original (un- 
hatted) shape function S{u!, /i). At the end, we formulate the resulting constraints 



2.5.1 Shape-function moments in the pole scheme 

Naively, ignoring renormalization effects, the moments Mj^ = duJUJ^S{uj) are 
given by hadronic parameters defined in terms of 5-meson matrix elements of 
local HQET operators [15]. In particular, Mq = 1 fixes the normalization of the 
shape function, Mi = vanishes by the HQET equation of motion, and M2 = 
— Ai/3 is determined by the matrix element of the kinetic-energy operator. The 
vanishing of the first moment is connected with the implicit definition of the heavy- 
quark pole mass built into the HQET Lagrangian via the equation of motion 
iv ■ D h = 0. These moment constraints have been implemented in various model 
paramcterizations for the shape function suggested in the literature [31, 32, 33]. 
Typically, one makes an ansatz for the shape function depending on a few HQET 
parameters such as A and Ai, and determines the values of these parameters from 
a fit to experimental data. 

Beyond tree level, all moments Mn with > receive ultra-violet divergences 
from the region uj —00 (or uj 00). However, as we have mentioned above, the 
values of uj needed for the description of physical decay rates are always restricted 
to a finite interval. It is thus sufficient for all purposes to define the moments of 
the renormalized shape function as 



The dependence of these moments on the renormalization scale // is controlled 
by the evolution equation (2.37). In addition, the moments depend on the lower 
cutoff on the u integral. The choice of Auv is a matter of convenience, and so we 
are free to pick a value that is numerically (if not parametrically) large compared 
with Aqcd- In this case, as we will now show, the dependence on Auv can also be 
controlled using short- distance methods. 

For sufficiently large values of Auv it is possible to expand the moments 
Mat (Auv, A*) in a series of 5-meson matrix elements of local HQET operators. If 
for simplicity we set F = 1 in the shape-function operator (which is legitimate. 



on the function S{lj, /i). 




(2.58) 
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since the Dirac structure is unaltered in HQET), the operators in question are 
Lorentz-scalar, "leading- twist" operators containing h . . .h [15, 17]. These are the 
operators that mix with h (in ■ D)^h under renormahzation. It is straightforward 
to find the corresponding operators of a given dimension. The unique dimension- 
3 operator is hh. The two operators of dimension 4 are hin • Dh (class-1) and 
hiv-Dh (class-2). The class-2 operator vanishes by the HQET equation of motion. 
The possible dimension-5 operators are 



class-1: h{in- D)'^h, h{iD±yh, 

class-2: h{iv-Dyh, hin-Div-Dh, hiv-Din-Dh, 



where again the class-2 operators vanish by the equation of motion. Moreover, 
it follows from the Feynman rules of HQET that the two class-1 operators do 
not mix under renormahzation, so the operator h {iD±yh can be ignored. Prom 
dimension 6 on the situation is more complicated, because several class-1 operators 
exist that can mix with h {in ■ D)^h. For N = 3 these arc of the form hiDGh 
or h . . . q q . . . h, where we omit Lorentz and color indices. We will restrict our 
discussion to operators of dimension less than 6. 

For the operator product expansion of the moments in (2.58) we need the 
forward matrix elements 

<o)^m^ ,.eo) 

of the leading-twist operators between i?-meson states in HQET. Using the equa- 
tion of motion, it can be shown that {hh) = 1, {hin-Dh) = 0, and {h {in-D)'^h) — 
— Ai/3 [15]. We can thus write an expansion of the form 



M;v(Auv, 1^) = Ajiv { (Auv, fi) + i^f^ (Auv, fi) ■ + O 
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Auv 

(2.61) 

This expansion is useful as long as the cutoff Auv is chosen much larger than the 
typical hadronic scale characterizing the matrix elements of the local operators. 
The matching coefficients A'n^'* in this relation can be calculated using on-shell 
external 6-quark states with residual momentum k. For operators of dimension up 
to 5 it suffices to calculate two-point functions. (Three and four-point functions 
would have to be considered at dimension 6.) We first evaluate the moments of 
the renormalized shape function in (2.34), finding at one-loop order 

Mr'°'^(Auv,/^) (2.62) 
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We then expand this result in powers of n ■ k/Auy. Keeping the first three terms 
in the expansion, we obtain 
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In the next step, we calculate the one-loop matrix elements of the local operators 
h {in • D)^h between heavy-quark states with residual momentum k. The relevant 
diagrams are the same as in Figure 2.2, where now the black square represents 
the local operators. Keeping v ■ k non-zero to regularize infra-red singularities, we 
obtain for the bare matrix elements 



{h{in-DYh) = {n-kyll 



^ / AT \ 
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(2.64) 



While the individual diagrams are infra-red divergent, taking the limit v • A; — > in 
the sum of all contributions is possible without encountering singularities. Then 
the one-loop contributions vanish, and the matrix elements simply reduce to their 
tree-level values. In other words, the one-loop contributions correspond to a mixing 
with class-2 operators, whose hadronic matrix elements vanish by the equations 
of motions. It follows that in (2.63) we must identify (n ■ A;)" {h{in ■ D)^h). 
Substituting the results for the HQET matrix elements given earlier, we obtain for 
the Wilson coefficients of the first three moments 
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At tree level, this reproduces the naive moment relations mentioned at the begin- 
ning of this section. Beyond tree level, the moments get corrected by calculable 
short-distance effects, which can be controlled using fixed-order perturbation the- 
ory as long as the ratio Auv/yU. is of 0(f). In particular, the renormalized first 
moment no longer vanishes, but is proportional to the cutoff Auv up to small 
power corrections. 

As mentioned earlier, the value of the first moment is connected with the defi- 
nition of the heavy-quark mass (see also [34, 35]). The first moment of the renor- 
malized shape function can be made to vanish to all orders in perturbation theory 
by choosing an appropriate scheme for the definition of m^. So far our calculations 
have assumed the definition of the heavy-quark mass as a pole mass, which 
is implied by the HQET equation of motion iv ■ D h = 0. Results such as (2.65) 
are valid in this particular scheme. A more general choice is to allow for a residual 
mass term Sm in HQET, such that iv ■ D h = Smh with 6m = O(Aqcd) [36]. 
It is well known that the pole mass is an ill-defined concept, which suffers from 
infra-red renormalon ambiguities [37, 38]. The parameter Apoie = Mb — ^1°^^, 
which determines the support of the shape function in the pole-mass scheme, in- 
herits the same ambiguities. It is therefore advantageous to eliminate the pole 
mass in favor of some short- distance mass. For the analysis of inclusive 5-meson 
decays, a proper choice is to use a so-called low-scale subtracted heavy-quark mass 
rribdif) [39], which is obtained from the pole mass by removing a long-distance 
contribution proportional to a subtraction scale /if — few x Aqcd, 



As long as mb{(JLf) is defined in a physical way, the resulting perturbative expres- 
sions after elimination of the pole mass are well-behaved and not plagued by renor- 
malon ambiguities. Replacing the pole mass by the physical mass shifts the values 
oin - k and lv by an amount 6m, since n ■ {m'^°^'^v + k) = mi,{iJ,f) + {n - k + 6m), and 
because the covariant derivative in the definition of the shape function in (2.12) 
must be replaced by in ■ D — 6m [36]. At the same time, Apoie = -^(a*/) ~ 
where A(/i/) = Mb — mb{fif) is a physical parameter. Note that this leaves the 
parameter lj = A — uj and hence the shape function S{i2i, fi) invariant. This follows 
since Apoie — tfpoie = A{nf) — (wpoie + 6m), where utpoie denotes the value in the 
pole-mass scheme used so far. 

2.5.2 Shape-function moments in a physical scheme 

From now on we will adopt a mass scheme defined by some specific choice of 6m. 
Let us denote by a; = o^poie + 6m the value of the light-cone momentum variable 
in that scheme and define "physical" moments Mf^^^ as in (2.58), but with all 
parameters replaced by their values in the new scheme, in particular A = A{fif). 
Then the expressions for the moments in (2.62) and (2.63) change according to the 
replacements n- k —>■ n- k + 6m and Ayy —>■ Auv ~ Sm everywhere. We now choose 




(2.66) 
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5m such that the first moment vanishes, thereby defining a low-scale subtracted 
heavy-quark mass (with fJ^f = Ayy) to all orders in perturbation theory. We will 
refer to this mass as the "shape-function mass" mf^. This is a "physical", short- 
distance mass in the sense that it is free of renormalon ambiguities. (However, the 
definition of the shape-function mass depends on the renormalization scheme used 
to define the shape function.) Prom (2.63) and (2.66), it follows that at one-loop 
order 
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Note that after introduction of the shape-function mass the coefficients Kn^ in 
the operator-product expansion for the moments in (2.61) vanish by definition. 
However, to first order in the values for the coefficients Kn^ and Kn^ of the 
zeroth and second moments given in (2.65) remain unchanged, since 5m — 0(0;^). 
This would no longer be true for the coefficients of higher moments. 

The shape-function mass can be related to any other short-distance mass using 
perturbation theory. For instance, at one-loop order its relations to the potential- 
subtracted mass introduced in [40] and to the kinetic mass defined in [41, 42] read 
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Note that, in addition to the dependence on the subtraction scale fif, the shape- 
function mass depends on the scale /j, at which the shape function is renormalized. 
While it is natural to set = as we did here, this is not necessary. Given a 
value for the shape-function mass for some choice of scales, we can solve (2.67) 
to obtain its value for any other choice, using the fact that the pole mass is scale 
independent. 

Proceeding in an analogous way, we can use the second moment to define a 
physical kinetic-energy parameter, commonly called //^. This quantity can be 
used to replace the HQET parameter Ai, which like the pole mass suffers from 
infra-red renormalon ambiguities [43, 44]. At one-loop order, we obtain 
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Taking the ratio of M^"^^^ and M^'^^^ has the advantage of efiminating the double 
logarithmic radiative corrections from this expression. Our definition is similar to 
the running parameter /x^ defined in the kinetic scheme [41, 42]. At one- loop order, 
the two parameters are related by 



/^^(/^/,/^/) = -/^/ ^^^'}^^^ + [/i^(/i/)]kin 
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Given a value for the kinetic energy in the shape-function scheme for some choice 
of scales, wc can solve (2.69) to obtain its value for any other choice, using that Ai 
is scale independent. 

Similarly, each new moment of the renormalized shape function can be used to 
define a new physical, scale-dependent parameter 

^^(Auv,/^) ^ MP'^^^(Auv,/^)/MP^"^(Auv,A^), 

which coincides with the corresponding HQET parameter Aj^ — (h {in ■ D)^h) at 
tree level, and which beyond tree level is related to HQET parameters through well- 
controlled perturbative expressions. Obviously, the presence of power divergences 
implies that higher moments are progressively less sensitive to HQET parameters, 
since they are dominated by the perturbative terms of order ct^Auv 



2.5.3 Moments of the scheme-independent function S'(a>, /x) 

It will be useful to rewrite the moment relations derived above in terms of the 
variable uj — A — a;, which is invariant under redefinitions of the heavy-quark 
mass. Defining a new set of scheme-independent moments 
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we obtain 
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where the parameters A(/x/, /x) = MB—mf^{iJ,f, fi) and l^i^fJ^fi A*) should be consid- 
ered as known physical quantities. Using the relations in the previous subsection, 
we have 



CFO^s{^J^*) 



TT 



TT 

21n^U^^%^ln^ 



1^ 



27r 



TT 



31n^ + 



3^2 CpOisifJ') 1^ f^f 



TT 



In^:^, 



(2.73) 



where denotes the scale at which initial values for the two parameters are 
obtained, for instance using relations such as (2.68) and (2.70). These relations 
are particularly simple if one chooses fif = fi. 

In the relations above we have eliminated the unphysical HQET parameter 
Ai in favor of the physical parameter /x^ defined in the shape-function scheme. 
At first sight, this seems to threaten the convergence of the operator product 
expansion. For instance, the term proportional to fi^ in the expression for the 
zeroth moment Mq in (2.72) contains a leading-power perturbative contribution 
of order a^(/i.), and similar contributions would arise from all other terms in the 
expansion. These contributions would have to be subtracted from the Wilson 
coefficient of the first term, if this coefficient were computed to two-loop order. 
The overall convergence of the operator product expansion is unaffected by this 
reorganization of perturbative corrections. 



2.5.4 Asymptotic behavior of the shape function 



The fact that for sufficiently large values of the cutoff the moments of the shape 
function can be calculated using an operator-product expansion implies that a 
similar expansion can be used to obtain a model-independent description of the 
asymptotic behavior of the shape function. Taking the derivative of the zeroth 
moment Mq in (2.71) with respect to ///, one obtains 
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This relation can be trusted as long as /ij ^ Aqcd- It allows us to determine 
the behavior of the shape function for large values of u. From (2.72) we find at 
one-loop order 
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The precise definitions of A and /^^ are not specified at this order. (Note that 
the shape function cannot depend on the value of the cutoff ///.) We have checked 
that this asymptotic behavior of the shape function is consistent with the evolution 
equation (2.57) when expanded to first order in 

Relation (2.75) is a model-independent result as long as ^ Aqcd- Wc stress 
the remarkable fact that this radiative tail of the shape function is negative, in 
contrast with the naive expectation based on a probabilistic interpretation of the 
shape function as a momentum distribution function. The point is that the def- 
inition of the renormalized shape function requires scheme-dependent ultra-violet 
subtractions. From (2.75) it follows that the shape function must have a zero, 
which for sufficiently large is located at a value luq ^ A-\- \xl \fe- 



A A 



2.6 Conclusions 

The hadronic physics governing the inclusive semileptonic decay B — > Xy^l'u is 
encoded in the structure functions Wi appearing in the Lorentz decomposition 
of the hadronic tensor W^^ in (2.16). In the shape-function region, only two 
combinations of these functions are required at leading order in Aqcd/'"^6- They 
follow from the factorization formula (2.13) using the explicit results for the hard 
functions Hij and the jet function J derived in this chapter. Explicitly, we obtain 
at next-to-leading order in renormalization-group improved perturbation theory 



where the dots represent power corrections in Aqcd/'"^6- 

In all our results, mf, denotes the heavy-quark mass defined in the "shape- 
function scheme" introduced in Section 2.5. The jet function J(|/m(,(P+ —Cj),iii) 
at an intermediate hard-coUinear scale //j ~ \A^^feAQCD can be calculated in fixed- 
order perturbation theory. The relevant expression valid at one-loop order is given 
in (2.27). The above results contain the renormalization-group function Uy{fih, fJ'i), 
which arises in the solution of evolution equations discussed in Section 2.4. Our 
results are formally independent of the precise choice of high-energy matching 
scale ///i ~ 777.6 . The numerical effect of the residual /ih dependence remaining after 
truncation of the perturbative expansion has been studied in [26] and was found 
to be small. 

The function S{uj,Hi) in (2.76) is the shape function after the transformation 
of variables from u to lj = A — u. The limits of integration for the variable a) 
(i.e., < a) < P+) are set by hadronic kinematics and are independent of the 
definition of the heavy-quark mass. The shape function is a non-perturbative 
object, which at present cannot be predicted from first principles. It enters our 
results (2.76) rcnormalizcd at the intermediate hard-collinear scale yUj. In (2.57), 
we have presented an analytic formula (valid at all orders in renormalization-group 
improved perturbation theory) that relates the shape function at a high scale to 
the shape function renormalized at a low hadronic scale. Many properties of the 
shape function that were so far unknown have been derived in Section 2.5. In 
particular, we have given explicit formulae relating the moments of the shape 
function to HQET parameters, and we have proved that the shape function has 
a negative tail for large values of a), whose explicit form can be calculated using 
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an operator product expansion. These new insights about the shape function will 
be very helpful in constructing a realistic model for the function S{uj,fii), which 
can then be refined by tuning it to experimental data such as the photon energy 
spectrum in inclusive 13 — > Xgj decays. 



Chapter 3 

Non-Perturbative Corrections 

3.1 Introduction 



In the previous chapter we have seen how to incorporate leading order perturbative 
corrections to i? — > Iv- Corrections suppressed by a power of Aqcd/'"^6 arc con- 
sidered the second largest source of uncertainty, for inclusive charmless B decays. 
The present chapter is devoted to a more thorough study of power corrections to 
inclusive B decays distributions in the shape-function region, using the two-step 
matching procedure. 

In a first step, hard fluctuations are integrated out by matching the QCD cur- 
rents onto soft-coUinear effective theory (SCET) [21, 23]. The current correlator 
is then expanded in terms of light-cone operators in heavy-quark effective theory 
(HQET) [11], thereby integrating out fluctuations at the hard-coUinear scale. The 
fact that the relevant HQET operators hve on the hght cone follows from the struc- 
ture of the multipolc expansion of soft fields in SCET. We carry out the matching 
procedure at tree level and to order Aqcd/'"^6 in the heavy-quark expansion. We 
indicate how our results would change if loop corrections were included. 

3.2 Short-distance expansion of the hadronic tensor 

SCET is the appropriate effective field theory for the description of the interactions 
among soft and hard-collinear degrees of freedom. Its Lagrangian is organized in 
an expansion in powers of \/A. The leading-order Lagrangian is 

^SCET = f ? f • + gn- As {x_ ) + ip±hc ipLhc] i 

z \ m ■ Ufic J 

+ qipsq + hivDsh + Cfy^, (3.1) 

where ^ is a hard-collinear quark field, g is a soft, massless quark field, /i is a heavy- 
quark field defined in HQET, Ag is a soft gluon field, and iZ)^^ = id'^ -\- gA^^ is the 
covariant derivative containing a hard-collinear gluon field. All fields in the above 
Lagrangian are evaluated at point except for the soft gluon field in the first 
term, which is evaluated at a;_ = ^ (n ■ x) n. The explicit form of the leading-order 
Yang- Mills Lagrangian can be found in [23, 24]. 

The terms up to second order in the expansion in -\/A are 

-^SCET — 'i-'^ -r l-'^q '-'h '^YM ' \'^-^) 

where 

1 r_ , r. 1 

(3.3) 
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is the next-to- leading term in the expansion of the HQET Lagrangian [11]. Ex- 
pressions for the remaining Lagrangian corrections have been presented in [23]. 

While it is consistent to apply a perturbative expansion at the hard and 
hard-coUinear scales, the soft-gluon couplings to hard-coUinear fields are non- 
perturbative and must be treated to all orders in the coupling constant. For 
instance, the hard-coUinear quark propagator derived from (3.1) should be taken 
to be the propagator in the background of the soft gluon field, summing up arbi- 
trarily many insertions of the field Ag. The most convenient way of achieving this 
summation is to decouple the leading-order interactions between the soft gluon 
field and hard-coUinear fields in (3.1) with the help of a field redefinition, under 
which [21, 45] 



ax) = S{x.) ^('\x) , Alix) = S{x.) 4f (a:) S^x.) , (3.4) 



where 



^ —oo ' 

is a soft Wilson line along the n direction. Introducing the new fields into the 
Lagrangian yields 

Lf = I L . + i^it , (3.6) 

and similarly all interactions between soft and hard-coUinear gluon fields are re- 
moved from the Yang-Mills Lagrangian Cyu- The propagator of the new hard- 
coUinear quark field is now given by the simple expression 

A,ix-y) = (0|T{e(°)(x)eW(y)}|0) = y (01 (3.7) 

The effect of soft-gluon attachments is taken into account by factors of the Wilson 
line S in the results below. 

We now list the expressions for the subleading corrections to the SCET La- 
grangian in terms of the redefined fields, using the formalism of gauge-invariant 
building blocks [46]. We define 

X = lyt^(o) ^ j^i^^^ ^ W^iiD'^^^W) , (3.8) 

where 

M. = Pexp(./«..£)(..«.)) (3.) 

^ — OO ' 

is a hard-collinear Wilson line. These "calligraphic" fields are invariant under both 
hard-coUinear and soft gauge transformations. Note that in the light-cone gauge. 



AS 



n ■ A^^l — 0, we simply have X = and A'^^ — gA^^l'^. In terms of these fields, 
the results compiled in [23] take the form 



+ ^YAS'9G,.S)^_-r^^1p^,)jX, 

4? = m,j'V>±hcX + h.c., (3.10) 

where iD'^^ — id'^ + A'^^, and we have dropped the subscript "s" on the soft co- 
variant derivative and field strength. The expression for C^^^ will not be needed 
for our analysis. The notation (. . . )x_ indicates that, in interactions with hard- 
collinear fields, soft fields are multipole expanded and live at position whereas 
hard-coUinear fields are always evaluated at position x. Because the SCET La- 
grangian is not renormalized [23], the above expressions are valid to all orders in 
perturbation theory. 

Next, we need the expressions for heavy-light current operators in SCET. In 
general, a QCD current qTb matches onto 

q(x) r b(x) = e-'^"^- + + + 4^ + J^^ + ...), (3.11) 

where we distinguish between type-74 "two-particle" operators and type-i? "three- 
particle" operators [47] . The operators arising at tree level are [23, 47, 48] 

j« = xrx^^{s^D,h)^ +xL*^^^^r{s^h)^ , 
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+ X^i%^^Txl{S^D^h)^ , (3.12) 
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+ S%ii^|r^4x..(5*/.). . (3.13) 

iii'b tn ■ o ^ ^ 

The expressions for the currents beyond tree level are more complicated, primarily 
because several new Dirac structures appear. The relevant formulae are known 
at leading [21] and next-to- leading order [47] in the power expansion. The corre- 
sponding results for the currents J\ and have not yet been derived. They 
would be needed if the analysis in this chapter should be extended beyond tree 
level. 

If perturbative corrections at the hard-coUinear scale are neglected, the hard- 
coUinear gluon fields can be dropped, and the above expressions for the effective 
Lagrangians and currents simplify. In this approximation X — > — > 0, 
iD^^ id^, and 0. While this leads to great simplifications in the calcu- 

lation, we stress that the structures of soft fields that arise do not simplify. We 
find operators containing S^gG^i^S, S'^lDp, gGi^i,]S, S'^h, S'^D^h, and S'^D^D^.h, 
and the same operators would arise if the calculation was extended beyond the 
tree approximation. The only exception is that we no longer retain Lagrangian 
corrections containing the soft quark field, because C^^^ — > in the limit where the 
hard-collinear gluon field is neglected. (The term qS i^±^^^^ is forbidden by mo- 
mentum conservation.) In Section 3.5, we analyze the subleading shape functions 
introduced at 0{as) by two insertions oi C^^^ . 

With all the definitions in place, we are now ready to evaluate the current 
correlator Tj^ including terms of up to second order in y/X, working at lowest order 
in as at the hard and hard-collinear scales. The leading term is readily found to 
be 

First-order corrections in vanish by rotational invariance in the transverse 
plane (provided we choose the coordinate system such that = and q±^ = 0), 
i.e., T^j^ — 0. At second order in the expansion the correlator receives several 
contributions, which can be represented symbolically as 

jt(2)j(o)^ jtWjW, Jt(0)j(2), (3.15) 




Figure 3.1: Representative examples of time-ordered products contributing to the 
current correlator Tij in SCET. Double lines show heavy-quark fields, dashed ones 
hard-collinear fields, and wavy lines denote the external currents. Lagrangian inser- 
tions and higher-order effective current operators are exemplified by soft gluons. 



and 



(3.16) 

Jt(0) J(0) J ^4^^(l) J ^4^£(1)^ (3.17) 

where at tree level only the type-A current operators appear. Examples of these 
time-ordered products are depicted in Figure 3.1. At 0{as), one must include a 
tree-level contribution of the form 



Jt{0) J{0) 



I dhcf^ J d'wLf^, (3.18) 



shown by the last diagram in the figure. Beyond tree level, one would also have to 
include the contributions from the type-S current operators. 

Because the external 5-mcson states in the definition of the hadronic tensor 
contain only soft constituents, and because the hard-collinear fields have been 
decoupled from the soft fields in the leading-order SCET Lagrangian (3.6), it is 
possible to contract all hard-collinear fields in the time-ordered products. At tree 
level, all we need is the hard-collinear quark propagator (3.7). Derivatives acting 
on hard-collinear fields give powers of p in momentum space, whereas components 
of x'^ appearing in the multipole-expanded expressions for the effective Lagrangians 
and currents can be turned into derivatives d/dp acting on the momentum-space 
amphtudes. Each insertion of a SCET Lagrangian correction introduces an inte- 
gral over soft fields located along the n light-cone. Consequently, the time-ordered 
products in (3.15) lead to expressions involving bi-local operators as in (3.14), while 
those in (3.16), (3.17), and (3.18) also lead to tri- and quadri-local operators. At 



first sight, this would seem to require the introduction of comphcated sublead- 
ing shape functions depending on up to three momentum variables uj,i, defined in 
terms of the Fourier transforms of the matrix elements of the non-local operators. 
However, the non-localities can be reduced using partial-fraction identities for the 
resulting hard-coUinear quark propagators. At tree level, it suffices to define shape 
functions of a single variable uj. 

To see how this works, consider the effect of an insertion of the Lagrangian 
Cf in (3.10). Since all hard-collinear fields must be contracted, we may consider 
without loss of generality the expression 

L{yi,y2) = I J d'zA^iy, - z)^ z^^n'' {S^gG,^S)^_ A^{z - y^) . (3.19) 

The fact that the soft fields live at position z^ implies that the two hard-collinear 
quark propagators carry the same momentum components n-p and p±. In analogy 
with the definition of the hard-collinear calligraphic gluon field, we now introduce 
the soft field [46] 



•^S//('^) 



u 

S\x) {iD^S){x) ^- j dtnT {S^gG^^S) {x + tn) , (3.20) 



which allows us to write {S'^ gG ^1,3) z_ as a derivative, —n-dzAsfx{z-.). Integrat- 
ing by parts in (3.19), we find 

[ ifi 
L{yi,y2) = i d^zA^{yi- z)-z'^As^,{z_)n-dzA^{z-y2) 

+ ijd^z[n- d, A^iyi - z)] | z'^Asf.{z.) A^{z - 1/2). (3.21) 

The hard-collinear propagator is a Green's function obeying the differential equa- 
tion 

{ji-d+^^A,{x)^p^'\x). (3.22) 

This allows us to write 

L{yu 2/2) = i y2±- As{y2-) - yi± ■ As{yi-) A^{yi - 7/2) 

-tjdh A^(yi -z)^z^- A,{z_) ^ A^{z - 2/2) 



d^z 
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z±-As{z-)A^{z-y2). (3.23) 



The terms involving transverse derivatives vanish at tree level, because they provide 
powers of transverse momenta of external lines, which are zero (recall that v± = 
and g± = 0). This leaves the terms shown in the first line, in which the z integral 



has been eliminated, and in which the product of two propagators has been reduced 
to a single propagator. In momentum space, these manipulations correspond to 
the partial-fraction identity 
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valid for two momenta that differ only in their n ■ p components. 

In the sum of all terms many cancellations and simplifications take place, and 
we find the rather simple result 
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In deriving these expressions, we have made use of the identity iD^ = S{id^+A'^)S^ 
with A'^ as defined in (3.20). Note that on the space of forward matrix elements the 
result for T-^ is hermitean, because we can integrate by parts and use translational 
invariance. 

In the last step, we can simplify the operator O2 by noting that the HQET 
equation of motion, iv ■ Dh — 0, along with — h, implies 



sHph = SUp^h + (^ - 



tn 



d{S^h) 



(3.27) 



It follows that all gauge-covariant derivatives of the heavy-quark fields are perpen- 
dicular derivatives. This fact restricts the number of subleading shape functions. 



3.3 Definition of subleading shape functions 

The Dirac structure of the operators 0„ can be simplified noting that the heavy- 
quark fields h of HQET are two-component spinor fields, so that between h . . . h 



the Dirac basis collapses to a set of four basis matrices (1, cr). In four-component 

notation, these are the upper-left 2x2 blocks of (1,775) (in the Dirac represen- 
tation). Instead of 775, we are free to take the matrices 7^75 and ^75. It follows 
that, between the = |(1 + ^) projectors supplied by the heavy-quark fields, any 
Dirac matrix F can be decomposed as 

r ^ 1 tr (r P,) 1 - ^ tr [r P„ a - 75] ^75 - ^ tr (r P, 7^^75) 7^^75 ■ (3.28) 
We denote by 

(S.../.>.<tM^-^Lffi^ (3.29) 

the forward P-meson matrix element of any HQET operator. Rotational invariance 
in the transverse plane implies that transverse indices can only be contracted using 
the symmetric and anti-symmetric tensors (we set 60123 = 1) 

= - ^^^^^4^ , 6f = 6-^"%.n,. (3.30) 
It follows that the only non- vanishing matrix elements are 

{i I d'zT{{hS)^{S^h)^ Lf{z)}) = ^ I due-i-^-s{u), 

J iii'b J 

{ihS)^Yi.7niSUD^^h)J = da;e-t---t(a;), 

n-x/2 

dt{{hS)^{SHD^^iD'iS)^^{S^h)J = ^- J du;e-'^-^-^u{u;), 
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(3.31) 
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If radiative corrections at the hard scale are included, it would be more appropriate 
to split up s{u!) — Skin ('^') + '-'mag Smagi^^), whcrc Cmag is the Wilsou coefRcieut of the 
chromo-magnetic operator in the subleading HQET Lagrangian. This ensures that 
the shape functions remain independent of the heavy-quark mass. The definitions 



of the functions t, u, v are chosen such that 
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ticu), 
u{u;) , 

v{uj) , 
(3.32) 



where [x,y] = S{x)S'^{y) is a product of two infinite-length soft Wilson lines, 
which on the light cone (i.e., for x,y \\ n) collapses to a straight Wilson line of 
finite length connecting x and y. 

Wc also need a variation of the first matrix element in (3.32), in which the 
derivative is located at position 0. Using hermitean conjugation, translational 
invariance, and the reality of tiuj), which follows from parity and time-reversal 
invariance of the strong interactions, we find that 



{{hiIj)A_) (0) i/i [0, x_] /i(x_)) = (/i(0) i/i [0, x_] {ipi_h) (x_)) , 



(3.33) 



implying that all terms containing a single insertion of can be related to the 
function t{uj). Prom this relation, it follows that 

j due-i-^-Hiu) = {{hS)^^[4,^{x^)-4,^{0)\ {S^h)J 
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dt {h{0) '- [0, tn\ ^^X9G^,.{tn) [tn, xJ\ h{x-)) , 



(3.34) 



which defines the function t{uj) in terms of a matrix element of the field-strength 
tensor. 

It is now straightforward to express the forward matrix element of the current 
correlator in terms of shape functions. The resulting traces of Dirac matrices 
can be simplified using identities for the e^^ tensor derived in [49]. Taking the 
imaginary part, we obtain for the hadronic tensor 
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where now p — rribV — q, and 
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(3.36) 



It follows from (3.35) that the subleading shape functions s(u;) and ^(a;) always 
come with the same trace as the leading shape function S{uj). However, u{ijj) is 
divided by the kinematic variable n-p, and so it does not enter in a universal (i.e., 

process- independent) combination with S{uj). 

Wc can now specialize our result to the case of semileptonic decay, for which 
Ti = 7^(1 — 75) and Tj = 7^(1 — 75). This yields 



= J dw d{n ■ p + cu) !^ {n^'v" 



+ n'^v'' - g 



le 



pua(3 



n. 



X 



1 + s(u;) + + + u{uj)-v{uj) 

nib J TTT'b ^ • P 



n • p 



UJ S{uj) t{ijj) t{ijj) ^ 



+ 



nib 



nib 



n ■ p 



(3.37) 



Similarly, for the contribution of the dipole operator Q^^ to B ^ Xsj decay, the 
Dirac structures are (up to prefactors) Fj = |[7^, — 75) and Tj = ^[fi, 7^](1 + 
75) , where the indices /i, u are contracted with the transverse polarization vector 
of the photon. In this case we obtain 
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(3.38) 



We stress, however, that while at leading power in Aqcd/^6 the dipole operator 
gives the only tree-level contribution to the B — > Xg^ decay rate, this is no longer 
the case when power corrections are included. For instance, interference terms 
of the dipole operator with current-current operators can lead to new subleading 
shape functions even at lowest order in perturbation theory. To derive these struc- 
tures, it would be necessary to match the entire effective weak Hamiltonian for 
B — > Xg^ decay onto SCET operators [50]. This task still has to be completed be- 
yond the leading order in A. Contrary to claims in [51, 52], a complete description 
of tree-level subleading shape-function effects in ^ ^ Xg^ decay is therefore still 
lacking. 



3.4 Moment relations and comparison with the literature 



Moments of the shape functions can be related to forward 5-meson matrix elements 
of local HQET operators [15] . In particular, setting x — in the defining relations 
(3.31) yields expressions for the normalization integrals of the shape functions. 
They are 



duj S{u!) 



dw {s{u!) , t{u!) , U{U!) , V (uj)} = . 



(3.39) 



The vanishing of the norm of all subleading shape functions is a consequence of 
Luke's theorem [53], and it ensures that there are no first-order AqcD/mi, correc- 
tions to total inclusive decay rates. For the functions t, u, v this is an obvious 
consequence of the fact that the integration domain in (3.34) and (3.32) shrinks 
to zero in the limit x — > 0. The interpretation of (3.39) is that subleading shape 
functions lead to local distortions of inclusive spectra, which cancel out when the 
spectra are integrated over a sufficiently large region in phase space. The first mo- 
ments characterize the strength of the distortions, while higher moments determine 
their shape. 

Taking a derivative in ■ dx in the definitions (3.31) brings down a factor of 
u; under the integrals on the right-hand side. Setting then x — > yields a set 
of relations for the first moments of the shape functions. The resulting matrix 
elements can be evaluated by means of the relations [54] 
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{^ J d^zT{CMD'^hmC^^\z)}) = -V^ {C'^\o)) = ^^^^^ ^ (3-41) 

where Ai and A2 are the familiar HQET parameters arising in the parameterization 
of second-order power corrections to inclusive decay spectra, and Cmag = 1 at tree 
level. We obtain 
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dujcu s{u!) 



Ai + 3A2 



J diUUu{uj) = : 



J dujujt{uj) — — A2 , J duj uj v{uj) — X2 ■ (3.42) 

We also recall that the first two moments of the leading shape function are 
/ du;u;S{Lu) = and / du;u;'^S{u;) = -Ai/3 [15]. 

3.5 Contributions from four-quark operators 



The last diagram in Figure 3.1 shows a tree- level contribution to the hadronic 
tensor involving two insertions of the subleading SCET Lagrangian £^^^ . The 



exchange of a hard-coUinear gluon implies that this graph is of order gi^, and so 
it vanishes in the hmit where ag is set to zero. The fact that the suppression 
factor is vra^ instead of ag/n reflects the phase-space enhancement of four-quark 
tree-level graphs compared with loop diagrams [55]. One might therefore expect 
that the four-quark contribution is numerically as important as the other tree-level 
subleading shape-function contributions. 

Applying the partial- fraction identity (3.24) twice, we find that the resulting 
contribution to the correlator T^"^^ can be written as in (3.25), adding a fifth oper- 
ator to the sum. It reads 

n-x/2 n-x/2 

I'' ' P J J 

tl 

(3.43) 

where ta are the generators of color SU(A^c)- Note that the field insertions are 
ordered according to "light-cone time" x_, just as they appear in the Feynman di- 
agram in Figure 3.1. This is a general result. Because the minus components n-phc 
of hard-collinear momenta are large, of order nih, hard-collinear fields always prop- 
agate forward in light-cone time. Turning, for instance, a forward-moving hard- 
collinear quark into a backward-moving hard-collinear anti-quark would require a 
hard quantum fiuctuation, which is already integrated out in SCET. As a result, 
Feynman amplitudes in SCET are ordered with respect to light-cone time, and 
that ordering is preserved in the matching onto HQET. This discussion explains 
why all our operators have the property that the coordinates z_ of soft fields range 
from to x_ in an ordered fashion. The results can therefore always be expressed 
in terms of bi- local operators depending only on and x_. However, at present 
we cannot exclude the possibility of non-trivial weight functions under these inte- 
grals, which could arise at higher orders in perturbation theory. If present, they 
may require a generalization of our definitions of subleading shape functions. 

Returning to the case of the four-quark operator in (3.43), we note that its 
contribution vanishes in the vacuum-insertion approximation due to the color- 
octet structure of the heavy-light quark bilinears. While this approximation is 
admittedly naive, phenomenological evidence based on studies of 5-meson lifetimes 
[55, 56, 57] and theoretical arguments based on lattice calculations [58, 59] and 
QCD sum rules [60] support the notion that matrix elements which vanish in 
the vacuum-insertion approximation are numerically suppressed, typically by an 
order of magnitude. It is thus very unlikely that the four-quark operator in (3.43) 
could give a larger contribution than loop-suppressed 0(0;^) corrections, which 
we have neglected. Note also that the contribution of the corresponding local 
operator to total inclusive rates, which is what remains when the corresponding 
shape functions are integrated over a sufficiently large domain, is bound to be 
tiny. The effect is of order (AgcD/^^b)^, and it is ctg-suppressed with respect to a 
four-quark contribution discussed by Voloshin [61], whose effect on the total decay 
rate is believed to be at most 3%. 



Nevertheless, it is interesting to study the structure of the operator O^lx^) in 
more detail and define corresponding subleading shape fTinctions. The decomposi- 
tion (3.28) implies that the Dirac structure can be rearranged in the form (omitting 
factors of S, 5^ and color indices for simplicity) 

U ^ \ « ,^ 1 + 




^iiip QQili^j = -2 [ i^j> I (iiii^j^^r- 1 (1^-1^)75 | ^iiip q, 

(3.44) 

where the notation implies a sum over the three rows in the equation. In the next 
step, we use that between ji . . .iji any Dirac matrix can be decomposed as 

r - I tr(r tr(r ^75) 1 75 - ^ tr(r i^i) | . (3.45) 

Only the first two terms survive after contraction of the index p in (3.44). At this 
stage we are left with four-quark operators with Dirac structures 

hhq'^{-f5)q, h-^-fshq-^iys) q , ^17^75 /i ^yi (75) g , (3.46) 

where (75) means 1 or 75. Lorentz and parity invariance imply that only the two 
Lorentz-scalar structures hhq-^q and h yi-75 h q ^75 q can have non-zero forward 
matrix elements between S-meson states. Putting back color factors and soft 
Wilson lines, we define the corresponding subleading shape functions 

n-x/2 n-x/2 
ti 

n-x/2 n-x/2 

2(-0^ / dh I ^t2([(/^^)ota],^75[ta(^^/^),JJ(g^),,J,#75[(^M*,n].) 
ti 

= y"du;e-t---/,(u;), (3.47) 

where k, I are color indices. Using the same notation as in (3.35), the corresponding 
contributions to the hadronic tensor are given by 
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(3.48) 



with the same traces as defined in (3.36). As far as the spin structure is con- 
cerned, the four-quark contributions can thus be absorbed into a redefinition of 
the subleading shape functions u and v, namely 

u{u;) = u{uj) - Tras fu{i^) , v{uj) = v{uj) + nag fv{i^) ■ (3.49) 



Note that the definitions (3.47) imply that the normahzation integrals as well as 
the first moments of the functions fu{uj) and /,t,,(cj) vanish (because the integration 
domain is of second order in n ■ x), and therefore the new functions u and v have 
the same first moments as the original ones, see (3.42). 

While the shape functions S, s, t, u, v are expected to be identical for charged 
and neutral B mesons up to tiny isospin-breaking corrections, this is no longer the 
case for the four-quark shape functions /„ and f^. The values of these functions will 
depend crucially on whether the light-quark fiavor q in the four-quark operators 
matches that of the i?-meson spectator quark. In the semileptonic decay B — > 
Xu the difference between the subleading shape functions /„ and /„ for B~ 
and B^ mesons is hkely to be one of the dominant sources of isospin-breaking 
effects on the decay distributions. 



3.6 Applications 

We are now ready to summarize our results. We absorb the contributions from 
four-quark operators into the functions u{u!) and v{u!) defined in (3.49). The 
moment constraints on the shape functions can be summarized by their expansions 
in distributions [15], which read 

S{uj) = S{uj)-^S"{u;) + ... , s{u;) = ^'^^' d'iu;) + ..., (3.50) 

2A 

t{u) = \25'{u) + ..., uiu;)^-^5'iu) + ... , viu) ^ -MS'iu) + . . . . 

These expressions allow us to test our results against existing predictions for in- 
clusive spectra obtained using a conventional heavy-quark expansion. 

The analytic properties of the shape functions are such that they have support 
for — oo < cu < Aqo, where is the asymptotic value of the mass difference 
{triB — "^fe)mi,-^oo in the heavy-quark limit. This parameter differs from the physical 
value of A by power-suppressed terms, 

A = ms - m(, = Aoo - + ... . (3.51) 

We would like the support in uj to extend over the physical interval — oo < a; < A, 
since this will ensure that the kinematic boundaries for decay distributions take 
their physical values set by the true S-meson mass. This can be achieved by 
shifting the arguments of all shape functions by a small amount Ac<j = |(Ai -|- 
3X2) /mb. For the subleading shape functions this changes nothing to the order we 
are working, since t{u! + Aa;) = t(a;) -|- . . . etc., where the dots represent terms 
of higher order in l/m?,. For the leading-order shape function, however, this shift 
produces a new 1 / nih correction: S{u) = S{u + Acu) — Acj S'{uj + Au) + . . . , where 
the prime denotes a derivative with respect to the argument. Using the fact that 



S{uj) and s{uj)/mh always appear together, we can absorb the extra term into a 
redefinition of the subleading-shape function s, defining a new function 



SoM = s{uj) - -(Ai + 3A2) S'{uj) . 



(3.52) 



Prom (3.42), it follows that the first moment of the function sq vanishes. In terms 
of these definitions. 



S{uj + Acu) = S{uj + Acu) + 



soiu + Ac) ^ ^^^^ ^ ^ ^ 



nib 



nib 



(3.53) 



In other words, our expressions for the hadronic tensors in (3.37) and (3.38) remain 
valid when in all shape functions the argument is shifted from a; to a; + Aa;, except 
that the subleading shape function s must be replaced with the redefined function 
So, whose norm and first moment vanish. Once this is done, the integrals over u 
extend from —00 up to the physical value of A. 

It is now convenient to express the shape functions as functions of a) = A — a; 
[63]. They have support in this variable for < a) < 00. The 5-functions in 
the tree-level expressions (3.37) and (3.38] for the hadronic tensors set a) = P+. 
We denote functions of a) by a hat, e.g. S{u;) = S{A — u; + Alo) — S{Aoo — '^), 
t{uj) = t(Aoo — (jj), and similarly for the other functions. In the equations below, 
A always refers to the physical parameter defined with the true 5-meson mass. 

Based on equation (3.37), we can write the structure functions Wi for B —>■ 
Xu iP decays as: 
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where y ^ n • p/nib- 

For B — > Xg^ we have from equation (3.38): 



wj: = -2 



nib 



-t(P+) + uiP+) - v{P+) 
nib 



(3.55) 



where for radiative decays n • p — nib- 



3.7 Conclusions 

In this chapter, we have used the formahsm of SCET to perform a systematic study 
of power-suppressed effects to charmless inclusive B decays. At tree level, the 
results can be expressed in terms of a set of subleading shape functions defined via 
the Fourier transforms of forward matrix elements of bi-local light-cone operators 



in heavy-quark effective tfieory. We fiave identified a new contribution arising 
from four-quark operators, wliicli was not considered previously. We liave sliown 
tliat, wlien sliape functions appearing in process-independent combinations are 
combined into single functions, then a total of three subleading shape functions 
are required to describe arbitrary current-induced decay distributions of B mesons 
into hght final-state particles. 

In the last part of this chapter, we have presented analytical expressions for 
the hadronic tensor of semileptonic and radiative charmless B decays. While this 
concludes the problem of tree-level power corrections in semileptonic decay, we have 
stressed that no complete (tree-level) analysis of power-suppressed corrections to 
the B — > Xg^ decay exists to date. The formalism developed in this chapter can, 
however, readily be extended to this case. 



Chapter 4 

Applications: Event Generator 



4.1 Introduction 



The theoretical tools for the calculation of inclusive B decays are QCD factorization 
on the one hand [15, 16, 17, 18, 25, 50, 62, 63, 64], and local operator product 
expansions (OPE) on the other [1, 13]. Both approaches perform a systematic 
separation of long-distance hadronic quantities from short-distance perturbative 
ones, while organizing the calculation in inverse powers of the heavy 6-quark mass 
mfj. The OPE is an appropriate tool for the calculation of total inclusive rates (for 
example in S — > Xcl~i^ decay) or for partial rates integrated over sufficiently large 
regions in phase space, where all components of the final-state hadronic momentum 
Px are large compared to Aqcd- QCD factorization, on the other hand, is better 
suited for the calculation of partial rates and spectra near kinematical boundaries, 
where typically some components of are large, while the invariant hadronic 
mass Mx = -J P'x small. 




It is important to note that the heavy-quark expansions valid in these two 
kinematical regions are not identical, because the power counting rules differ in 
the two regimes. Also the nature of the non-perturbative inputs is different. In the 
OPE region, non-perturbative physics is encoded in a few hadronic parameters, and 
the heavy-quark expansion is the usual Wilsonian expansion in local operators. In 
the endpoint (or shape- function) region, the presence of multiple scales complicates 
the power counting, and the interplay between soft and collinear modes gives rise to 
large non-localities. As a result, non-perturbative physics is described by hadronic 
structure functions called "shape functions" , and the heavy-quark expansion is an 
expansion in non-local string operators defined on the light-cone. The connections 
between the two regimes is that moments of the shape functions can be expressed 
in terms of local operators. 

In this chapter we develop a formalism that smoothly interpolates between 
the two kinematical regimes (see [65] for a related discussion, which is however 
restricted to the tree approximation). This is essential for building an event gen- 
erator for inclusive B X„ l^v and B Xg'y decays, which can be used to study 
partial and differential decay rates in different kinematical domains. In the shape- 
function region, our approach relies on exact QCD factorization theorems, which 
exist in every order of power counting. They allow us to systematically disentangle 
short- and long-distance physics and, in the process, resum parametrically large 
logarithms order by order in perturbation theory. This factorization can be done 
with high accuracy for the terms of leading power in l/m;,, and with somewhat 
less sophistication for the first-order power corrections. For the second-order power 
corrections, we only include contributions that do not vanish when integrated over 
all phase space. This is a safe approximation; the effects of the remaining 
terms can to a large extent be absorbed by a redefinition of the subleading shape 
functions arising at order l/m^. 




Our formalism is "optimized" for the shape-function region in the sense that 
sophisticated theoretical technology is applied in this regime. However, when our 
expressions for the differential decay rates are integrated over sufficiently wide 
domains, they automatically reduce to the simpler results that can be derived using 
the OPE approach, up to yet unknown terms of 0(0^). The moment relations for 
the shape functions are crucial in this context. Note that local I /ml corrections 
in the OPE receive contributions from terms of leading power {1/m^), subleading 
power {l/rrih), and sub-subleading power (1/m^) in the shape- function region, so 
the transition is highly non-trivial. In implementing the program outlined here, we 
include all presently known information on the triple differential B — > Xu I'P decay 
rate and on the differential B — > decay rate in a single, unified framework. 
We neglect, for simplicity, hadronic power corrections of order l/mf and higher, 
which are known to have a negligible effect on the observables considered here. 
The only possible exception is contributions from "weak annihilation", which are 
estimated as part of our error analysis. We also ignore the existing results on 
0{Poal) radiative corrections for some single-differential distributions, because the 
corresponding corrections are not known for the double or triple differential B 
Xu I'v decay spectra. While these 0(/5octs) terms are sometimes found to be large 
when naive perturbation theory in asirrib) is used, their effects are expected to be 
small in our scheme, which is based on a complete scale separation using QCD 
factorization. We see no reason why the PqoI terms should be enhanced compared 
to other, unknown corrections of 0{af). 

A technical complication in realizing the approach described here has to do 
with the treatment of phase-space factors. The heavy-quark expansion of the 
hadronic tensor for B ^ X^ I'v decay gives rise to expressions that are singular 
at certain points in phase space. One way to avoid these singularities is to also 
expand phase-space factors order by order in l/mt (see, e.g., the treatment in 
[66]). However, since this expansion depends on the kinematical cuts of any given 
analysis, it cannot be implemented in a straightforward way in an event genera- 
tor. An alternative is to reorganize the heavy-quark expansion in such a way that 
the expansion parameter is related to hadronic (as opposed to partonic) kinemat- 
ical variables, in which case kinematical singularities are always canceled by exact 
phase-space factors. Following this strategy, we obtain expressions for decay dis- 
tributions and partial decay rates which are free of explicit reference to partonic 
quantities such as the 6-quark mass. A dependence on mi, enters only implicitly 
via the first moment of the leading-order shape function 5'(tD). The philosophy of 
our approach is that this function^ is extracted experimentally from a fit to the 
B — > photon spectrum, which has been measured with good precision in the 
region where = Mb — 2E^ ~ Aqcd- This is analogous to the extraction of 
parton distribution functions from deep inelastic scattering. The photon spectrum 

^More precisely, we define a new shape function 5 (a)) by the combination of 
leading and subleading shape functions contributing to S — > Xgj decay, and we 
will use the same function to make predictions for B — > X^ I'V decay distributions. 



is experimentally accessible to energies as low as 1.8 GeV, which corresponds to a 
sampling of the shape function for values of Co up to about 1.7 GcV. Once the shape 
function has been extracted over this range, we can use it to obtain predictions 
for arbitrary partial B — > Xul^v decay rates with cuts. In doing so, the residual 
hadronic uncertainties in the extraction of \Vub\ only enter at the level of power 
corrections. 

We emphasize that the program outlined above is equivalent to an approach 
put forward in [16] and later refined in [67, 68, 69], in which \Vuh\ is extracted with 
the help of shape-function independent relations between weighted integrals over 
differential decay distributions in S — > Xg^ and B — > I'V. The experimental 
error in the results for these weighted integrals corresponds, in our approach, to the 
error in the prediction of S — > X^ partial rates resulting from the experimental 
uncertainty in the extraction of the shape function from the B photon 
spectrum. While the shape-function independent relations are very elegant, it is 
more convenient for the construction of a generator to have a formulation where the 
shape function is used as an input. In this way, it is possible to impose arbitrary 
cuts on kinematical variables without having to recompute the weight functions in 
each case. 

The chapter is structured as follows: In Section 4.2 we collect the relevant 
formulae for the calculation of the S ^ X57 photon spectrum. These expres- 
sions can be used to extract the leading non-perturbative structure function from 
experiment. An analogous presentation for the triple differential decay rate in 
B Xul^v decays is presented in Section 4.3. In order to perform a numerical 
analysis one needs to rely on parameterizations of the shape functions. A collection 
of several useful functional forms is given in Section 4.4. In Section 4.5 we present 
a full error analysis of partial B — > X^ decay rates for a variety of experimen- 
tal cuts. We also explore the sensitivity of the results to the &-quark mass and 
to the functional forms adopted for the shape functions. Section 4.6 contains our 
conclusions. 

4.2 Inclusive radiative decays 

The decay process B — > Xs7, while more complex in its short-distance physics, is 
considerably simpler in its kinematics than the semileptonic process B — > X^ I v . 
Since the radiated photon is on-shell, the hadronic variables P± that describe the 
momentum of the Xg system are trivially related to the photon energy E'^ by 
Pj^ — Mb — 2Ej and P_ = Mb- In the crudest approximation, namely at tree 
level and leading power, the photon-energy spectrum is directly proportional to the 
leading shape function, dTg/dE^ oc S{P+). In this section we collect all relevant 
formulae needed to compute the B Xg^ photon spectrum or, equivalently, the 
invariant hadronic mass distribution. It is implicitly assumed that these spectra 
are sufficiently "smeared" (e.g., by experimental resolution) to wash out any sharp 
hadronic structures. In cases where the resolution is such that the K* resonance 



fit; 



peak is observed, it can be accounted for by combining the formulae in this section 
with the prescription for subtracting the K* peak proposed in [32]. 
The differential B —>■ Xgj decay rate can be written as 

^ = ^ iVM'm'M [Cf^innW U{n,, 3^,(P+) , (4.1) 

where the structure function depends on the photon energy via P+ = Mb — 2E^. 
The prefactor contains the electromagnetic fine-structure constant a normalized 
at — 0, two powers of the running 6-quark mass (defined in the MS scheme) 
originating from the electromagnetic dipole operator Qj^ in the effective weak 
Hamiltonian, and the square of the corresponding Wilson coefficient Cj^, which is 
needed at next-to-leading order in renormalization-group improved perturbation 
theory [10]. Renormalization-group running from the hard scale /ih ~ iTib to the 
intermediate scale yUj ~ a/^6Aqcd gives rise to the evolution factor U{iik,l^i), 
whose explicit form is discussed in Appendix B. We keep U and (Cy^)^ outside 
of the structure function 3^^; it is understood that when combining the various 
terms in (4.1) all perturbative quantities should be expanded for consistency to 
the required order in a^. 

4.2.1 Leading-power factorization formula 

At leading order in l/nib the structure function factorizes as [50] 

(P+) = \Hs{i^h)\'' / dwmbJ{mb{P+-u;),iii)S{u;,iii). (4.2) 
Jo 

At this order a single non-perturbative parton distribution function arises, called 
the leading shape function [16] and denoted by S{uj,fii). Our notation is adopted 
from [63, 66]: hatted shape functions have support for > 0. The function S 
is defined in terms of a non-local matrix element in heavy-quark effective theory 
(HQET). Renormalization-group running between the intermediate scale and a 
low hadronic scale is avoided when using the shape functions renormalized at the 
intermediate scale /ij. Evolution effects below this scale are universal (i.e., process 
independent) and so can be absorbed into the renormalized shape function. Short- 
distance contributions from scales above /x/j ~ nib are included in the hard function 
Hs, which in practice is obtained by matching the effective weak Hamiltonian onto 
a current operator in soft-coUinear effective theory (SCET). At next-to-leading 



order in perturbation theory, the result reads 
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where the variable z = {rric/mbY denotes the ratio of quark masses relevant to 
charm-loop penguin diagrams, and the "penguin function" g{z) can be approxi- 
mated by the first few terms of its Taylor expansion. 
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The Wilson coefficients Ci and CgJ in (4.3) multiply the current-current oper- 
ators Q\'^ and the chromo-magnetic dipole operator Qsg in the effective weak 
Hamiltonian. The quantities e^w ~ —1.5% and Spcng ~ —0.6% account for small 
electroweak corrections and the effects of penguin contractions of operators other 
than Qi'*^, respectively. The differential decay rate (4.1) is formally independent 
of the matching scales and /ij. The /i^ dependence of the evolution factor 
U{nh,Hi) cancels the scale dependence of the product YnKnh) [C^^{nhW l^sil^h)]'^, 
while its dependence compensates the scale dependence of the convolution in- 
tegral J(/ii) S{fii). 

Finally let us discuss the jet function J, which appears as the hard-scattering 
kernel in the convolution integral in (4.2). It can be written in terms of distributions 
that act on the shape function S. At one-loop order, the jet function is given by 
[62, 63] 
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where the star distributions have the following effect on a function / when inte- 



grated over a domain [2] : 
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(4.6) 



4.2.2 Kinematical power corrections 

There exists a class of power corrections to (4.2) that do not involve new hadronic 
quantities. Instead, the power suppression results from the restriction of certain 
variables (P+ in the present case) to a region where they are kinematically sup- 
pressed (here P+ <S Mb). The corresponding terms are known in fixed-order 
perturbation theory, without scale separation and renormalizat ion-group resum- 
mation [70, 71] (see also [32]). To perform a complete RG analysis of even the 
first-order terms in l/rrib is beyond the scope of this work. Since, as we will see 
later, power corrections only account for small corrections to the decay rates, an 
approximate treatment will suffice. To motivate it, we note the following two facts 
[50]: First, while the anomalous dimensions of the relevant subleading SCET and 
HQET operators are only known for a few cases [47], the leading Sudakov double 
logarithms are the same as for the terms of leading power, because they have a 
geometric origin in terms of Wilson lines [72] . The leading Sudakov double log- 
arithms are therefore the same as those resummed into the function U in (4.1). 
Secondly, the kinematical power corrections m B ^ Xg'-y decay are associated with 
gluon emission into the hadronic final state Xg. Because of the kinematical restric- 
tion to low-mass final states, i.e. Mx ~ MbAqcd, wc associate a coupling as^Jl) 
with these terms, where typically ft fj^i. Strictly speaking, however, the scale 
ambiguity associated with the choice of p, could only be resolved by computing the 
relevant anomalous dimensions. 

Within this approximation, the kinematical power corrections to the structure 
function 3^^ can be extracted from [32, 50]. We find it convenient to express the 
result in terms of the variable 



which in the shape- function region scales like Aqcd/?7t.&. We obtain 
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The coefficient functions hij{x) are 
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In the shape-function region the expressions for hij{x) could, if desired, be ex- 
panded in a power series in x = 0(AQCD/''^fe), and this would generate a series 
of power-suppressed terms 9^^™^"^(P+) with n > 1, where the superscript "n" in- 
dicates the order in the l/nib expansion. Note that this expansion would contain 
single logarithms Ina; ~ ln(AQCD/'"^6)- These are precisely the logarithms that 
would be resummed in a more proper treatment using effective field-theory meth- 
ods. 

Outside the shape-function region the variable x can take on arbitrarily large 
positive values, and ?'!^™(P+) is no longer power suppressed. Note that for P_|_ — > 
Mb (corresponding to x ^ cxd and — > 0) most functions hij{x) grow like x"^ or 
weaker, so that the spectrum tends to a constant. The only (well known) exception 
is hssix), which grows like x^, giving rise to a l/E^ soft-photon singularity [71]. 
The main effect of the kinematical power corrections (4.8) to the photon spectrum 
is to add a radiative tail extending into the region of small photon energies. These 
corrections therefore become the more significant the larger the integration domain 
over is. 



4.2.3 Subleading shape- function contributions 

At order l/m^, in power counting, different combinations of subleading shape func- 
tions enter the B ^ Xgj and B — > Xul~u decay distributions [51, 52, 73, 74]. 
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They provide the dominant hadronic power corrections, which must be combined 
with the kinematical power corrections discussed in the previous section. We in- 
clude their effects using the rcsuhs of recent calculations in [66, 75, 76]. Little is 
known about the subleading shape functions apart from expressions for their first 
few moments. In particular, the norms of these functions vanish at tree level, while 
their first moments are determined by the HQET parameters Ai and A2, which are 
defined via the forward S-meson matrix elements of the kinetic-energy and the 
chromo- magnetic operators, respectively [54]. 

For the case of i? — > Xg"") decay, subleading shape-function contributions are 
currently only known for the matrix elements of the dipole operator Qy-y, and the 
corresponding hard and jet functions have been computed at tree level. Adopting 
the notations of [66], the relevant subleading shape functions arc t(cj), m(cj), and 
v{Co). An additional function, called Sq, has been absorbed by a redefinition of 
the leading shape function, and it is included in our definition of S{(li). Roughly 
speaking, ■u(a)) is the "light-cone generahzation" of the local HQET kinetic-energy 
operator. The functions v{Lb) and i{Cj) are both generahzations of the local chromo- 
magnetic HQET operator, but t{Cj) contains also a light-cone chromo-electric op- 
erator, which has no equivalent in the local OPE expansion. (Such a contribution 
arises since there are two external 4-vectors in the SCET expansion, n and 
while there is only v in the HQET expansion.) The contribution of subleading 
shape functions to the B — > Xg^ photon spectrum is 

Compared to [66], we have replaced l/m^, with 1/{Mb — P+) in the prefactor, which 
is legitimate at this order. (The form of the shape functions restricts P+ to be of 
order Aqcd-) The appearance of the HQET parameter A = {Mb — rnb)mb^oo is 
peculiar to the subleading shape-function contributions. This quantity is defined 
via the first moment of the leading-order shape function [63]. 

The formula given above can be modified to suit the purpose of extracting the 
shape function from the photon spectrum better. To this end, we absorb a linear 
combination of the subleading shape functions into a redefinition of the leading 
shape function, in such a way that the moment relations for this function remain 
unchanged to the order we are working. This is accomplished by defining 

Siif, ^ S{^) + - ^<''> - + - . (4.12) 

rub 

When using S instead of 5* in the leading-power formula (4.2), the subleading 
shape-function contribution becomes 

^r^^)(P,) = -^^^5(P+). (4.13) 
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The hatted shape functions used in this chapter are related to the original 
definitions in [66] by 

oicu) — S{A — u!)-\ ^ -, 

nib 

i{u}) = t{A - u) , u{u}) = u{A - u) , v{u}) = v{A - u) , (4.14) 

where the unhatted functions have support on the interval between — oo and A. It 
is convenient to rewrite A — cu — cu + Acu, where 

Au = A-{MB-nib)^ 'y' +... (4.15) 

2777,6 

accounts for the mismatch between the HQET parameter A and the difference 
(Mb — nib) due to power-suppressed terms in the 1/nib expansion [11]. It follows 
that the variable uj = {Mb — nib) — ^ I'uns from — oo to [Mb — n^-b)- The moment 
relations for the leading and subleading shape functions derived in [16] and [51, 66] 
can be summarized as 

S (u) = Siu; + Au) H ^ ^ d(uj) - — 6 (u) -\ (u) + . . . , 

nib O 27776 

i{u;) = t{u + Au) = A2 5'{ijo) + . . . , 

2Ai 

u{Cj) = u{u + Au) = — — S'{u) + . . . , 

E v{uj + Au) ^ -X2S'{uj) + ... . (4.16) 
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The function S has the same moment expansion as S. The hadronic parameter A2 
determines the leading contribution to the hyperfine splitting between the masses 
of B and B* mesons through m^, — n^'s — 4A2 + 0(l/mft) [54], from which it follows 
that A2 ~ 0.12 GeV^. The value of the parameter Ai is more uncertain. In much the 
same way as the 6-quark pole mass, it is affected by infrared renormalon ambiguities 
[43, 44]. It is therefore better to eliminate Ai in favor of some observable, for which 
we will choose the width of the leading shape function. 



4.2.4 Residual hadronic power corrections 

At order l/7?7^ a new set of sub-subleading shape functions enter, which so far 
have not been classified completely in the literature. Since the functional form of 
even the subleading shape functions is rather uncertain, there is no need to worry 
too much about the precise form of sub-subleading shape functions. Most of their 
effects can be absorbed into the subleading functions. An exception, however, are 
terms that survive when the sub-subleading shape functions are integrated over 
a wide domain. Whereas the norms of all subleading (~ 1/nib) shape functions 
vanish, the norms of the sub-subleading shape functions (~ l/iTil) are in general 
non-zero and given in terms of the heavy-quark parameters Ai and A2. (At tree 
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level, the class of functions with non-zero norm has been studied in [65].) Our 
strategy in this chapter will be as follows: Wc start from the well-known expres- 
sions for the (tree-level) second-order power corrections to the B X^'-f photon 
spectrum [19] (and similarly for the triple-differential B decay distribu- 

tion [13, 1], see Section 4.3.4). They are of the form Xi/ml times one of the singular 
distributions 5(p^), 5'(p^), or S"{p'^), where — {mi,v — qY is the invariant partonic 
mass squared of the final-state jet. As mentioned earlier, the power counting in 
the shape-function region is different from the one used in OPE calculations, and 
indeed a good portion of the l/ml terms in the OPE is already accounted for by 
the contributions proportional to the leading and subleading shape functions in 
(4.2) and (4.11). We identify the corresponding terms using the moment relations 
for the shape functions in (4.16). In particular, this reproduces all terms at order 
l/mf in the OPE which contain derivatives of S{p^). We include the remaining 
terms of the form (Xi/ml) S{p^) by replacing 

S{p'^) = S{p+p_) ^ / duj6(p+ + u) 6{lj) 

P--P+ J 

- du5{P^-u)S{u)^^^^. (4.17) 

Here p± are the light-cone projections of the partonic momentum p^, which are 
related to the hadronic quantities P± by P± — p± -\- {Mb — rrib). Similarly, cu — 
{Mb - rub) - OJ. 

The result of these manipulations is 

Together with (4.2) and (4.11) this accounts for all known first- and second-order 
power corrections to the B Xg'-y photon spectrum, both in the shape-function 
region and in the OPE region. The redefinition (4.12) of the leading shape function 
from S io S leaves the form of the second-order power corrections unaffected. 

In Section 4.5 we study the numerical impact of second-order power correc- 
tions on various B — > X^l^v partial rates and find their effects to be tiny. It is 
therefore a safe approximation to neglect hadronic power corrections of order l/m\ 
or higher. The only possible exception to this conclusion relates to the so-called 
weak annihilation terms in 5 — > X^ decay, which will be included in our error 
analysis. 
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4.3 Inclusive semileptonic decays 



In terms of the Wi functions introduced in chapter 1, the triple differential decay 
rate for B X, l~i> reads 



Gp\Vub\ 



dP+ dP- dPi 



Uy{llh,l^i) {Mb-P+) 



167r3 

(P_ - Pi){Mb -P- + Pi- P+) ?i 



+ (Mb - P-)(P- - P+) 3^2 + {P- - Pi){Pi - P+) 3^3] , 



(4.19) 

where we have collected the relevant combinations of Wi into the three functions 

Wo 



Uy{l^h,l^i) 3^3 



4 y 



and defined a new kinematical variable 

P--P-^ 



y 



Mb-P^ 



(4.20) 



(4.21) 



which can take values < y < 1. The leading evolution factor Uy{iJ,h, fJ-i) has 
been factored out in (4.19) for convenience, as we have done earlier in (4.1). The 
function Uy{jih, l^i) differs from the corresponding function in P — > decay by 
a |/-dependent factor, 



(4.22) 



where the function or in the exponent is related to the cusp anomalous dimension 
and is given in Appendix B. 

Eq. (4.19) for the triple differential rate is exact. Note that there is no refer- 
ence to the 6-quark mass at this point. The only dependence on is through 
the structure functions 3^i{P+,y) (via hard matching corrections and via the mo- 
ment constraints on the shape function S), which are independent of the leptonic 
variable P;. The fact that the total decay rate r,u is proportional to is not in 
contradiction with (4.19). It is instructive to demonstrate how these five powers 
of nib are recovered in our approach. At tree level and leading power the functions 
3^2 and 3^3 vanish, while 3^i — S{P+). Integrating over the full range of Pi and P_ 
builds up five powers of {Mb — P+). For the purpose of illustration, let us rename 
the P+ variable to cD in the last integration, so that the total decay rate is given as 



G% I Vub I 
1927r3 

1927r3 



Mb 



du {Mb - (of S{Co) 



GWVub^^ 
1927r3 



Mb -mi, 



rrib 



{mb + {uj)f 



1 + 01 — 



duo {rrib + a;) S{lj) 



(4.23) 
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Figure 4.1: The hadronic phase space in P+ and P_. The hght gray region contains 
background from B ^ Xc decays, while the dark gray region is only populated 
by 5 ^ Xu events. The line separating the two regions is the contour where 
Mx = P+P- = M^. Each point represents a, B ^ Xul^u event in a Monte-Carlo 
simulation using the results of this chapter. While the shape-function region of 
large P_ and small P+ is highly populated, there is not a single event with P+ 
larger than 3 GeV out of the 1300 events generated. 



At tree level, the first moment of the shape function S{uj) vanishes. Beyond tree 
level this is no longer the case, and the average (uj) depends on the size of the 
integration domain. The above observation motivates the use of the shape-function 
scheme [63] , in which the 6-quark mass is defined as mf = mf^' + {u) + 0{l/mb). 
After this is done, (4.23) recovers the form of the conventional OPE result. 

Eq. (4.19) and the above argument tell us that the differential rate is a priori 
rather insensitive to the 6-quark mass in the endpoint region, where P+ (and 
therefore {lj)) is a small quantity. Only when the rates are integrated over a 
sufficiently wide domain, so that shape-function integrals can be approximated 
using a moment expansion, a dependence on enters indirectly via the first 
moment of the leading-order shape function. Likewise, a dependence on other 
HQET parameters such as Ai enters via the sensitivity to higher moments. 

In the remainder of this section we present the various contributions to the 
structure functions Hi, following the same line of presentation as we did in the case 
of S — > decay in Section 4.2. As before, while the resulting expressions are 
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"optimized" for the shape-function region, they can be used over the entire phase 
space and give the correct result for the total decay rate up to corrections of 0(«^). 
In the shape-function region, where P+ is a small quantity, one may organize each 
as a series in inverse powers of 1/{Mb — P+)- No assumption about the variable 
y is made, which is treated as an 0(1) quantity.^ A preview of the results of our 
calculation is depicted in Figure 4.1, which shows an illustration of our prediction 
for the distribution of events in the plane (P+, P_). 



4.3.1 Leading-power factorization formula 



The leading-power expressions for the hadronic structure functions Wi have been 
calculated in [63] at one-loop order in rcnormalization-group improved perturba- 
tion theory. At this level 3^2 docs not obtain a contribution, whereas and 3^^ do. 
Symbolically, they take the factorized form Hui J <^ S, consisting of hard functions 
Hui and the convolution of the jet function J with the leading shape function S. 
More precisely, 

^f\P+,y)^Hui{y,i^h) I dLbymbJ(ymb(P+-u),iJ,i)S(u,iii), (4.24) 
where the hard functions are given by 



Hui{y,Hh) 



14 



-41n^^ + 101n^ 



2\ny 



-4L2(1-|/)- — -12 

y 6 



47r 1 — y ' 



(4.25) 



and Hu2 

scales fih 



0. As before, the differential decay rate is independent of the matching 

' mf, and 



rrih Aqcy). The jet function J has already been given 
in (4.5). Note that the 6-quark mass appears only as the argument of logarithms, 
where it plays the role of setting the renormahzation scale. 



4.3.2 Kinematical power corrections 

As in the case of B — > Xsj decay, there is a class of power corrections to the 13 — > 
Xu I'P decay distributions which are small only because of the restriction to certain 
regions in phase space, but which are not associated with new hadronic parameters. 

In the present case, these terms can be extracted from the one-loop expressions 
derived in [2]. They are then convoluted with the leading shape function. As 
previously, the scale separation that can be achieved for these power-suppressed 

^In the shape- function region, where P+ ^ P_, we have y ~ p_/m{,, which is 
the variable used in the leading-power analysis in [63]. 
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terms is only approximate, and we thus assign a coupling as{p) with them, where 

the scale fi is expected to be of order /Xj ~ y'TrvVqcD- 

The resulting expressions for the structure functions can be written in a com- 
pact form in terms of the variables x and y defined in (4.7) and (4.21). We find 



1 



Mb 

X 



P+ Air 



duj S{uj, jii) 
'^9i{x,y) 



(1 + xyy{x + y) x{l + xyy'^{x + y) 
V X/ X \ X/ 



P+,y) 



Mb - P+ 47r 
f2{x,y) 



X 



2xg2{x,y) 



(1 + xyy'^{x + y) {l + xyy^{x + y) 

1 CVn„(/i) 



In 



Mb - P+ 47r 

Mx,y) 



X 



doj S{lij, Hi) 
'^93{x,y) 



{l + xyy^{x + y) {1 + x)'^y^{x + y) 



(4.26) 



where the functions /j, gi are given by 



fi{x,y) 
9i{x,y) 
f2{x,y) 
92{x,y) 
h{x,y) 

93{x,y) 



-9y + lOy^ + x{-16 + 12y + 6y^) + x'^{13y - 12) , 

-2y^ - 2xy'^{A + y) - x^y{l2 + Ay + y"^) - Ax^{y + 2) + 3x'^{y - 2) , 

1/2 + xy{^ + 4|/ + y2) + 3x^1/(10 + y) + ;x'^(12 + 19y) + lOx^ , 

2|/2 + 4:xy{l + 2y) + x2y(18 + Sy) + 6x-'^(l + 2y) + , 



2l/'(2y-ll) 



xy 



-94 + 29y + 2y'^) + 2x2|/(-72 + 18|/ + IS^/^) 



+ a;-'^(-72 - 42y + l^y"^ - 3y^) - 10x^(6 - 6y + y^) , 
4y* - 6x{y - 5)y^ - 4x'^y'^{-20 + 6y + y^) 
+x^y{90 - lOy - 28y^ + y^) + x^(36 + 36|/ - 507/^ + Ay^) 
+5x^(6 -6?/ + 1/2). 



(4.27) 



The above formulae are the exact 0(0;^) corrections to the leading-power expres- 
sion. This means that, when integrated over the entire phase space, they will give 
rise to the correct result for the total rate up to that order. In the shape-function 
region (where P+ <S P_) the integrands in (4.26) can be expanded in powers of 
l/mfe by counting y — 0(1) and x — 0{l/mb). Note that this organizes the l/nib 
expansion as an expansion in powers of the hadronic variable 1/{Mb — P+). The 



leading terms read 



1 CpoisiiA 

Mb - P+ 47r 

1 CpOisilj) 



p+ 



Mb - P+ 47r 

1 CpOLsiP) 



Mb - P+ 47r 



dCo S{lu, Hi 



dCj S{u;, Hi 



y \y J ^ 

i 
y 



\ 22 8, y' 
4 + - In - 

y y X. 



. (4.28) 



Further accuracy can be achieved by adding the next-order corrections, for which 
we obtain 



1 Cpasilj) 



{Mb - P+Y 47r 

X 

y y^ \i 

(Mb - P+f 47r 



duj (P+ — uj) S{uj, Hi) 



_ 16 3 /12 20 A, y 
-12 + — + — + — + 61n^ 

y r \y y J x 



dCj (P+ — Cj) S{u;, Hi) 



X 



^f''^'\P^,y) 



\ 2 7 4 v' 
l + - + ^- ^ln- 
y X 



1 Cpasin) 



{Mb - P+)2 47r ,,o 

X 



du (P+ — a)) 5'(ci), Hi) 



6 + 5?_!lW5?-?? 

y r Vr y J ^. 



(4.29) 



In the various phase-space regions of interest to the determination of iKibI, the 
above terms (4.28) and (4.29) approximate the full result (4.26) very well (see 
Section 4.5 below). 

Let us comment here on a technical point already mentioned in section 4.1. 
When combined with the phase-space factors in (4.19), the exact expressions for 
grkin (4.26) arc rcgTilar in the limit P_ — > P+, corresponding to ?/ ^ 0. However, 
this feature is not automatically ensured when the structure functions, but not the 
phase-space factors, are expanded about the heavy-quark limit. With our choice 
of the variables x and y, we encounter terms as singular as 1/y" at n-th order 
in the expansion, as is obvious from the explicit expressions above. Phase space 
scales like y^ in the limit y — (note that Pi — > P+ as P_ — > P+ because of 
(1.5)), so that the results (4.28) and (4.29) can be applied without encountering 
any kinematical singularities. In order to achieve this, it was crucial to define 
the variable y in the way we did in (4.21). We emphasize this point because 
straightforward application of the technology of SCET and HQET developed in 
[66, 75, 76] would give an expansion of the structure functions J'j in powers of l/p_, 
whereas phase space is proportional to 4p^ = (p_ — p+)^ oc y^ . In the kinematical 
region where p+ < 0, which is allowed due to off-shell effects in the B meson, this 
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leads to singularities as p_ — > 0. In order to avoid these singularities, we have 
reorganized the SCET expansion as an expansion in — P+) instead of 

where |p+| <C P- in the shape-function region. 



4.3.3 Subleading shape- function contributions 

The contributions from subleading shape functions to arbitrary B X„ l^u decay 
distributions have been derived (at tree level) in [66, 75, 76]. The results involve 
the same set of subleading shape functions as previously discussed in Section 4.2.3. 
Again, the structure function 3^2 does not obtain a contribution, while 



1 



Mn-P^ 



iiiPA - viP^ 



i-.hadr(l) 



1 



2 

'Ply 



-(A - P+) S{P^) - 2i{P^) + 



t(P+) + v{P^ 



(4.30) 



At this point we recall the discussion of Section 4.2.3, where we have argued that 
the B — > Xsj photon spectrum should be used to fit the function S of (4.12), 
which is defined to be a linear combination of the leading shape function S and 
the subleading shape functions t, u, v. When the above results are rewritten in 
terms of the new function S nothing changes in the expressions for !jf''' except for 
the simple replacement S ^ S, which we from now on assume. At the level of 
subleading shape functions j^^^^^^^ — q and j^^^'^^^^^ remain unchanged, while 



1 



Mn-P^ 



-{K-P+)S{P+) + 2i{P^ 



1 



+ [u{P^) - v{P+)] ( - - 1 



(4.31) 



It follows that there reside some linear combinations of subleading shape functions 
in the triple differential decay rate that cannot be extracted from information on 
the photon spectrum va. B ^ decays. In the end, this dependence gives rise 
to a theoretical uncertainty. 



4.3.4 Residual hadronic power corrections 

In analogy with our treatment for the case of i? — > X57 decay, we start from the 
expressions for the 1/m^ corrections to the triple differential B — > X^l'V decay 
rate obtained by applying the OPE to the hadronic tensor [1, 13]. We saw in 
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chapter 1 that: 



+^f!^ + 4AA. (4.32) 



P- \ 3 

The desired 1/{Mb — P+Y corrections to the structure functions can then be 
extracted by expanding the leading and subleading contributions 3^f^ and '^^'^^^^^^ 
in terms of their moments in (4.16), and by subtracting the results from (4.32). 
Following the same procedure as in Section 4.2.4 to express the remaining power 
corrections in terms of the leading shape function, we obtain 



(Mb - P+y \ 3</2 3 

These expressions remain unchanged when the shape function S is used instead of 
S. 



4.3.5 Weak annihilation contributions 

In the OPE calculation several contributions appear at third order in the power 
expansion: l/m^ corrections to the kinetic and chromo-magnetic operators, the 
Darwin and spin-orbit terms, and weak annihilation contributions. The Darwin 
and spin-orbit terms correspond to the forward S-meson matrix elements of (light) 
flavor- singlet operators [77]. The corresponding HQET parameters and can 
in principle be extracted from moments of inclusive B — > Xf.l'i' decay spectra. 
They are insensitive to the flavor of the spectator quark inside the B meson. The 
weak annihilation contribution, on the other hand, results from four-quark opera- 
tors with flavor non-singlet structure. Graphically, this contribution corresponds 
to a process in which the b and u quark annihilate into a . Weak annihilation 
terms come with a phase-space enhancement factor of IGvr^ and so are potentially 
more important than other power corrections of order Because of the flavor 
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dependence, these contributions can effect neutral and charged B mesons differ- 
ently [78] . One choice of basis for the corresponding four-quark operators is [55] 



{BlbL^^ULULrhlB) = ^-^^ Bi , {B\bRULULbR\B) = B2 , (4.34) 

where Jb is the 5-meson decay constant, and Bi are hadronic parameters. In the 
vacuum saturation approximation they are given by Bi — B2 — 1 for charged 
B mesons and Bi = B2 = for neutral ones. The total semileptonic rate is 
proportional to the difference {B2 — Bi), which implies that the weak annihilation 
contribution would vanish in this approximation. Currently, only rough estimates 
are available for the magnitude of the deviation of this difference from zero. The 
resulting effect on the total branching ratio is [61] 

SB(B X. « 3.9 (^) " (^) . (4.35) 

Again, we expect this effect to be different for charged and neutral B mesons. The 
most important feature of weak annihilation is that it is formally concentrated at 
the kinematical point where all the momentum of the heavy quark is transferred 
to the lepton pair [78]. At the parton level this implies that the corresponding 
contribution is proportional to S{q'^ — ml). It is therefore included in every cut 
that includes the endpoint, and its effect is independent of the specific form of 
the cut. 

We suggest two different strategies to control this effect. The first is to include 
it in the error estimate as a constant contribution proportional to the total rate. 
A recent study [79] puts a limit on this effect of ±1.8% on the total rate (at 68% 
confidence level) by analyzing CLEO data. The second one is to impose a cut 
9^ < 9max) thus avoiding the region where the weak annihilation contribution is 
concentrated. The maximal value of is {Mb — M^)^, but one must exclude a 
larger region of phase space, such that the excluded contribution to the decay rate 
at large (corresponding to a region near the origin in the (P_,P+) plane) can 
be reliably calculated. In our numerical analysis, we will study the effect of a cut 
q^ < {Mb — MuY , which satisfies this criterion. 

For completeness, we note that even after the weak annihilation contribution 
near maximum q^ has been removed, there could in principle exist other, flavor- 
specific contributions to the semileptonic decay amplitudes that are different for 
charged and neutral B mesons. The leading terms of this kind contribute at order 
l/m^ in the shape-function region and are parameterized by a set of four-quark 
subleading shape functions [66, 75, 76]. Model estimates of these contributions 
show that they are very small for all observables considered for an extraction of 
\Vub\ [76, 80]. If only flavor- averaged decay rates are measured, the effects of four- 
quark subleading shape functions can be absorbed entirely by a redefinition of the 
functions u{Cj) and v{uj) [66], without affecting the moment relations in (4.16). 



4.4 Shape-function parameterizations 



Hadronic-physics effects enter the description of inclusive decay rates via non- 
perturbative shape functions. Perturbation theory cannot tell us much about the 
local form of these functions, but moments of them are calculable provided that 
the domain of integration is much larger than Aqcd- Since the shape functions 
contain information about the internal structure of the B meson, knowledge of 
them relates directly to the determination of the 6-quark mass mb, the kinetic- 
energy parameter Ai, and in principle the matrix elements of higher-dimensional 
operators. Improved measurements of the shape of the B — > Xsj photon spectrum 
will therefore lead directly to a more precise determination of HQET parameters. 
This argument can be turned around to constrain the leading shape function using 
knowledge of and Ai from other physical processes such as a 6 — > c moment 
analysis [81]. We emphasize, however, that there are obviously infinitely many 
locally different functions that have identical first few moments. In this section we 
present a few functional forms that can be used to model the shape functions and 
to fit the current experimental data. 

To achieve stringent constraints on the leading shape function a precise defini- 
tion of the HQET parameters is required. It is a well-known fact that the pole- 
mass scheme introduces uncontrollable ambiguities. To avoid these uncertainties 
several short-distance definitions have been proposed, such as the MS scheme, the 
potential-subtraction scheme [40], the T(IS') scheme [82], the kinetic scheme [39], 
or the shape- function scheme [63] . While the decay rates are of course independent 
of the particular choice, it is advantageous to use a mass scheme that is designed for 
the physics problem at hand. In the case of inclusive B decays into light particles, 
this is the shape-function scheme. 

4.4.1 Models for the leading shape function 

Model-independent constraints on the shape function S{u;,iii) can be derived by 
analyzing moments defined with an upper limit of integration ujq, i.e. 



For practical applications, ujq should be taken of order the size of the window where 
the B Xs7 photon spectrum is experimentally accessible, ujq = Mb — 2E!^"^ with 
^min j-^ 1.8 GeV. These moments can be expanded in terms of matrix elements of 
local operators as long as cuq is large compared to Aqcd- In the shape-function 
scheme, HQET parameters are defined to all orders in perturbation theory through 
ratios of such moments, e.g. [63] 




(4.36) 
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Here, the factorization scale jif ^ Aqcd is related to the size of the integration 
domain via the implicit equation ujq = fi f + A{fi f , fii) . In practice fif is close to the 
intermediate scale /ij. At tree level, the relations between parameters in the shape- 
function scheme and the pole scheme are A(/x/,/ii) = Ap°'® and /^^(/x/,/^^) = — Ai. 
The corresponding relations at one- and two-loop order have been worked out in 
[63] and [83] , respectively. These relations allow us to obtain precise determinations 
of A(/x/,/ij) and fi'^{fif, fii) from other physical processes. 

For reference purposes, it is helpful to quote values for A and /x^ using only 
a single scale /x* instead of two independent scales Hf and /Xj. To one-loop order, 
these parameters can be related to those determined from the moments via 
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where we have neglected higher- dimensional operator matrix elements that are sup- 
pressed by inverse powers of /x/. A typical choice for the scale h* is 1-5 GeV, which 
we will use as the reference scale throughout this chapter. It will be convenient 
to connect the parameter A extracted from the first moment of the shape func- 
tion with a low-scale subtracted quark-mass definition referred to as the "shape- 
function" mass. Following [63], we define 



rribiiJLf, Hi) = Mb- A(/x/, Hi) ■ 



(4.39) 



The general procedure for modeling the leading shape function S{uj,Hi) from 
a given functional form F{uj) is as follows. The shape of -F(a') is assumed to be 
tunable so that it can be used to fit the B — > Xg^ photon spectrum. Only the 
norm of the shape function is fixed theoretically. Note that the moment relations 
(4.37) are insensitive to the norm, so that formulae for A and follow directly 
from the functional form of F{Cj). Examples of such formulae will be given below. 
We define moments Mjy (a)o) of F in analogy with (4.36). The first relation in 
(4.37) implies that for a given Cjq the factorization scale is 
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Now that /i/ is known, the norm is determined by requiring that the zeroth moment 



QO 



of the shape function is [63] 
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It follows that [Mo(a)o; /^j)/^d^'('^o)] -^('^) serves as a model of S{u), Hi) oi' S{Cj, Hi)- 
We now suggest three two-parameter models for the leading-order shape func- 
tion based on an exponential-type function F*^®^p\ a gaussian-type function F'^^'^^^^\ 
and hyperbolic-type function F^^'^^p). We use two parameters that can be tuned 
to fit the photon spectrum: a dimensionful quantity A which coincides with the 
position of the average (a)), and a positive number h which governs the behavior 
for small a). The functions we propose are 

F^-^\^:^.l>) = ^(f)"exp(-.,.4), 
Ffe--)(a;;A,6) = exp 

For convenience, we normalize these functions to unity. The parameters d{i) are 
determined by the choice A = (a)) . We find 
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where a) = X^fcLo(^ + '^)~'' generalized Riemann zeta function. An illus- 

tration of the different functional forms is given on the left-hand side in Figure 4.2. 

We show a plot with the choice 6 = 2, corresponding to a linear onset for small Co. 

For the first two models, analytic expressions for the HQET parameters A and 
Hi are available. Following the discussion above, we compute the moments on the 
interval [O,ci)o] and find for the exponential form F(^^P^(a); A, 6) 

A r(i + 6) -r(i + 6,^) 



h Tib) - Tib, ^) ' 

A2 r(2 + 6) -r(2 + 6,^; 
6^ r(6) -r(6,^) 



(4.44) 
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Figure 4.2: ToP: Different functional forms for the leading shape function. We 
show F(<=^P)(c2;,A,2) (solid), F'^^''''^^\Cj , k,2) (dotted), and F(^yp)(cO, A, 2) (dash- 
dotted) as functions of the ratio Cj/K. Bottom: The same functions with the 
parameters A and h tuned such that m6(/i*,/i*) = 4.61 GeV and /^^(/i,,,, /i*) = 0.2 
GeV^. See text for explanation. 
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where nj =Uq- A{iJ,f,iii). 

A similar calculation for the gaussian form F(s'iuss)(^. ^-j yjgyg 
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The corresponding relations for ^'^'^^^^(u); A, 6) must be obtained numerically. 

Ultimately the shape function should be fitted to the B Xg^ photon spec- 
trum, and the above equations then determine A and /x^. On the other hand, 
these formulae can be inverted to determine A and b from the current values of the 
HQET parameters. For example, if we adopt the values mfe(/i*,/i*) = 4.61 GeV 
and yU^,) = 0.20 GeV^ for the parameters in (4.38) at fi^, = 1.5 GeV, then 

we find the parameter pair A ^ 0.72 GeV, b ^ 3.95 for the exponential model, 
A fti 0.71 GeV, b ^ 2.36 for the gaussian model, and A ^ 0.73 GeV, b ^ 3.81 for 
the hyperbolic model. On the right-hand side of Figure 4.2 we show these three 
different functions plotted on the interval [0, cuq] over which the moment constraints 
are imposed. While the exponential (solid) and hyperbolic (dash-dotted) curves 
are barely distinguishable, the gaussian model has quite different characteristics. 
It is broader, steeper at the onset, faster to fall off, and the maximum is shifted 
toward larger a). 

An important comment is that, once a two-parameter ansatz is employed, the 
shape-function parameters (i.e., and yU^) can either be determined from a fit 
to the entire photon spectrum, or to the first two moments of the spectrum. Both 
methods are equivalent and should yield consistent results. If they do not, it would 
be necessary to refine the ansatz for the functional form of the shape function. 

In most applications shape functions are needed for arguments a) of order Aqcd- 
However, in some cases, like the ideal cut on hadronic invariant mass, uj is required 
to be as large as M/j, which is much larger than Aqcd- The large-cj behavior of the 
shape functions can be computed in a model-independent way using short-distance 
methods. For the leading shape function, one finds [63] 

» A«oD. /..) = 1 A i„ ^ + A + . . . . (4.46) 

TT u-Ay Hi J 

Note that this radiative tail is negative, implying that the shape function must go 
through zero somewhere near a) ~ few Aqcd- For practical purposes, we "glue" 
the above expression onto models of the non-perturbative shape function starting 
at LJ = A + iJ-i/ y/e ~ 1.6 GeV, where the tail piece vanishes. In this way we 
obtain a continuous shape-function model with the correct asymptotic behavior. 
We stress that for applications with a maximal P+ not larger than about 1.6 GeV 
the radiative tail of the shape function is never required. This includes all methods 
for extracting \Vub\ discussed later in this chapter, except for the case of a cut on 
hadronic invariant mass, Mx < Mq, if Mq is above 1.6 GeV. 



4.4.2 Models for subleading shape functions 



In the last section we have been guided by the fact that the B — > Xg^ photon 
spectrum is at leading power directly determined by the leading shape function. 
This helped in finding models that have roughly the same shape as the photon 
spectrum. At the subleading level considered here, however, no such guidance 
is provided to us. The available information is limited to the tree-level moment 
relations (4.16), stating that the norms of the subleading shape functions vanish 
while their first moments do not. In [66], two classes of models have been proposed, 
in which the subleading shape functions are "derived" from the leading shape 
function. A particularly simple choice is 



t(a)) = —X2S'{lj) , u{uj) 
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(4.47) 



Below, we will sometimes refer to this set of functions as the "default choice" . We 
choose the parameter — Ai in the expression for ■u(a)) (as well as in the expressions 
for the second-order hadronic power corrections) to coincide with the quantity 
li^{lif,lii) given in (4.44) and (4.45). However, for consistency with the tree- 
level moment relations, we identity the parameter A in (4.13) and (4.31) with the 
quantity A(//j,//j) evaluated in the limit where a;o ^ oo. This implies A = A for 
all three types of functions and ensures that the subleading shape functions have 
zero norm when integrated over < Co < oo. 

There are of course infinitely many possibilities to find models for subleading 
shape functions that are in accordance with (4.16). Any function with vanishing 
norm and first moment can be arbitrarily added to any model for a subleading 
shape function without violating the moment relations. Several such functions have 
been proposed in recent work on subleading shape functions, see e.g. [66, 74, 76, 80]. 
Specifically, we define the functions 



hi{uj) 

/l4(tD) 



^a-l ^-az f a{2-z)] , 



Nfll 2r(a 

M2 

— e 



Mo 



l-2az + 



z 
a^z^ 



M2 



— e 

7r-2 

TT^ 2x/7ra 



\-2z 



-2e-^ + 2z e-^^ Ei (z 



4 n-2 
8 



+ 



[l + z^ 



l-2z^|- 

TT 



z\nz+'-{l + z')-'^{l-z') 



(4.48) 



where z = uj/Qq, and the reference quantity Qq = O(Aqcd) depends on the type of 
function, namely f2o = A for hi and /i2, ^0 = | A for h^, and fio = ^ A for h^. The 
quantity a is a free parameter. The functions (4.48) have by construction vanishing 




Figure 4.3: Nine models for the subleading shape function u{uj) obtained by adding 
or subtracting one of the four functions hn{uj) to the defauh model in (4.47), shown 
as a thick line. See text for explanation. 



norm and first moment. Their second moments are given by the parameter M2, 
provided the normalization constants are chosen as = 1 for hi and /i2, and 

4 - vr 1 -n'^iA - tt) 1 

N = l- , \ - , N = l- ; / - (4.49) 
2(7r - 2) a 8(7r - 2) a ^ ^ 

for /z.3 and /14, respectively. The values for the parameters a and M2 should be 
chosen such that the following characteristics of subleading shape functions are re- 
spected: First, they are dimensionless functions, so that their values are naturally 
of 0(1) for ci) ~ Aqcd- Secondly, when integrated over a sufficiently large domain, 
their contributions are determined in terms of their first few moments. In partic- 
ular, this implies that for values of u) S> Aqcd the integrals over the subleading 
shape functions must approach zero. Taking these considerations into account, we 
use M2 = (0.3 GeV)^ in all cases and choose a = 3.5 for hi, a = 5 for /i2, and 
a = 10 for hs and h^. 

Given the four functions (4.48), we can construct several new models for the 
subleading shape functions t(a)), u(ct'), and {)(a)). For each function, we construct 
a set of 9 models by adding or subtracting any of the functions hn{oj) to the 
default choice in (4.47). Together, this method yields 9^ = 729 different sets 
{ti{u) , Uj{u) , Vk{u})} with i,j,k = 1, . . . , 9. This large collection of functions will 
be used to estimate the hadronic uncertainties in our predictions for partial decay 
rates. Note that for most of these sets we no longer have ti{u) = —Vk{uj), which 
was an "accidental" feature of the default model (4.47). The fact that the two 
functions have equal (but opposite in sign) first moments does not imply that 
their higher moments should also be related to each other. 

For the case of u{(2}) the resulting functions are shown in Figure 4.3, where 
we have used the exponential model (4.42) with parameters A = 0.72 GeV and 
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b — 3.95 for the leading shape function. In the region cu ~ Aqcd they differ 
dramatically from each other, while the large cu dependence is dominated by the 
moment relations (4.16). 



4.4.3 Illustrative studies 

We stressed several times that the calculation of the hadronic tensor is "optimized" 
for the shape-function region of large P_ and small P+, while it can smoothly be 
extended over the entire phase space. The notions "large P_" and "small P+" are 
to be understood as the sizes of integration domains for P_ and P+. Only when 
the differential distributions are integrated over a sufficiently large region in phase 
space, global quark-hadron duality ensures that the partonic description used in 
the present chapter matches the true, hadronic distributions with good accuracy. 
A more ambitious goal would be to calculate the differential decay rate point by 
point in the (P+, P_) plane. This can be done invoking local quark-hadron duality, 
as long as there is a sufficiently large number of hadronic final states contributing 
to the rate at any given point in phase space. 

It is instructive to integrate the triple differential decay rate (4.19) over the 
leptonic variable Pi in the range P+ < Pi < P-, which yields the exact formula 



Our theoretical prediction for the double differential decay rate (4.50) is shown 
on the left-hand side of Figure 4.4. We use the exponential model for the leading 
shape function with parameters mb{ii*, //*) = 4.61 GeV and iJ>^{iJ>*, //*) = 0.2 GeV^, 
as well as the default choice (4.47) for the subleading shape functions. For very 
small P_ values the rate turns negative (to the left of the gray line in the figure), 
signaling a breakdown of quark-hadron duality. It is reassuring that the only region 
where this happens is the "resonance region" , where the hadronic invariant mass 
is of order Aqcd, and local duality breaks down. 

Another useful quantity to consider is the differential P+ rate, which is obtained 
by integrating the double differential rate over P_ in the range P+ < P_ < Mb- 
The resulting P+ spectrum is shown on the right-hand side of Figure 4.4. In the 
plot we also disentangle the contributions from different orders in power counting. 

4.5 Predictions and error estimates for partial rates 

Before discussing predictions for partial B — > X„/~z/ rates for various kinematical 
cuts, let us recapitulate the ingredients of the calculation and general procedure. 
We have presented expressions for the triple differential decay rate, which can be 
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Figure 4.4: ToP: Theoretical prediction for the double differential decay rate. 
The light area represents a large decay rate. Black regions denote areas where 
the decay rate is close to zero. The dotted line is given by P+P_ = Ml, which 
means that charm background is located in the upper wedge. See text for further 
explanation. Bottom: The P+ spectrum extended to large values of P+. The thin 
solid line denotes the leading-power prediction, the dashed line depicts first-order 
power corrections, the dash-dotted line shows second-order power corrections, and 
the thick solid line is their sum. 
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organized in an expansion in inverse powers of {Mb — P+)- The leading-power con- 
tribution is given at next-to-leading order in renormalization-group improved per- 
turbation theory. At first subleading power two contributions arise. The first type 
involves subleading shape functions and is included at tree level, while the second 
type contributes perturbative corrections of order ag that come with the leading 
shape function. Further contributions enter at second subleading power and are 
again of the two types: perturbative corrections of order and non-perturbative 
structures at tree level. In summary, then, partial rates can be computed term by 
term in an expansion of the form 

r„ = r(°) + (rl^'^^^i) + r^^'^i^^) + (rl^''^^^) + r^^'^i^^) + .... (4.5i) 

The goal of this section is to test the convergence of this series expansion and 
to perform a thorough analysis of uncertainties. For the kinematical corrections 
the sum of all terms is known and given by the expressions in (4.26), while 
the first two terms in the series correspond to the expanded results in (4.28) and 
(4.29). We will find that in all cases of interest the first two terms give an excellent 
approximation to the exact result for T^^. 

For the purpose of illustration, we adopt the exponential model for the shape 
function and present numerical results for two sets of input parameters, which are 
biased by the results deduced from fits to -B ^ Xcl'v moments [83]. Specifically, 
we use mb{n^,n^) = 4.61 GeV, /i^(/i*,/i*) = 0.2 GeV^ (set 1) and mfo(/i.*, /i.*) = 
4.55 GeV, /x^(/x*, /x^,) = 0.3 GeV^ (set 2). The values of the 6-quark mass coincide 
with those obtained at two-loop and one-loop order in [83] (see also the discussion 
below), while the values of /i^ are close to the corresponding values in that refer- 
ence. As was mentioned before, in the future the leading shape function S{i2j,fii) 
should be extracted from a fit to the B — > Xgj photon spectrum, in which case 
the uncertainty in its shape becomes an experimental error, which can be system- 
atically reduced with improved data. In the process, the "theoretically preferred" 
parameter values used in the present work will be replaced with the "true" val- 
ues extracted directly from data. While this will change the central values for the 
partial rates, our estimates of the theoretical errors will only be affected marginally. 

The different sources of theoretical uncertainties are as follows: First, there 
are uncertainties associated with the functional forms of the subleading shape 
functions. To estimate them, we take the spread of results obtained when using 
the large set of different models described in Section 4.4.2, while the central value 
for a partial decay rate corresponds to the default model (4.47). Secondly, there 
are perturbative uncertainties associated with the choice of the matching scales 
/ih, IJ'i, and p,. Decay rates are formally independent of these scales, but a residual 
dependence remains because of the truncation of the perturbative series. Our error 
analysis is as follows: 

• The hard scale fih is of order m^. In perturbative logarithms the scale appears 
in the combination {yrrih/ fi^), see e.g. (4.25). To set a central value for fi^ we 
are guided by the average {y)mb. The leading term for the double differential 
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decay rate (PTu/dP+dy is proportional to 2y^(3 — 2y). It follows that the 
average y on the interval [0, 1] is 0.7. However, in some applications y is not 
integrated over the full domain. Also, there are large negative constants in 
the matching correction if^i in (4.25), whose effect can be ameliorated by 
lowering the scale further. In the error analysis we use the central value of 
lifi — mi,/2 2.3 GeV and vary the scale by a factor between l/-\/2 and -\/2. 
For the central value as{iJih) ~ 0.286. 

• The intermediate scale /li ~ -y/mfeAQCD serves as the renormalization point 
for the jet and shape functions. We fix this scale to /Xj = 1.5 GeV. Variations 
of fii would affect both the normalization and the functional form of the 
shape function, as determined by the solution to the renormalization-group 
equation for the shape function discussed in [63, 50]. In practice, effects on 
the shape are irrelevant because the shape function is fitted to data. The only 
place where the intermediate scale has a direct impact on the extraction of 
\Vub\ is through the normalization of the shape function (4.41). In the analysis 
we therefore estimate the uncertainty by assigning the value ±( "''^^''' )^ as a 
relative error, where Q;s(A*i) ~ 0.354. 

• The scale jl appears as the argument of ag in the perturbative contributions 
rj^™. We vary /2 from /jj/ y/2 to V^/ii with the central value p, — jii — l.b GeV. 

These three errors are added in quadrature and assigned as the total perturbative 
uncertainty. Finally, we need to estimate the effects from higher- dimensional op- 
erators at third and higher-order in power counting. If the considered cut includes 
the region of phase space near the origin (P+ ~ P_ ~ Aqcd), then the domi- 
nant such contributions are weak annihilation effects, which we have discussed in 
Section 4.3.5. From the analysis in [79] one can derive a bound on the weak an- 
nihilation contribution that is ±1.8% of the total decay rate, for which we take 
Vu ~ 70 iKipps^"*^ (see below). The resulting uncertainty dV"^^ = ±1.3 |K(fePps^^ 
affects all partial rates which include the region near the origin in the (P+, P_) 
plane. The uncertainty from weak annihilation can be avoided by imposing a 
cut < (see Section 4.5.6). For all observables considered in the present 
work, other power corrections of order 1/m^ can be safely neglected. This can 
be seen by multiplying the contributions from second-order hadronic power cor- 
rections to the various decay rates (called additional suppression 
factor Aqcd /rub ~ 0.1. 

The following subsection contains a discussion of the total decay rate. In the 
remainder of this section we then present predictions for a variety of kinematical 
cuts designed to eliminate (or reduce) the charm background. These partial rates 
can be computed either numerically or, in many cases, semi-analytically. In Ap- 
pendix C we discuss how to perform the integrations over the kinematical variables 
P; and P_ analytically. 



4.5.1 Total decay rate 



Before presenting our predictions for the various partial decay rates, it is useful 
to have an expression for the total B — > X^l^P decay rate expressed in terms of 
the heavy-quark parameters defined in the shape-function scheme. We start from 
the exact two- loop expression for the total rate derived in [84], add the second- 
order hadronic power corrections, which are known at tree level [13, 1], and finally 
convert the parameters rrif, and Ai from the pole scheme to the shape-function 
scheme. The relevant replacements at two-loop order can be taken from [83] and 
read 
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where here and from now on rrib = "f^bil^*-! 1^*) ^-nd //^ = //^(/i*,//*) are defined 
in the shape-function scheme. At a reference scale = 1.5 GeV the values of 
these parameters have been determined to be rrib = (4.61 ± 0.08) GeV and = 
(0.15 ± 0.07) GeV^ [83],^ where we account for the small l/rrib correction to the 
relation for the pole mass in the above formula (corresponding to a shift of about 
—0.02 GeV in rrib), which was not included in that paper. 
The resulting expression for the total decay rate is 
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We observe that for ^ 1.5 GeV and // = 0{mb), the perturbative expansion coef- 
ficients are strongly reduced compared to their values in the pole scheme (—0.768 
and —2.158, respectively), indicating a vastly improved convergence of the per- 
turbative expansion. For rrib = 11 = 4.61 GeV, and fi^ = 0.15 GeV^ we obtain 
for the one-loop, two-loop, and power corrections inside the brackets in (4.53): 
{1 — 0.017 — 0.030 — 0.004}. All of these are very small corrections to the leading 
term. 

Including the uncertainties in the values of rrib and quoted above, and varying 
the rcnormalization scale /j, between rrib and mb/2 (with a central value of mb/ a/2). 



^Thc values obtained from a one-loop analysis are rrib — (4.55 ± 0.08) GeV and 
= (0.34 + 0.07) GeV^. 
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Table 4.1: Partial decay rate Tu{Eq) for a cut on charged-lepton energy Ei > Eq 
in the S-meson rest frame, given in units of |l^fepps~^. Predictions are based 
on the shape-function parameter values mi, — 4.61 GeV, /.t^ = 0.2 GeV^ (top) and 
rrib = 4.55 GeV, nl = 0.3 GeV^ (bottom). 



Eo [GeV] 


Mean 


Subl. SF 


Pert. 


Total 


1.9 
2.0 
2.1 
2.2 
2.3 
2.4 


24.79 
18.92 

13.07 
7.59 
3.12 
0.42 


±0.53 
±0.60 

±0.71 
±0.81 
±0.89 
±1.05 


+i.yu 

-1.66 
+1.35 
-1.20 
+0.82 
-0.75 
+0.38 
-0.34 
+0.15 
-0.16 
+0.16 
-0.22 


+2.34 
-2.15 
+1.95 
-1.84 
+1.66 
-1.63 
+1.55 
-1.54 
+1.55 
-1.55 
+1.65 
-1.65 


1.9 
2.0 

2.1 
2.2 
2.3 
2.4 


21.10 
15.83 

10.73 
6.12 
2.47 
0.29 


±0.53 
±0.60 

±0.68 
±0.74 
±0.84 
±0.99 


+1.57 
-1.35 
+1.08 
-0.94 
+0.64 
-0.55 
+0.31 
-0.23 
+0.17 
-0.22 
+0.18 
-0.24 


+2.08 
-1.92 

+1.77 
-1.68 
+1.57 
-1.54 
+1.50 
-1.48 
+1.53 
-1.53 
+1.61 
-1.62 



we get 

1^ - 68.0 1L'K]± 0.7 

I ^ub I P° 

= (68.0±0.7[;.^]^-[;.])(^^^)"\ (4.54) 

Here and below, we quote values for decay rates in units of |Ki6pps^^- To convert 
these results to partial branching fractions the numbers need to be multiplied by 
the average 5-meson lifetime. Without including the two-loop corrections, the 
central value in the above estimate increases to 70.6. For comparison, with the 
same set of input parameters our new approach based on (4.19) predicts a total 
decay rate of r„ = (71.4^5q ± 0.5) |Ki6pps~^, where the first error accounts for 
perturbative uncertainties while the second one refers to the modeling of subleading 
shape functions (to which there is essentially no sensitivity at all in the total rate). 
The fact that this is in excellent agreement with the direct calculation using (4.53) 
supports the notion that the formalism developed in this chapter can be used to 
describe arbitrary B Xul'v decay distributions, both in the shape-function 
region and in the OPE region of phase space. 

4.5.2 Cut on charged-lepton energy 

Traditionally, the most common variable to discriminate against the charm back- 
ground is the charged-lepton energy Ei. As long as one requires that Ei is bigger 
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Table 4.2: Same as Table 4.1, but for the partial decay rate (^'o) for a cut 

on lepton energy Ei > Eq in the T(4S') rest frame. 



E, [GeV] 


Mean 


Sul)l. SF 


Pert. 


Total 


1.9 
2.0 
2.1 
2.2 
2.3 
2.4 


24.82 

19.00 

13.25 

7.99 

3.83 

1.31 


±0.54 
±0.61 
±0.71 
±0.78 
±0.86 
±0.99 


-1.66 
+1.37 
-1.21 
+0.85 
-0.76 
+0.42 
-0.37 
+0.18 
-0.13 
+0.10 
-0.14 


+2.35 
-2.15 
+1.96 
-1.85 
+1.68 
-1.63 
+1.54 
-1.53 
+1.54 
-1.53 
+1.61 
-1.61 


1.9 
2.0 
2.1 
2.2 
2.3 
2.4 


21.16 

15.94 

10.94 

6.49 

3.05 

0.98 


±0.54 
±0.60 
±0.68 
±0.74 
±0.84 
±0.92 


+1.58 
-1.35 
+1.10 
-0.95 
+0.66 
-0.57 
+0.34 
-0.26 
+0.17 
-0.18 
+0.13 
-0.18 


+2.09 
-1.93 
+1.78 
-1.69 
+1.58 
-1.54 
+1.50 
-1.48 
+1.53 
-1.53 
+1.56 
-1.57 



than (M^ — M'^)/2Mb ~ 2.31 GeV, the final hadronic state cannot have an invari- 
ant mass larger than M^. For this ideal cut, and using the default set of subleading 
shape functions, we find 

I /'Fkm(l) -phadr(l)\ /-pkm(2) . -phadr(2)\ 

i M + 1,1 u + i « j + (,i it + i « ) gg\ 

= [ 6.810 ±(0.444 - 3.967) ±(0.042 - 0.555)] iKbpps"^ 

The corrections from subleading shape functions are quite sizable, in accordance 
with the findings in [52, 73, 74]. Note that the sum Fu"*'^'' ± fJI"^*-^^ = 0.486 is an 
excellent approximation to the exact result Fj^'" = 0.482 (all in units of |V^Mfepps~^) 
obtained using (4.26), indicating that the expansion of the kinematical power cor- 
rections is converging rapidly. The same will be true for all other observables 
considered below. 

In practice, the cut on Ei can be relaxed to some extent because the back- 
ground is well understood, thereby increasing the efficiency and reducing the im- 
pact of theoretical uncertainties. Our findings for different values of the cut Eq 
are summarized in Table 4.1. Here and below, the columns have the following 
meaning: "Mean" denotes the prediction for the partial decay rate, "Subl. SF" the 
uncertainty from subleading shape functions, and "Pert." the total perturbative 
uncertainty. In the column "Total" we add the stated errors plus the uncertainty 
from weak annihilation in quadrature. 

Experiments often do not measure the partial rates in the S-meson rest frame, 
but in the rest frame of the T(4S') resonance produced in e+e" collisions. Boosting 
to the T(45') frame with (3 — v/c^ 0.064 has a small effect on the spectrum and 
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rates. The exact formula for this boost is [32] 
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Mb/2 j-nyj 
(3-Eo "-^ 



(3+ -max ( P-,^ 



(4.56) 



where /3± = VlTp/^/T^, and the factor 7 = 1/^1- 



hand side takes the time dilation of the B-meson lifetime 



1.002 on the left- 

•jtb into account. 

(In other words, branching fractions are Lorentz invariant.) The above formula 
can be accurately approximated by the first term in an expansion in which 
yields [32] 



d 



1 c/ri^) 



dE E dE 



+ 0(/3^), 



(4.57) 



E=Eo 



as long as Eq is not too close to the kinematical endpoint (i.e., Eq < P_Mb/2 a; 
2.47 GeV). The numerical results for the partial decay rate 'yT^\Eo) in the rest 
frame of the T(4S') resonance are given in Table 4.2. 



4.5.3 Cut on hadronic P+ 

Cutting on P+ samples the same hadronic phase space as a cut on the charged- 
lepton energy, but with much better efficiency [63, 85]. The phase space P+ < Ap 
with the ideal separator Ap = M|,/Mp ^ 0.66 GeV contains well over half of all 
B Xu events. Here we find with the default settings 

p(0) _|_ ^pkin(l) _|_ phadr(l)-^ _|_ ^pkin(2) _|_ phadr(2)^ 

= [ 53.225 +(4.646 - 11.862) +(0.328 - 0.227)] |K&Pps-^ 

(4.58) 

Wc sec a much better convergence of the power series than in the case of a cut on 
the charged-lepton energy, namely 53.225 — 7.216 — 0.100 when grouping the above 
numbers according to their power counting. Once again, the sum F^^^^^ + rj^™^^-' = 
4.973 is very close to the full kinematical correction T^"^ — 4.959 (in units of 

iKbpps-^). 

Often times it is required to impose an additional cut on the charged-lepton 
energy, as leptons that are too soft are difficult to detect. In Table 4.3 we list 
results for both Ei > and Ei > 1.0 GeV. For the ideal cut we find that the 
prediction is quite precise, as the total theoretical uncertainty is only about 6.8%. 
For comparison, the ideal cut for the lepton energy is uncertain by about 50%, but 
rapidly improving as the energy cut is relaxed. 



4.5.4 Cut on hadronic invariant mass and 

The most efficient separator for the discrimination of B XJ'u events is a cut 
on the invariant mass Mx of the hadronic final state, Mx < [86, 87]. It has 



Table 4.3: Partial decay rate r„(Ap,£'o) for a cut on the hadronic variable P+ < 
Ap and lepton energy Ei > Eq, given in units of |Kibpps^^. Predictions are based 
on the shape-function parameter values rrib = 4.61 GeV, ii^ = 0.2 GeV^ (top) and 
rub = 4.55 GeV, nl = 0.3 GeV^ (bottom). 



Ap [GeV] 


Eo [GeV] 


Mean 


Subl. SF 


Pert. 


Total 


0.70 
0.65 


0.0 
0.0 


48.90 
45.34 


±1.15 
±1.46 


+2.83 
-2.65 
+2.55 
-2.41 


+3.30 
-3.15 
+3.20 
-3.09 


0.60 


0.0 


41.34 


±1.76 


+2.26 
-2.15 


+3.13 
-3.05 


0.55 


0.0 


36.91 


±2.01 


+1.95 
-1.87 


+3.08 
-3.02 


0.50 


0.0 


32.09 


±2.34 


+1.64 
-1.58 


+3.12 

-3.09 


0.70 


1.0 


43.36 


±1.02 


+2.54 
-2.39 


+3.01 
-2.88 


0.65 


1.0 


40.18 


±1.30 


+2.28 
-2.16 


+2.92 
-2.82 


0.60 


1.0 


36.59 


±1.59 


+2.01 
-1.92 


+2.86 
-2.80 


0.55 


1.0 


32.61 


±1.86 


+1.73 

-1.67 


+2.84 
-2.80 


0.50 


1.0 


28.29 


±2.19 


+ 1,41 
-1.10 


+2.91 
-2.89 


0.70 


0.0 


39.95 


±1.19 


+2.18 
-2.06 


+2.79 
-2.70 


0.65 


0.0 


36.94 


±1.42 


+1.95 
-1.86 


+2.72 
-2.66 


0.60 


0.0 


33.67 


±1.65 


+1.71 
-1.65 


+2.69 
-2.65 


0.55 


0.0 


30.15 


±1.88 


+1.47 
-1.43 


+2.70 
-2.68 


0.50 


0.0 


26.40 


±2.09 


+1.22 
-1.21 


+2.73 
-2.72 


0.70 


1.0 


35.42 


±1.13 


+1.95 

-1.85 


+2.59 

-2.51 


0.65 


1.0 


32.73 


±1.34 


+1.74 
-1.66 


+2.53 
-2.48 


0.60 


1.0 


29.81 


±1.55 


+1.52 
-1.47 


+2.51 
-2.48 


0.55 


1.0 


26.65 


±1.76 


+1.29 
-1.27 


+2.52 
-2.51 


0.50 


1.0 


23.29 


±1.95 


+1.07 
-1.06 


+2.56 
-2.55 



Table 4.4: Partial decay rate r„(Mo, q^) for combined cuts Mx < Mq on hadronic 
invariant mass, > q^ on leptonic invariant mass, given in units of iKttpps^^. 
Predictions are based on the shape-function parameter vahies rrih = 4.61 GeV, 
fil = 0.2 GeV^ (top) and nib = 4.55 GeV, i^l = 0.3 GeV^ (bottom). 



Mo [GeV] 


ql [GeV^] 


Eo [GeV] 


Mean 


Subl. SF 


Pert. 


Total 


Md 


0.0 


0.0 


59.30 


±0.36 


+A:2'2 
-3.73 


+4.42 
-3.96 


1.70 


0.0 


0.0 


53.13 


±0.73 


+3.67 
-3.31 


+3.95 
-3.61 


1.55 


0.0 


0.0 


45.72 


±1.16 


+3.11 
-2.84 


+3.55 
-3.32 


Md 


6.0 


0.0 


34.37 


±0.37 


+2.97 
-2.58 


+3.25 
-2.89 


1.70 


8.0 


0.0 


24.80 


±0.36 


+2.24 
-1.98 


+2.59 
-2.37 


Md 


{Mb - Md? 


0.0 


12.55 


±0.49 


+1.41 
-1.24 


+1.95 
-1.83 


Md 
1.70 


0.0 
0.0 


0.0 
0.0 


50.08 
44.20 


±0.54 
±0.86 


+3.52 

-3.11 
+2.98 
-2.69 


+3.78 
-3.40 
+3.35 
-3.09 


1.55 


0.0 


0.0 


37.76 


±1.22 


+2.46 
-2.26 


+3.03 
-2.86 


Md 


6.0 


0.0 


29.42 


±0.35 


+2.50 
-2.16 


+2.82 
-2.52 


1.70 


8.0 


0.0 


20.87 


±0.39 


+1.84 
-1.61 


+2.26 
-2.08 


Md 


{Mb - MdT 


0.0 


10.49 


±0.48 


+1.16 
-1.00 


+1.76 
-1.68 



Table 4.5: Partial decay rate r„(Mo, -^o) for combined cuts Mx < Mq on 
hadronic invariant mass, q"^ > q^ on leptonic invariant mass, and Ei > Eq on 
charged- lepton energy, given in units of iT^bl^ps"^- Predictions are based on the 
shape-function parameter values rrib ~ 4.61 GeV, /i^ — 0.2 GeV^ (top) and = 
4.55 GeV, fxl = 0.3 GeV^ (bottom). 



Mo [GeV] 


ql [GeV^] 


Eo [GeV] 


Mean 


Subl. SF 


Pert. 


Total 


Md 
1.70 


0.0 

0.0 


1.0 

1.0 


53.49 

48.25 


±0.36 

±0.63 


+3.91 
-3.45 
+3.42 
-3.08 


+4.13 
-3.69 
+3.70 
-3.38 


1.55 


0.0 


1.0 


41.81 


±1.03 


+2.91 
-2.66 


+3.34 
-3.12 


Md 
1.70 


6.0 
8.0 


1.0 
1.0 


33.88 
24.74 


±0.37 
±0.36 


+2.94 
-2.55 
+2.23 
-1.97 


+3.22 
-2.87 
+2.59 
-2.37 


Md 


{Mb - MdY 


1.0 


12.55 


±0.49 


+1.41 
-1.24 


+1.95 
-1.83 


Md 


0.0 


1.0 


45.29 


±0.50 


+3.27 

-2.88 


+3.54 
-3.18 


1.70 


0.0 


1.0 


40.22 


±0.77 


+2.78 
-2.50 


+3.15 
-2.90 


1.55 


0.0 


1.0 


34.55 


±1.09 


+2.31 
-2.11 


+2.85 
-2.69 


Md 
1.70 


6.0 
8.0 


1.0 
1.0 


28.99 
20.82 


±0.34 
±0.39 


+2.48 
-2.13 
+1.83 
-1.60 


+2.80 
-2.50 
+2.26 
-2.08 


Md 


{Mb - MdT 


1.0 


10.49 


±0.48 


+1.16 
-1.00 


+1.78 
-1.68 



also been argued [88] that a cut on (f can reduce the shape-function sensitivity, 
since it avoids the colhncar region in phase space where P ^ P^. In order to 
optimize signal efficiency and theoretical uncertainties, it was suggested in [89] to 
combine a cut with a cut on hadronic invariant mass. 

The theoretical predictions obtained in [88, 89] were based on a conventional 
OPE calculation, which was assumed to be valid for these cuts. The assessment 
of the shape-function sensitivity was based on convolving the tree-level decay rate 
with a "tree-level shape function" , for which two models (a realistic model similar 
to the ones considered here, and an unrealistic 5-function model) were employed. 
The shape-function sensitivity was then inferred from the comparison of the re- 
sults obtained with the two models. The sensitivity to subleading shape functions 
was not considered, since it was assumed to be very small. Since our formal- 
ism smoothly interpolates between the "shape-function" and "OPE" regions, and 
since we include radiative corrections as well as power corrections as far as they 
are known, we can estimate the sensitivity of a combined Mx-(f cut to the leading 
and subleading shape functions much more accurately. Contrary to [89] , we do not 
find a significant reduction of the shape-function sensitivity when adding the (f 
cut to a cut on hadronic invariant mass. 

In Tables 4.4 and 4.5 we give results for typical cuts on Mx and g^, with and 
without including an additional cut on charged-lepton energy. Let us study the 
contributions for the optimal cut Mx < Md in detail. We find with the default 
settings 

r^f) + {V^^''^'' + V^''^^) + (r!l™^^^ + phadr(2)^ 

= [ 58.541 +(8.027 - 9.048) +(2.100 - 0.318)] iK^pps-i. 

(4.59) 

Note the almost perfect (accidental) cancellation of the two terms at order l/m(,. 
The resulting power series, 58.541 — 1.022 + 1.782, again exhibits good convergence. 
As previously, the sum V^^^^^ + P^™*^^' = 10.127 is a good approximation to the 
exact value P|^'° = 9.753 (in units of |Kj6pps~^)- The analogous analysis for a 
combined cut Mx < 1.7 GeV and q'^ > 8.0 GeV^ reads 

-p(O) _|_ ^pkin(l) -phadr(l)-^ ^-pkin(2) -phadr(2)-j 

= [ 25.880 +(4.049 - 6.358) +(1.399 - 0.171)] |146pps"S 

(4.60) 

which means that the power series is 25.880 — 2.309 + 1.228. Here we have P^"^^^ + 



4.5.5 Cut on s™^ and Ei 

In [90], the BaBar collaboration employed a cut on both Ei > Eq and a new 
kinematical variable s^^^ < Sq, where the definition for involves both hadronic 
and leptonic variables. In the S-meson rest frame, it is 

s]r ^Ml + q^- 2Mb + . (4.61) 



Table 4.6: Partial decay rate r„(so, Eq) for combined cuts s^^^^ < sq and Ei > Eq, 
given in units of iKi^P ps~^. Predictions are based on the shape- function parameter 
values rub = 4.61 GeV, nl = 0.2 GeV^ (top) and = 4.55 GeV, = 0.3 GeV^ 
(bottom) . 



So [GeV] 


Eo [GeV] 


Mean 


Subl. SF 


Pert. 


Total 


3.5 
3.5 

3.5 
3.5 


1.8 
1.9 

2.0 
2.1 


17.39 
15.86 

13.70 
10.78 


±0.62 
±0.63 

±0.66 
±0.73 


+ 1.54 
-1.36 
+1.33 
-1.18 
+1.05 
-0.94 
+0.71 
-0.64 


+2.08 
-1.96 
+1.94 
-1.84 
+1.77 
-1.71 
+1.62 
-1.59 


3.5 
3.5 
3.5 

3.5 


1.8 
1.9 
2.0 

2.1 


14.57 
13.18 
11.28 

8.77 


±0.60 
±0.61 
±0.64 

±0.69 


+1.25 
-1.09 
+1.06 
-0.92 
+0.82 
-0.71 
+0,51 
-0.46 


+1.87 
-1.77 
+1.76 
-1.68 
+1.63 
-1.58 
+ 1.54 
-1.51 



Rewriting the phase space of this cut in the variables P+, P_, Pi, we find 

0<P+ < mm{MB-2Eo,y/^), 
P+<P- < min(^-^,MB^, 

P+<Pi < mm{MB-2Eo,P-), (4.62) 

where it is understood that if = (M^ - P+){Mb - P-) < {Mb - v^)^ then 
the interval P™'"^ < Pi < Pl"^ must be excluded from the Pi integration. Here 

pmax/min^p^ p -J _ 




+ 2Mb ) ^'-'^^ 

A summary of our findings is given in Table 4.6. When compared to the 
pure chargcd-lcpton energy cut in Table 4.1, the additional cut on s^^^ eliminates 
roughly another 20-30% of events. However, the hope is that this cut also reduces 
the sensitivity to the leading shape function, which we expect to be sizable for the 
pure El cut. The uncertainty from subleading shape functions, however, is almost 
unaffected by the cut. 

4.5.6 Eliminating weak annihilation contributions 

In Section 4.3.5 we have argued that a cut on high g^, i.e., < q^, will eliminate 
the effect of weak annihilation and remove the uncertainty associated with this 
contribution. The cutoff should be small enough to exclude the region around 



Table 4.7: Examples of partial decay rates with a cut on < (Mb — Md)^ imposed 
to eliminate the weak annihilation contribution. We consider an additional cut on 
the hadronic variable P+ < Ap (top) , or on the hadronic invariant mass Mx < Mq 
(bottom). As before, decay rates are given in units of |Vubpps~^. Predictions are 
based on the shape-function parameters rrib — 4.61 GeV and /i^ — 0.2 GeV^. 



Ap [GeV] 


Mean 


Subl. SF 


Pert. 


Total 


0.70 

0.65 
0.60 
0.55 
0.50 


39.96 
37.18 
34.05 
30.61 
26.86 


±1.27 

±1.50 
±1.71 
±1.89 
±1.97 


+2.16 
-2.01 
+1.99 
-1.85 
+1.82 
-1.69 
+1.63 
-1.52 
+1.44 
-1.33 


+2.51 
-2.38 
+2.49 
-2.38 
+2.50 
-2.41 
+2.49 
-2.42 
+2.44 
-2.38 


Mo [GeV] 


Mean 


Subl. SF 


Pert. 


Total 


Mn 
1.70 
1.55 


46.75 
40.70 
33.69 


±0.65 
±1.12 
±1.56 


+2.82 
-2.50 
+2.32 
-2.11 
+1.88 
-1.73 


+2.89 
-2.58 
+2.58 
-2.39 
+2.44 
-2.32 



= ml, where this contribution is concentrated. It is instructive to assess the 
"cost" of such an additional cut in terms of the loss of efficiency and, more impor- 
tantly, the behavior of the remaining uncertainties. In order to do this, we combine 
the cut < {Mb — MdY with either a cut on P+ or on Mx- While this particular 
choice for q^ still leaves some room to improve the efficiency by increasing q^, it is 
not desirable to raise the cut much further, since this would threaten the validity 
of quark-hadron duality. 

The results are summarized in Table 4.7 and can be compared to the previous 
"pure" P+ and Mx cuts in Tables 4.3, 4.4, and 4.5. As an example, let us consider 
the case < 0.65 GeV, which is close to the charm threshold. Without the 
additional q^ cut we found that the total theoretical uncertainty (including the 
weak annihilation error) is When cutting in addition on q^ < (Mb — MdY, 

the efficiency decreases by about 20% as expected. However, due to the absence 
of the weak annihilation uncertainty, the overall uncertainty decreases to 11q'1%- 
Therefore both strategies result in comparable relative uncertainties, with a slight 
favor for imposing the additional cut from the theoretical point of view. 

While the small reduction of theoretical errors hardly seems worth the effort 
of imposing the q^ cut, performing an analysis of the type outlined here and com- 
paring its results with those obtained without the additional cut may help to 
corroborate the expectation that the weak annihilation contribution is indeed not 
much larger than what has been found in [79]. 
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Table 4.8: Values of the exponent a{mb) for different kinematical cuts. The pa- 
rameter ji^ — 0.2 GeV^ is kept fixed. Also quoted is the sensitivity of the partial 
decay rates to the functional form of the shape functions. See text for explanation. 





nib [GeV] 


4.50 


4.55 


4.60 


4.65 


4.70 


Mx < Md 


a 

Functional Form 


9.5 
1.4% 


8.8 
1.1% 


8.2 
0.8% 


7.7 
0.5% 


7.3 
0.4% 


Mx < 1.7 GeV 


a 

Functional Form 


12.5 

2.9% 


11.5 

2.6% 


10.5 

2.2% 


9.7 

1.9% 


8.9 

1.6% 


Mx < 1.7 GeV 
q^>8 GeV2 


a 

Functional Form 


10.3 
2.0% 


9.8 
1.7% 


9.3 
1.5% 


9.0 
1.4% 


8.7 
1.4% 


g2 > (Mb - Mof 


a 

Functional Form 


11.4 
5.0% 


11.1 
4.4% 


10.9 
4.0% 


10.8 
3.6% 


10.6 
3.2% 


P+ < mi/Mb 


a 

Functional Form 


16.7 

5.3% 


15.0 
4.8% 


13.6 
4.4% 


12.2 
4.0% 


11.1 
3.6% 


El > 2.2 GeV 


a 

Functional Form 


22.6 
16.2% 


21.0 
13.1% 


19.7 
11.0% 


18.5 
9.3% 


17.4 
7.9% 



4.5.7 Dependence on rrib and shape-function sensitivity 

Non-perturbative hadronic physics enters in our approach via the form of the 
leading and subleading shape functions. The strongest sensitivity by far is to the 
first moment of the leading shape function, which determines the HQET parameter 
A and with it the 6-quark mass. Given that the value of = m5(/i^,;U*) can 
be determined with good precision from other sources (such as moments of the 
leptonic or hadronic invariant mass spectra in B ^ Xc l^v decays), it is instructive 
to disentangle this dependence from the sensitivity to higher moments or, more 
generally, to the functional form of the shape functions for fixed m^. 
To explore the dependence on mi, we define the exponent 

which means that ~ Table 4.8 shows the values of this exponent over 

a wide range of values of for a variety of experimental cuts. To estimate the 
sensitivity to the functional form we scan over a large set of models for the sublead- 
ing shape functions, and we also study the difference between the results obtained 
using the exponential or the gaussian ansatz for the leading shape function. The 
corresponding variations are added in quadrature and given as a relative change 
in the corresponding partial decay rates (labeled "Functional Form"). In all cases, 
/i^ = 0.2 GeV^ is kept fixed. Because we restrict ourselves to only two functional 
forms for the leading shape function in this study, the resulting sensitivities should 
be interpreted with caution. 
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The entries in the table are listed in roughly the order of increasing sensitivity 
to rriiy and to the functional form of the shape functions, with the hadronic invari- 
ant mass cut showing the least sensitivity and the lepton energy cut exhibiting 
the largest one. To some extent this reflects the different efficiencies (or "inclu- 
siveness") of the various cuts. It is reassuring that a ~ 10 for the pure cut, in 
accordance with the findings of [91, 92]. Perhaps somewhat surprisingly, for this 
cut a substantial sensitivity to shape-function effects remains even for fixed rrih 
and /i^. It is well known that the partial rate with a cut > {Mb — M£,Y can 
be calculated using a local OPE in powers of Aqcd/'^c [88, 91], thereby avoiding 
the notion of shape-function sensitivity. Differences between the functional forms 
of the shape functions in our approach correspond to effects that are formally of 
order 1/m^ and higher. It is not unreasonable that these effects should be of order 
3-5%. 

We also checked that for much more relaxed cuts the value of a(mf,) tends to 4.8, 
as stated in (4.54). For example, for a cut P+ < Ap we find (with = 4.61 GeV 
and nl = 0.2 GeV^): 



Ap [GeV] 


0.6 0.8 1.0 1.2 1.6 2.0 3.0 Mp 


a 


15.4 9.8 7.0 5.8 5.1 5.0 4.9 4.8 



4.6 Conclusions 

A high-precision measurement of the parameters of the unitarity triangle is an 
ongoing quest, which necessitates the close cooperation of theory and experiment. 
The determination of \Vuh\ from inclusive B — > X„Z~P decay requires the mea- 
surement of partial decay rates with kinematical cuts that eliminate the large 
background from B X(.l~u decay, as well as theoretical predictions for such 
quantities. To this end, it is desirable to have a theoretical description of the triple 
differential decay rate, which can be used for predicting arbitrary partial rates 
obtained after integrating over certain regions of phase space. One problem in 
providing such a description is that the power-counting rules of the heavy-quark 
expansion are different in different kinematical domains. In this chapter we have 
overcome this difficulty. 

In the shape-function region, our results are in agreement with QCD factoriza- 
tion theorems, and perturbative effects have been separated from non-perturbative 
shape functions. When the allowed phase space extends over a large domain, our 
results smoothly reduce to the expressions obtained from the local operator prod- 
uct expansion. We have presented a formalism in which event distributions and 
partial decay rates are expressed without explicit reference to partonic quantities 
such as the 6-quark mass. The sensitivity to such hadronic parameters enters indi- 
rectly, via the moments of shape functions. The most important non-perturbative 
object, namely the leading-order shape function, can be extracted from the photon 
spectrum in B ^ Xg'-f decay. This is analogous to extractions of parton distri- 
bution functions from fits to data on deep inelastic scattering. In this way, the 
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dominant uncertainty from our ignorance about bound-state effects in the B me- 
son is turned into an experimental uncertainty, which will reduce with increasing 
accuracy of the experimental data on the photon spectrum. Residual hadronic 
uncertainties are power suppressed in the heavy-quark expansion. 

One goal of this chapter was to present a detailed framework in which this 
program can be carried out. We have given formulae that can be readily used 
for the construction of an event generator, as well as to estimate the remaining 
theoretical uncertainties in a robust and automated fashion. 

In practice the leading shape function needs to be parameterized. We have 
suggested three different functional forms, which can be used to fit the data of 
the B — > Xg^ photon spectrum. Once the data is accurately described by a 
choice of the shape functions, this function can be used in the predictions for 
partial B Xul~iy rates and spectra. Subleading shape functions give rise to 
theoretical uncertainties starting at the level of l/m^ power corrections. We have 
estimated these uncertainties using a large set of models, each of which obeys the 
known tree-level moment relations, but which are very different in their functional 
form. A second error estimate is determined by the residual renormalization- 
scale dependence. We also considered uncertainties from weak annihilation effects, 
which in principle can be avoided by cutting away the region of phase space in 
which they contribute. We have suggested a cut on high leptonic invariant mass, 
which accomplishes just that. 

The second half of this chapter contains detailed numerical predictions for a 
variety of partial B — > X„ l~u decay rates with different kinematical cuts, including 
cuts on the charged-lepton energy (both in the rest frame of the B meson and of the 
T(45') resonance), on the hadronic quantity P+ = Ex — on Mx, on g^, and on 
various combinations of these variables. Along with our predictions for the rates we 
have presented a complete analysis of theoretical uncertainties. Once the data on 
the B photon spectrum are sufficiently precise to accurately determine the 

leading-order shape function, a determination of \ Vub\ with theoretical uncertainties 
at the 5-10% level now seems feasible. 
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Chapter 5 

Applications: Weight Function 

5.1 Introduction 

The kincmatical variable P+ is common to semileptonic and radiative charmless 
inclusive B decays. In B ^ X^'-f decays P+ is related to the 5-meson mass and 
the photon energy, = Mb — ^E^, and the measurement of its spectrum leads 
directly to the extraction of the leading hadronic structure function, called the 
shape function [15, 16, 17]. The P+ spectrum in semileptonic B — > Xul~P decays, 
on the other hand, enables us to determine \Vub\ [34, 85, 93], but this requires 
a precise knowledge of the shape function. One approach for measuring |V^5| is 
to first extract the shape function from the B — Xs'~f photon spectrum, and 
then to use this information for predictions of event distributions in P — > X^ I'l'- 
A comprehensive description of this program has been presented in chapter 4. 
Equivalently, it is possible to eliminate the shape function in B ^ X^ I'P decay 
rates in favor of the B photon-energy spectrum. This idea was first put 

forward in [16] and later refined in [67, 68, 69, 94]. Partial B — > X„ decay rates 
are then given as weighted integrals over the B — > Xsj photon-energy spectrum, 

exp. input theory exp. input 

where the weight function W{A, Pj^) is perturbatively calculable at leading power 
in Aqcd/''^6- a comparison of both sides of the equation determines the CKM 
matrix element \Vub\ directly. For the measurement of the left-hand side to be free 
of charm background, A must be less than M^/Mb ~ 0.66 GeV. However, the P+ 
spectrum in P — > X^ l~y decays displays many of the features of the charged- lepton 
energy spectrum, so that it is not inconceivable that the cut can be further relaxed 
for the same reasons that experimenters are able to relax the lepton cut beyond the 
charm threshold. We stress that for an application of relation (5.1) a measurement 
of the B — > Xg^ photon spectrum is needed only for > |(Ms — A) 2.3 GeV 
(or slightly lower, if the cut is relaxed into the charm region). This high-energy 
part of the spectrum has already been measured with good precision. 

Previous authors [16, 67, 68, 69, 94] have considered relations such as (5.1) in 
the slightly different form 

exp. input theory exp. input 

Normalizing the photon spectrum by the total^ rate TgiE^) as done in (5.1) has 

^Due to an unphysical soft-photon singularity, the total decay rate is commonly 
defined to include all events with photon energies above E^, = m{,/20 [32]. 
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several advantages. Firstly, it is a known fact that event fractions in 5 — Xg'-y 
decay can be calculated with better accuracy than partial decay rates (see [50] 
for a recent discussion), and likewise the normalized rate does not suffer from the 
relatively large experimental error on the total branching ratio. Secondly, relation 
(5.1) is independent of the CKM factor |V^t;,V^*|. Thirdly, unlike the total B Xg^ 
decay rate, the shape of the photon spectrum is rather insensitive to possible New 
Physics contributions [32], which could distort the outcome of a \ Vuh\ measurement 
via relation (5.2). Lastly, as we will see below, the weight function W{A,Pj^) 
possesses a much better perturbative expansion than the function H^(A, P+) = 
\VtbV^J\^W{^,P+)/Ts{E^). This last point can be traced back to the fact that 
most of the very large contribution from the Oi — O'j^ operator mixing in the 
effective weak Hamiltonian cancels in the theoretical expression for the normalized 
photon spectrum. 

In principle, any partial B — > X„/~z/ decay rate can be brought into the form 
(5.1), with complicated weight functions. The relation between the two P+ spectra 
is particularly simple, because the leading-power weight function is a constant at 
tree level. Experiments typically reject semileptonic i?-decay events with very low 
lepton energy. The effect of such an additional cut can be determined from [95]. 
Alternatively, it is possible to modify the weight function so as to account for a 
lepton cut, however at the expense of a significant increase in complexity. We will 
not pursue this option in this chapter. 

The weight function depends on the kinematical variable P+ and on the size 
A of the integration domain. It possesses integrable singularities of the form 
Q;"(yu) ln'^[m5(A — P+)//i^], with /c < n, in perturbation theory. Different strategies 
can be found in the literature concerning these logarithms. Leibovich et al. re- 
summed them by identifying jji^ with mb(A — P+) [67, 68]. The P+ dependence of 
the weight function then enters via the running coupling as{-\/ mh{^ — P+))- This 
is a legitimate choice of scale as long as (A — P+) has a generic value of order Aqcd; 
however, it is not a valid choice in the small region where mb{A — P_|_) ~ Aqqj-,. 
A key result underlying relation (5.1) is that, by construction, the weight function 
is insensitive to soft physics. Quark-hadron duality ensures that the region near 
the point P+ = A, which is without any physical significance, does not require 
special consideration after integration over P+. In the approach of [67, 68], the 
attempt to resum the above logarithms near the endpoint of the integral leads 
to integrals over unphysical Landau singularities of the running coupling in the 
nonperturbative domain. Hoang et al. chose to calculate the weight function in 
fixed-order perturbation theory at the scale fj, = rrib [94]. This leads to paramet- 
rically large logarithms, since A — P+ ^ rrih. In the present chapter we separate 
physics effects from two parametrically distinct scales, a hard scale fih ~ n^-b and 
an intermediate scale /Xj ~ ^/mbAqcD, so that we neither encounter Landau sin- 
gularities nor introduce parametrically large logarithms. The shape of the weight 
function is then governed by a perturbative expansion at the intermediate scale, 
a^iHi) ln^[mb(A — P+)//x^]. As will be explained later, the coefficients in this series. 



1 nc; 



as well as the overall normalization, possess themselves an expansion in as{iih)- 

The calculation of the weight function starts with the theoretical expressions 
for the P+ spectra in 5 — > and B —>■ Xgj decays, which are given as [95] 

+6(1 - y) 3^2{P+, y)+y 3'3{P+, y) } , (5.3) 



where 



(5.4) 



— * 

with P_ = Ex + \Px\ is a kinematical variable that is integrated over the available 
phase space. Expressions for the structure functions valid at next-to-leading 
order (NLO) in renormalizat ion-group (RG) improved perturbation theory and 
including first- and second-order power corrections can be found in [95] (see also 
[50, 62, 63]). Symbolically, they are written as H{iih) ■ J{iii) ®S{iii), where H[iih) 
contains matching corrections at the hard scale /z^. The jet function J{iJ,i), which 
is a perturbative quantity at the intermediate scale //j, is convoluted with a non- 
perturbative shape function renormalized at that same scale. Separation of the 
two scales fih and fii allows for the logarithms in matching corrections to be small, 
while logarithms of the form In Uhl i^ii which appear at every order in perturbation 
theory, are resummed in a systematic fashion and give rise to the RG evolution 
functions U{fih,IJ'i) and ar{fJ-h- fJ-i) [63]. 

The leading-power jet function J{p^,fii) entering the expressions for is uni- 
versal and has been computed at one-loop order in [62, 63]. The two- loop expres- 
sion for J has been obtained apart from a single unknown constant in [96], which is 
the two-loop coefficient of the local 5(p^) term. (This term was recently calculated 
in [97]). This constant does not enter in the two- loop result for the weight function 
W{A,P+) in (5.1). Due to the universality of the leading-power jet function, it 
is possible to calculate the complete 0{a'^{ni)) corrections to the weight function. 
However, the extraction of hard corrections at two-loop order would require multi- 
loop calculations for both decay processes, which are unavailable at present. As a 
result, we will be able to predict the A and P+ dependence of the weight function 
W{A, P+) at next-to-next-to-leading order (NNLO) in RG-improved perturbation 
theory, including exact two-loop matching contributions and three-loop running 
effects. However, the overall normalization of the weight function will have an 
uncertainty of 0{a'^{nh)) from yet unknown hard matching corrections. 

The total rate rs{E^) has been calculated in a local operator product expansion 
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and reads (including only the leading non-perturbative corrections) [50] 



UE^) = 32^4 mlml{nn)[Cf,inn)r x 

A2 Ci{fih) 



9ml C^^M 



(5.6) 



where Hg is the hard function of 3^^, and Hj- contains the remaining radiative 
corrections. We will present an explicit expression for this quantity at the end of 
Section 5.3 below. The hadronic correction proportional to A2/m^ cancels against 
an identical term in 3"^. Apart from two powers of the running 6-quark mass defined 
in the MS scheme, which is part of the electromagnetic dipole operator Oj^ in 
the effective weak Hamiltonian, three more powers of mi, emerge from phase-space 
integrations. To avoid the renormalon ambiguities of the pole scheme we use a low- 
scale subtracted quark- mass definition for m^. Specifically, we adopt the shape- 
function mass m5^(^*,/i*) [63, 83] defined at a subtraction scale = 1.5 GeV, 
which relates to the pole mass as 

mf' - mfil^., 1^.) + 1^. ^^^^ + .... (5.7) 

TT 

Throughout this chapter we will use m^^(/i*, //*) as the 6-quark mass and refer to 
it as mi, for brevity. The present value of this parameter is = (4.61 ± 0.06) GeV 
[96]. 



5.2 Calculation of the weight function 
5.2.1 Leading power 

The key strategy for the calculation of the weight function is to make use of 
QCD factorization theorems for the decay distributions on both sides of (5.1) and 
to arrange the resulting, factorized expressions such that they are both given as 
integrals over the shape function, J duj S{uj) gi{cu), with different functions gi for the 
left-hand and right-hand sides. Relation (5.1) can then be enforced by matching 
fl'LHS to g-RHS- Following this procedure, we find for the integrated P+ spectrum in 
B —>■ Xu decay after a series of integration interchanges 

A 1 A P+ 

dP+^oc Jdyy-^'^HM j dP+{MB - P+f j duymi,J{ymi,{P+ - u)) S{u) 
0^0 n 

A Mb -C^ 1 

^dujSiuj) J dq5q' Jdyy-'^H^{y)j (^In '^'^^y^^ +ln|/^ , 
* 

(5.8) 
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where = min(A, Mb — q), and [96] 

j In— = / dp^J{p^,iii) (5.9) 
\ f^i / Jo 

is the integral over the jet function. For the sake of transparency, we often sup- 
press the exphcit dependence on /^j and iih when it is clear at which scales the 
relevant quantities are defined. The function is a linear combination of the 
hard functions entering the structures Hi in (5.3), which in the notation of [95] is 
given by 

Hu{y, fih) = '2y\3-2y) H^iiy, ^ih)+l2y\l-y) Hu2{y, f^h)+2y^ H^ziy, f^h) • (5.10) 

RG resummation effects build up the factor in (5.8), where a = ardj^h, l^i) is 
the value of the RG-evolution function 

f^^du /-"^Mi) Yens («) 

avilJ'h, IJ'i) ^ — ^cusp{as{l^)) ^ - da-^^f—, (5.11) 

which depends only on the cusp anomalous dimension [27, 28]. The quantity a has 
its origin in the geometry of time-like and light-like Wilson lines underlying the 
kinematics of inclusive B decays into light particles. Our definition is such that 
a is a positive number for Hh > Hi and vanishes in the limit //^ — > //j. We find it 
convenient to treat the function ar{iJ>h, IJ'i) as a running "physical" quantity, much 
like OLsilj) or rnb(iJ.). Since the cusp anomalous dimension is known to three-loop 
order [98], the value of a can be determined very accurately. Note that three- loop 
accuracy in a (as well as in the running coupling as) is required for a consistent 
calculation of the weight function at NNLO. The corresponding expression is 

[To (3o 

where the expansion coefficients r„ and Pn oi the cusp anomalous dimension and 
/3-function can be found in appendix B. 

Instead of the jet function J itself, we need its integral j (In / fif, fii) in the 
second line of (5.8). Since the jet function has a perturbative expansion in terms 
of "star distributions", which are logarithmically sensitive to the upper limit of 
integration [2], it follows that j(L, /Xj) is a simple polynomial in L at each order 
in perturbation theory. The two-loop result for this quantity has recently been 
computed by solving the integro-differential evolution equation for the jet function 
[96]. An unknown integration constant of 0{alL^) does not enter the expression 
for the weight function. 
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We now turn to the right-hand side of (5.1) and follow the same steps that lead 
to (5.8). It is helpful to make an ansatz for the leading-power contribution to the 
weight function, W'^'^\A, P^), where the dependence on A is solely given via an 
upper limit of integration. To this end, we define a function f{k) through 



{Mb - P+r Jo 



A-P4. 



dkf{k) {Mb 



(5.13) 



This allows us to express the weighted integral over the B Xg^ photon spectrum 
as 



JdP^^W('\A,P^)^ 


A P+ 

JdP+{MB-P+fW^''\A,P+) J dujiUf, J{m.b{P+ — a))) S{u;) 

mb(Ag — cu — k) 



oc 



oc diu S{(2} 



dq5q / dkf{k)j ( In 



(5.14) 



Note that the jet function J (and with it j) is the same in semileptonic and radiative 
decays. The difference is that the argument of the jet function in (5.8) contains an 
extra factor of y, which is absent in (5.14). Comparing these two relations leads 
us to the matching condition 



/' 

^0 



(5.15) 



which holds to all orders in perturbation theory and allows for the calculation 
of M/(°)(A,P+) via (5.13). The main feature of this important relation is that 
the particular vahic of Q is irrelevant for the determination of f{k). It follows 
that, as was the case for the jet function J, the perturbative expansion of f{k) 
in as{ijLi) at the intermediate scale involves star distributions, and W^^\A, P+) 
depends logarithmically on (A — P+). At two-loop order it suffices to make the 
ansatz 



f{k) oc S{k) + Cf 



(^s{l^i) 



c^^ S{k) + c\ 



(^s{l^i) 

An 



An 

4'^(5(A;) + cf) 



(1) 
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(5.16) 



+ 2c?^ 
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where the star distributions have the following effect when integrated with some 
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smooth function 0(A;) over an interval 



rU /i\[M?/fn6] r 



In 



rribk 



m 



k 



0(0) 



A; 
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A sensitivity to the hard scale enters into f{k) via the appearance of 
Hu{y,f^h) in (5.15). Because of the polynomial nature of j{L,fii), all we ever 
need are moments of the hard function with respect to Iny. We thus define the 
master integrals 



Tn{a,lih)= / dyy ^" if„(y, /z^) In'' y ; Kia, /i^ 
Jo 



Tn{a,fih) 



(5.18) 



which can be calculated order by order in as{j^h). Therefore, the coefficients c^"' of 
the perturbative expansion in (5.16) at the intermediate scale have the (somewhat 
unusual) feature that they possess themselves an expansion in as{^h)- This is a 
consequence of the fact that, unlike the differential decay rates (5.3) and (5.4), 
the weight function itself does not obey a simple factorization formula, in which 
the hard correction can be factored out. Rather, as can be seen from (5.15), it is 
a convolution of the type W — H{^h) ® J{iJ>i)- To one-loop accuracy, the hard 



function Hy, reads 



Hyiy,lih) = 2y\3-2y) 



1 + 
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-41n^^ + 101n^-41ny 
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-4L2(l-y)---12 



TT 



Iny. 



(5.19) 



Explicit expressions for the quantities Tq, cj^\ and /i„ entering the distribution 
function f{k) will be given below. 



5.2.2 Sub leading power 

Power corrections to the weight function can be extracted from the corresponding 
contributions to the two P+ spectra in (5.3) and (5.4). There exists a class of 
power corrections associated with the phase-space pref actors {Mb — P+)'^ in these 
relations, whose effects are treated exactly in our approach, see e.g. (5.13). This 
is important, because these phase-space corrections increase in magnitude as the 
kinematical range A over which the two spectra are integrated is enlarged. One 
wants to make A as large as experimentally possible so as to increase statistics 
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and justify the assumption of quark-hadron duality, which underhes the theory of 
inclusive B decays. 

The remaining power corrections fall into two distinct classes: kinematical 
corrections that start at order as and come with the leading shape function [32, 2], 
and hadronic power corrections that start at tree level and involve new, subleading 
shape functions [51, 52, 66, 73, 74, 75, 76, 99]. Because different combinations 
of these hadronic functions enter in B ^ X„ l^u and B Xg'-f decays, it is 
impossible to eliminate their contributions in relations such as (5.1). As a result, at 
0{AQCY)/mi,) there are non-perturbative hadronic uncertainties in the calculation 
of the weight function W{A,P^), which need to be estimated before a rehable 
extraction of \Vub\ can be performed. For the case of the charged-lepton energy 
spectrum and the hadronic invariant mass spectrum, this aspect has been discussed 
previously in [52, 74] and [99], respectively. 

Below, we will include power corrections to first order in Aqcd/''^()- Schemat- 
ically, the subleading corrections to the right-hand side of (5.1) are computed 
according to ri^^ ~ W^°^ ® dr'-^^/dP^ + W^^^ ® dr''^^/dP+, where the superscripts 
indicate the order in l/rrii, power counting. The power corrections to the weight 
function, denoted by W^^\ are derived from the mismatch in the power corrections 
to the two decay spectra. The kinematical power corrections to the two spectra 
are known at 0{as), without scale separation. We assign a coupling as{fi) to these 
terms, where the scale p, will be chosen of order the intermediate scale [95]. At 
first subleading power the leading shape function is convoluted with either a con- 
stant or a single logarithm of the form ln[(P+ —6j)/[Mb — -?+)], and we have (with 
n = 0,1) 



f dv r r + 

dP+-^ 3 asifi) dP+{MB-P+f dujS{uj)hi 

Jo Jo Jo 



Mb-P+ 

4i„n 



= a,(/2) / dujS{Cj) / dk{MB-Cj -kfhi 
Jo Jo 



MB-LO-k 

(5.20) 



and similarly for the photon spectrum. On the other hand, the weighted integral 
in (5.1) also contains terms where the photon spectrum is of leading power and 
the weight function of subleading power, 

rA 7-p pA 

/ dP+^W'^''^^\A,P+)3 / dP+{MB-P+fS{P+)W'^''^^\A,P+). (5.21) 
Jo Jo 

Therefore the kinematical corrections to the weight function must have the form 
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and a straightforward calculation determines the coefficients A and B. The 
hadronic power corrections to the weight function, W^^'^^^^^ can be expressed in 
terms of the subleading shape functions t{uj), u{lj), and v{uj) defined in [66]. These 
terms are known at tree level only, and at this order their contribution to the weight 
function can be derived using the results of [95]. 



5.3 Results 



Including the first-order power corrections and the exact phase-space factors, the 
weight function takes the form 
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2pi{d) In 



Mn-P^ 



47r 
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Mb-P+ 2(l-a)(3-a) 



(5.23) 



+ 2p2(^) 



4(l-a)(A-P+) + 2(4-3a) 



i{P^ 



+ (4-a; 



SiP+) 



+ (8-13a + 4a2)5^ 
S{P+) 



m 



S{P^ 



Mb - P+ S{P+) 



where A = MB—Tni, is the familiar mass parameter of heavy-quark effective theory, 
and (5 = (A — P^)/{Mb — -P+). The first line denotes an overall normalization, 
the next three lines contain the leading-power contributions, and the remaining 
expressions enter at subleading power. The different terms in this result will be 
discussed in the remainder of this section. 

For the leading-power terms in the above result we have accomplished a com- 
plete separation of hard and intermediate (hard-collinear) contributions to the 
weight function in a way consistent with the factorization formula W = H ^ J 
mentioned in the previous section. The universality of the shape function, which 
encodes the soft physics in both B — > Xg^ and B l~P decays, implies 

that the weight function is insensitive to physics below the intermediate scale 
/^j ~ ^^foAqcD- In particular, quark-hadron duality ensures that the small re- 
gion in phase space where the argument rrih^A — P+) of the logarithms scales as 
^QCD smaller does not need special consideration. At a technical level, this can 
be seen by noting that the jet function is the discontinuity of the coUinear quark 
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propagator in soft-coUinear effective theory [22, 63], and so the P+ integrals can 
be rewritten as a contour integral in the complex plane along a circle of radius 



5.3.1 Leading power 



The leading-power corrections in the curly brackets in (5.23) are determined com- 
pletely at NNLO in RG-improved perturbation theory, including three-loop run- 
ning effects via the quantity a in (5.11), and two-loop matching corrections at the 
scale as indicated above. To this end we need expressions for the one-loop co- 
efficients Cn^ including terms of 0{as{iJih})-, while the two- loop coefficients Cn are 
needed at leading order only. We find 
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(5.25) 



As always (7^^ = 4/3, (7^ = 3, and Po — ll—2nf/Sis the first coefficient of the QCD 
/3-function. The term proportional to in the expression for c^^ arises because 
of the elimination of the pole mass in favor of the shape-function mass, see (5.7). 
Since the logarithms ln['mh{A — P^)/ij'^] in (5.23) contain mj, all coefficients except 
Cn""* receive such contributions. However, to two-loop order only Cq^"* is affected. 
Next, the corresponding expressions for the hard matching coefficients hi are 
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calculated from (5.18). To the required order they read 
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(5.26) 

where the n-th derivative of the polygamma function. Because of the 

exact treatment of the phase space there are corrections to the logarithms in (5.23), 
which are finite-order polynomials in the small ratio 5 = (A — P+)/(Mb — P+). 
Explicitly, 
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(5.27) 



This concludes the discussion of the leading-power expression for the weight func- 
tion. 



5.3.2 Subleading power 

The procedure for obtaining the kinematical power corrections to the weight func- 
tion has been discussed in Section 5.2.2. For the coefficients A and B in (5.23) we 
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(5.28) 



Here Ci{fih) denote the (effective) Wilson coefficients of tlie relevant operators in 
the effective wealc Hamiltonian, which are real functions in the Standard Model. 
The variable z — {mc/mbY enters via charm-loop penguin contributions to the 
hard function of the B — > Xg^ photon spectrum [32] , and 
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Furthermore we need the integrals over k in (5.22), which encode the phase-space 
corrections. They give rise to the functions 

Ia{5) = 1-25 + 25'' -5^ + \5\ 

5 



Ib{5) = /A(5)ln 
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55 5 5 15 20 ^ ' 



The hadronic power corrections come from subleading shape functions in the 
theoretical expressions for the two decay rates. We give their tree-level contribu- 
tions to the weight function in the last two hues of (5.23), where 5' denotes the 
leading shape function, and u, v are subleading shape functions as defined in 
[66]. For completeness, we also include a contribution proportional to resulting 
from finite-mass effects in the strange-quark propagator in 5 — ^ Xs'~f decays. For 
rus = O(Aqcd) these effects are formally of the same order as other subleading 
shape-function contributions [100], although numerically they are strongly sup- 
pressed. The appearance of subleading shape functions introduces an irreducible 
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hadronic uncertainty to a \\ub\ determination via (5.1). In practice, this uncer- 
tainty can be estimated by adopting different models for tlie subleading sliape 
functions. Tliis will be discussed in detail in Section 5.4.4 below. Until then, let us 
use a "default model" , in which we assume the functional forms of the subleading 
shape functions t{uj), u{uj), and 'O(a)) to be particular linear combinations of the 
functions S'{uj) and (A — cu) S{u;). These combinations are chosen in such a way 
that the results satisfy the moment relations derived in [66], and that all terms 
involving the parameter A cancel in the expression (5.23) for the weight function 
for any value of a. These requirements yield 

(5.32) 

and the last two hnes inside the large bracket in the expression (5.23) for the weight 
function simplify to 

AisF(a) S'{P+) ^ A^sf(«) an 4Aisp(a) 

J{E^ ~ ~ J^^^ ^^^-^ {Ms - P^Y ' ^^-^^^ 

where 

2 (2 + a)Ai a(7-4a)A2 , 

AssF(a)- 3(i_«)(3_„)+2(l-a)(3-a) ""^ " ^^"^^^ 

Here Ai and A2 = \{M'^t — M^) are hadronic parameters describing certain B- 
meson matrix elements in heavy- quark effective theory [54]. The strange- quark 
mass is a running mass evaluated at a scale typical for the final-state hadronic 
jet, for which we take 1.5 GeV. As mentioned above, the numerical effect of the 
strange-quark mass correction is small. For typical values of the parameters, it 
reduces the result for Ag3p by about 10% or less. The expression on the right-hand 
side in (5.33) is equivalent to that on the left-hand side after the integration with 
the photon spectrum in (5.1) has been performed. It has been derived using the 
fact that the normalized photon spectrum is proportional to the shape function 
S{P+) at leading order. Note that the second term in the final formula is power 
suppressed with respect to the first one. It results from our exact treatment of 
phase-space factors and thus is kept for consistency. 

5.3.3 Normalization 

Finally, let us present explicit formulae for the overall normalization factor in 
(5.23). The new ingredient here is the factor Tq, which is defined in (5.18). At 
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one-loop order we find 
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When the product of Tq with the quantity [50] 
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(5.36) 



from the total B — > Xgj decay rate is consistently expanded to 0{as{fJ,h)), the 
double logarithm cancels out. Here 5* = 1 — 2E^/mb — 0.9, and the functions 
fij{S*) capture effects from operator mixing. 



5.4 Numerical results 



We are now in a position to explore the phenomenological implications of our 
results. We need as inputs the heavy-quark parameters A2 = 0.12GeV^, Ai = 
(-0.25 ± 0.10) GeVl and the quark masses rub = (4.61 ± 0.06) GeV [96], m, = 
(90 ± 25)MeV [101, 102], and rric/mb = 0.222 ± 0.027 [50]. Here rrib is defined 
in the shape-function scheme at a scale //* = 1.5 GeV, is the running mass 
in the MS scheme evaluated at 1.5 GeV, and rric/mb is a scale invariant ratio of 
running masses. Throughout, we use the 3- loop running coupling normalized to 
as{Mz) = 0.1187, matched to a 4-fiavor theory at 4.25 GeV. For the matching 
scales, we pick the default values fi'^^^ = mbj \pl and ixf^^ = p,^^^ = 1.5 GeV, which 
are motivated by the underlying dynamics of inclusive processes in the shape- 
function region [63, 95]. 

In the remainder of this section we present results for the partial decay rate 
r„(A) computed by evaluating the right-hand side of relation (5.1). This is more 
informative than to focus on the value of the weight function for a particular choice 
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Figure 5.1: Residual scale dependence of r„(A) for A = 0.65 GeV. The dashed hne 

depicts the variation of about its default value mf,/\/2, the dash-dotted line the 
variation of /ij about 1.5 GeV, and the dotted line the variation of /2 also about 
1.5 GeV. The highlighted area shows the combined perturbative uncertainty. 



of P+. For the purpose of our discussion we use a simple model for the normalized 
photon spectrum that describes the experimental data reasonably well, namely 

(PJ-expf-.^) (5.37) 



r, dP+ T{h)A>^ ^ ^' K 
with A = 0.77 GeV and h = 2.5. 

5.4.1 Studies of the perturbative expansion 

The purpose of this section is to investigate the individual contributions to r„(A) 
that result from the corresponding terms in the weight function, as well as their 
residual dependence on the matching scales. For A = 0.65 GeV we find numerically 

r^(0.65 GeV) ^ 43.5 (1 + 0.158 [a.(w)]- 0.095 [a.(M,)] + 0.076 [a^(Mi)] 

- 0.037 {as{lJ.i)as{tih)\ + 0.009 [kin] - 0.043 [hadr]) = 46.5 . 

(5.38) 

The terms in parenthesis correspond to the contributions to the weight function 
arising at different orders in perturbation theory and in the l/rrib expansion, as 
indicated by the subscripts. Note that the perturbative contributions from the in- 
termediate scale are typically twice as large as the ones from the hard scale, which 
is also the naive expectation. Indeed, the two-loop Q;g(/ij) correction is numerically 
of comparable size to the one- loop as{^h) contribution. This confirms the impor- 
tance of separating the scales fii and fih- The contributions from kinematical and 
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Figure 5.2: Convergence of the perturbative series and residual scale dependence. 
Left: RG-improved results at LO (dotted), NLO (dashed), and NNLO (sohd) 
as a function of the scales /ij, fj,h, p, which are varied simultaneously about their 
default values. RIGHT: Fixed-order results at tree- level (dotted) and one-loop 
order (dashed). 



hadronic power corrections turn out to be numerically small, comparable to the 
two-loop corrections. 

The weight function (5.23) is formally independent of the matching scales /x/j. 
Hi, and ft. In Figure 5.1 we plot the residual scale dependence resulting from 
the truncation of the perturbative series. Each of the three scales is varied inde- 
pendently by a factor between l/\/2 and \/2 about its default value. The scale 
variation of Hi is still as significant as the variation of /ih, even though the former 
is known at NNLO and the latter only at NLO. We have checked analytically that 
the result (5.23) is independent of /ij through two-loop order, i.e. the residual scale 
dependence is an 0(Q;^(/ij)) efi:ect. In order to obtain a conservative estimate of 
the perturbative uncertainty in our predictions we add the individual scale depen- 
dencies in quadrature. This gives the gray band shown in the figure. 

Figure 5.2 displays the result for r„(0.65 GeV) at different orders in RG- 
improved perturbation theory. At LO, we dismiss all terms including the kine- 
matical power corrections; however, leading logarithms are still resummed and give 
rise to a non-trivial dependence of Tq on the coefficient a. At NLO, we include the 
0{as{iJ,h)), 0{as{iJ,i)), and 0{as{p,)) contributions, but drop terms of order a^{iJ,i) 
or as{i-ii) as{fj.h). At NNLO, we include all terms shown in (5.23). In studying the 
different perturbative approximations we vary the matching scales simultaneously 
(and in a correlated way) about their default values. Compared with Figure 5.1 
this leads to a reduced scale variation. The gray bands in Figure 5.2 show the total 
perturbative uncertainty as determined above. While the two-loop NNLO contri- 
butions arc sizable, we observe a good convergence of the perturbative expansion 
and a reduction of the scale sensitivity in higher orders. The right-hand plot in 
the figure contrasts these findings with the corresponding results in fixed-order 
perturbation theory, which are obtained from (5.23) by setting /ih = /li = p, = /j, 
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and truncating the series at 0{as{iJ>)) for consistency. We see that the fixed-order 
results are also rather insensitive to the value of fi unless this scale is chosen to 
be small; yet, the predicted values for are significantly below those obtained in 
RG-improved perturbation theory. We conclude that the small scale dependence 
observed in the fixed-order calculation does not provide a reliable estimator of the 
true perturbative uncertainty. In our opinion, a fixed-order calculation at a high 
scale is not only inappropriate in terms of the underlying dynamics of inclusive 
decay processes in the shape-function region, it is also misleading as a basis for 
estimating higher-order terms in the perturbative expansion. 



5.4.2 Comments on the normalization of the photon spec- 
trum 

We mentioned in section 5.1 that the use of the normalized photon spectrum is 
advantageous because event fractions in i? — >• Xsj decay can be calculated more 
reliably than partial decay rates. In this section we point out another important 
advantage, namely that the perturbative series for the weight function VF(A, P+) 
is much better behaved than that for W^(A, P+). The difference of the two weight 
functions lies in their normalizations, which are 

W{A,P+)^mlTo{a,i,,)Hr{ii,), W{A,P+) « ^"^^'^'^^ 



[C^^{piHWfn'M\HsM\'' 

(5.39) 

Here Hg is the hard function in the factorized expression for the structure function 
3^^ in (5.4), which has been derived in [50]. Note that the two weight functions 
have a different dependence on the 6-quark mass. In the case of W, three powers of 
rrib enter through phase-space integrations in the total decay rate Ts{E^:), and it is 
therefore appropriate to use a low-scale subtracted quark-mass definition, such as 
the shape-function mass. In the case of W, on the other hand, two powers of the 
running quark mass ffih{f^h) enter through the definition of the dipole operator Oy^, 
and it is appropriate to use a short-distance mass definition such as that provided 
by the MS scheme. In practice, we write frib{iih) as mb{mf,) times a perturbative 
series in as{nh)- 

The rnost pronounced eff'ect of the difference in normalization is that the weight 
function W receives very large radiative corrections at order as{nh), which range 
between —68% and —43% when the scale /ih is varied between mi, and mb/2. This 
contrasts the well-behaved perturbative expansion of the weight function W, for 
which the corresponding correct ionsjvary between —11% and —7%. In other words, 
the hard matching corrections for W are about six times larger than those for W. 
Indeed, these corrections are so large that in our opinion relation (5.2) should not 
be used for phenomenological purposes. 

The different perturbative behavior of the hard matching corrections to the 
weight functions is mostly due to the mixing of the dipole operator O^^ with other 
operators in the effective weak Hamiltonian for B — > decay. In order to 



1 on 



60 




0.2 0.4 0.6 0.8 

A [GeV] 



Figure 5.3: Perturbative uncertainties on r„(A) encountered when using the weight 
function W{A,P^) (dark gray) or W{A,Pj^) (hght gray). In the latter case, the 
normaUzation of the photon spectrum is chosen such that the two predictions agree 
at A = 0.65 GeV for central values of the matching scales. 



illustrate this fact, consider the one-loop hard matching coefficients defined as 

Wo,l + k^^ + ..., Wo,l + k^^ + .... (5.40) 

TT TT 

With our default scale choices we have k = —2.32 + 1.13 = —1.19, where the second 
contribution (+1.13) comes from operator mixing, which gives rise to the terms 
in the last two lines in (5.36). For the weight function W, this contribution has 
the opposite sign than the othei^terms, so that the combined jvalue of k is rather 
small. For the weight function W, on the other hand, we find k = —1.58 — 5.52 = 
—7.10. Here the contribution from operator mixing is dominant and has the same 
sign as the remaining terms, thus yielding a very large value of k. Such a large 
0{as) correction was not observed in [94], because these authors chose to omit 
the contribution from operators mixing. Note that at a higher scale = f^b, as 
was adopted in this reference, the situation is even worse. In that case we find 
k = -2.81 + 1.19 = -1.62 and k = -0.66 - 9.13 = -9.79. 

A visualization of the perturbative uncertainty is depicted in Figure 5.3, where 
predictions for r„(A) are shown using either (5.1) or (5.2). In each case, the 
error band is obtained by varying the different scales about their default values, 
/i„ G [fif^y V^, fi^^^ V2], and adding the resulting uncertainties in quadrature. The 
dark-gray band bordered by solid lines denotes the perturbative uncertainty of 
predictions when using the normalized photon spectrum, as in (5.1). (At the point 
A = 0.65 GeV this uncertainty is identical to the gray band depicted in Figures 5.1 
and 5.2). The light-gray band bordered by dashed lines corresponds to the use of 
the absolute photon spectrum, as in (5.2). The difference in precision between the 
two methods would be even more pronounced if we used the higher default value 
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I^h — iT^b for the hard matching scale. Obviously, the use of the normalized photon 
spectrum will result in a more precise determination of \Vub\- 



5.4.3 Comments on /3oa^ terms and scale separation 

The separation of different momentum scales using RG techniques, which is one of 
the key ingredients of our approach, is well motivated by the dynamics of charm- 
less inclusive B decays in the shape-function region. Factorized expressions for the 
S-decay spectra involve hard functions renormalized at /ih multiplied by jet and 
shape functions defined at a lower scale. While physics at or below the interme- 
diate scale is very similar for B X„ l~P and B — > Xg'y (as is manifested by the 
fact that the leading shape and jet functions are universal), the physics at the hard 
scale in B ^ Xg^ decay is considerably more complicated than in scmilcptonic 
decay, and it might even contain effects of New Physics. Therefore it is natural to 
respect the hierarchy /ih ^ A*i and disentangle the various contributions, as done 
in this chapter. In fact, our ability to calculate the dominant two- loop corrections 
is a direct result of this scale separation. Nevertheless, at a technical level we can 
reproduce the results of a fixed-order calculation by simply setting all matching 
scales equal to a common scale, fih = fii = Ji = fi. In this limit, the expressions 
derived in this chapter smoothly reduce to those obtained in conventional pertur- 
bation theory. While factorized expressions for the decay rates are superior to 
fixed-order results whenever there are widely separated scales in the problem, they 
remain valid in the limit where the different scales become of the same order. 

In [94], the 0(/3oas) BLM corrections [103] to the weight function 1^(A,P+) 
in (5.2) were calculated in fixed-order perturbation theory. For simplicity, only 
the contribution of the operator O-j^ to the B — > Xs^ decay rate was included 
in this work. We note that without the contributions from other operators the 
expression for W is not renormalization-scale and -scheme invariant. Neglecting 
operator mixing in the calculation of W is therefore not a theoretically consistent 
approximation. However, having calculated the exact NNLO corrections at the 
intermediate scale allows us to examine some of the terms proportional to (5Qa1{iJLi) 
and compare them to the findings of [94]. In this way we confirm their results for 
the coefficients multiplying the logarithms ln"[mfe(A — P+)/|U,^] with n = 1,2 in 
(5.23). While the /SqQ;^ terms approximate the full two-loop coefficients of these 
logarithms arguably well, we stress that the two-loop constant at the intermediate 
scale is not dominated by terms proportional to Numerically we find 
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-31.8 + 38.8 ^° 
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+ 0(a2), (5.41) 



which means that the approximation of keeping only the BLM terms would over- 
estimate this coefficient by almost an order of magnitude and give the wrong sign. 
This shows the importance of a complete two-loop calculation, as performed in 
this chapter. 
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We believe that the perturbative approximations adopted in this chapter, i.e. 
working to NNLO at the intermediate scale and to NLO at the hard scale, are suffi- 
cient for practical purposes in the sense that the residual perturbative uncertainty 
is smaller than other uncertainties encountered in the application of relation (5.1). 
Still, one may ask what calculations would be required to determine the missing 
CK^{lJ>h) terms in the normalization of the weight function in (5.23), or at least the 
terms of order /9o<^s(/"/i)- ^'^^ ^^'^ of B ^ Xg'y decay, the contribution of 
the operator Oj^ to the normalized photon spectrum was recently calculated at 
two- loop order [104], while the contributions from other operators are known to 
0{/3oal) [105]. What is still needed are the two-loop corrections to the double 
differential (in P+ and y) S — > X„Z~P decay rate in (5.3). 

5.4.4 Subleading corrections from hadronic structures 

Due to the fact that different linear combinations of the subleading shape functions 
t{u;), u{cj), and v{co) enter the theoretical description of radiative and semileptonic 
decays starting at order Aqcd/'^6! the weight function cannot be free of such 
hadronic structure functions. Consequently, we found in (5.23) all of the above 
subleading shape functions, divided by the leading shape function S{uj). Our 
default model (5.32) for the subleading shape functions was chosen such that the 
combined effect of all hadronic power corrections could be absorbed into a single 
hadronic parameter Aggp. More generally, we define a function 5hadr(^) via (a 
factor 2 is inserted for later convenience) 

r„(A) = [r,(A)]def [1 + 25hadr(A)] , (5.42) 

where [r„(A)]def denotes the result obtained with the default model for the sub- 
leading shape functions. From (5.23), one finds that 

W(A) = I J\p+ {Mb - P+)'[2(4 - 3a) ht{P+) + (4 - a) h^{P+) + 

+(8 - 13a + 4a2) K{P+)] | |4(1 - a) (3 - a) ^ dP+ {Mb - P+f S{P+) 

(5.43) 

where we have used that, at leading order in as and Aqcb /iTib, the B — > Xg'y 
photon spectrum is proportional to {Mb — P+)^S{P+). In the relation above, 
ht{u}) = t{i2i) — [t(ci))]dcf etc. denote the differences between the true subleading 
shape functions and the functions adopted in our default model. By construction, 
these are functions with vanishing normalization and first moment. 

The above expression for 5hadr(A) is exact to the order we are working; however, 
in practice we do not know the precise form of the functions hi{Qj). Our goal is then 
to find a conservative bound, |5hadr(A)| < ehadr(A), and to interpret the function 
£hadr(A) as the relative hadronic uncertainty on the value of extracted using 
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Figure 5.4: Estimates for the hadronic uncertainty ehadr(A) obtained from a scan 
over models for the subleading shape functions. The dashed hues correspond to 
the individual results for the four hi{uj) functions suggested in [95]. The thick solid 
line, which covers one of the dashed lines, shows the maximum effect. 




relation (5.1). To obtain the bound we scan over a large set of realistic models 
for the subleading shape functions. In [95], four different functions /ij(ci)) were 
suggested, which can be added or subtracted (in different combinations) to each 
of the subleading shape functions. Together, this provides a large set of different 
models for these functions. To be conservative, we pick from this set the model 
which leads to the largest value of |5hadr(^)|- The integrand in the numerator in 
(5.43) is maximized if all three /ij(u)) functions are equal to a single function, whose 
choice depends on the value of A. In the denominator, we find it convenient to 
eliminate the shape function S{P^) in favor of the normalized photon spectrum. 
Working consistently to leading order, we then obtain 



, , 5-5a + a^ U{fih,fJ'ij 
ehadr(A) = — 3 max 



/■A 

/ dP+{MB-P+Yhi{P+) 
Jo 



l-a)(3-a) ml i , 1 dV, 

(5.44) 

where as before a = aY^Hh, fJ-i) ~ 0.12 for the default choice of matching scales, 
and U{iJ,h,iJ,i) fa 1.11 [95]. 

The result for the function ehadr(A) obtained this way is shown in Figure 5.4. 
We set the matching scales to their default values and use the model (5.37) for 
the photon spectrum, which is a good enough approximation for our purposes. 
From this estimate it is apparent that the effects of subleading shape functions 
are negligible for large values A 3> Aqcd and moderate for A ~ Aqcd, which is 
the region of interest for the determination of |V^b|. In the region A < 0.3 GeV 
the accuracy of the calculation deteriorates. For example, we find ehadr = 2.0% for 
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A = 0.65 GeV, ehadr = 4.8% for A = 0.5 GeV, and ehadr = 8.2% for A = 0.4 GeV. 
5.5 Conclusions 

Model-independent relations between weighted integrals of i? — and 13 — > 
l~P decay distributions, in which all reference to the leading non-perturbative 
shape function is avoided, offer one of the most promising avenues to a high- 
precision determination of the CKM matrix element \Vub\. In order to achieve a 
theoretical precision of better than 10%, it is necessary to include higher-order 
corrections in and Aqcd/'"T'6 in this approach. 

In this chapter, we have calculated the weight function VF(A, P+) in the relation 

— * 

between the hadronic P+ = Ex — \Px\ spectra in the two processes, integrated over 
the interval < P+ < A. Based on QCD factorization theorems for the differential 
decay rates, we have derived an exact formula (5.15) that allows for the calculation 
of the leading-power weight function to any order in perturbation theory. We have 
calculated the A- and P+-dependent terms in the weight function exactly at next- 
to-next-to-leading order (NNLO) in renormalization-group improved perturbation 
theory, including two-loop matching corrections at the intermediate scale Hi ~ 
y^m^AqcD and three-loop running between the intermediate scale and the hard 
scale fih ~ iTT-b- The only piece missing for a complete prediction at NNLO is 
the two-loop hard matching correction to the overall normalization of the weight 
function. A calculation of the a'^d-ih) term would require the knowledge of both 
decay spectra at two-loop order, which is currently still lacking. We also include 
various sources of power corrections. Power corrections from phase-space factors 
are treated exactly The remaining hadronic and kinematical power corrections are 
given to first order in Aqcd/ nib and to the order in perturbation theory to which 
they are known. 

A dedicated study of the perturbative behavior of our result for the weight 
function has been performed for the partial B — * X„/^z/ decay rate r„(A) as 
obtained from the right-hand side of relation (5.1). It exhibits good convergence 
of the expansion and reduced scale sensitivity in higher orders. We find that 
corrections of order Q;^(/Xi) at the intermediate scale are typically as important 
as first-order agdj^h) corrections at the hard scale. We have also seen that fixed- 
order perturbation theory significantly underestimates the value of r„(A), even 
though the apparent stability with respect to scale variations would suggest a 
good perturbative convergence. In order to obtain a well-behaved expansion in 
powers of as{nh), it is important to use the normalized photon spectrum in relation 
(5.1). A similar relation involving the differential B — > Xgj decay rate receives 
uncontrollably large matching corrections at the hard scale and is thus not suitable 
for phcnomcnological applications. At next-to-leading order in the l/rrib expansion, 
the weight function receives terms involving non-perturbative subleading shape 
functions, which cannot be eliminated. Our current ignorance about the functional 
forms of these functions leads to a hadronic uncertainty, which we have estimated 
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by scanning over a large set of models. We believe that a reasonable estimate of 
the corresponding relative uncertainty ehadr on \Vub\ is given by the solid line in 
Figure 5.4. 

Let us summarize our main result for the partial 13 — > Xul^v decay rate with 
a cut P+ < 0.65 GeV, which is close to the charm threshold M^/Mb, and present 
a detailed hst of the various sources of theoretical uncertainties. We find 



where the central value is derived assuming that the B — > Xgj photon spectrum 
can be accurately described by the function (5.37). The errors refer to the per- 
turbative uncertainty as estimated in Section 5.4.1, the uncertainty due to the 
ignorance about subleading shape functions as discussed in Section 5.4.4, the error 
in the value of the 6-quark mass, other parametric uncertainties from variations 
of rric, rUsi and Ai, and finally a 6% uncertainty in the calculation of the normal- 
ization of the photon spectrum [50]. To a good approximation the errors scale 
with the central value. The above numbers translate into a combined theoretical 
uncertainty of 4.4% on \Vuh\ when added in quadrature. 



r„(0.65GeV) = 

(46.5 ± 1.4 [pert] ±1.8 [hadr] ± 1.8 [mt] ± 0.8 [pars] ± 2.8 [norm] ) \Vub\^ pS 



-1 




(5.45) 
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Chapter 6 
Conclusions 

Inclusive charmless B decays offer the most precise methods for measuring \Vub\- 
The physics of these decays is rich and involves a delicate interplay of theory and 
experiment. In this work we have developed and described an elaborate theoretical 
machinery that can be applied to these decays. 

We have started our discussion by explaining how to relate the differential 
decay rates of ^ ^ Xsj and B —>■ IP to the hadronic tensor. We have reviewed 
the traditional (approximate) methods of calculating the hadronic tensor, namely, 
an expansion in the strong coupling constant and HQET based OPE. We have 
explained the problems that arise from the use of these methods, and argued that 
they can be overcome by the use of SCET, the appropriate effective field theory for 
charmless inclusive B decays. We therefore reviewed some of its basic elements. 
Using SCET one can prove that the hadronic tensor factorizes, at each order in 
I /nib, as a sum of products of a calculable hard function multiplied by a calculable 
jet function convoluted with a non perturbative shape function. 

At leading order in 1/ all these functions are unique. Wc have calculated the 
leading order hard and jet function by performing a two-step matching calculation 
QCD— >SCET— >HQET. We have resuumed large logarithms that arise in the cal- 
culation at next-to-leading order in renormalization-group improved perturbation 
theory. We have employed the operator product expansion to relate moments of 
the rcnormalized shape function with HQET parameters such as m^, A and /i^ 
defined in a new physical subtraction scheme, called the "shape function scheme". 
We have obtained an analytic expression for the asymptotic behavior of the shape 
function, which reveals that it is not positive definite. 

Having analyzed the leading power terms wc have considered the contributions 
of subleading shape functions. At tree-level, the results can be expressed in terms 
of forward matrix elements of bi-local light-cone operators. Four-quark operators, 
which arise at 0{as), were also included. We have shown that at tree level only 
three independent subleading shape functions are needed. 

Based on these calculations we turned, in the second part of this work, to ap- 
plications. In the first application, we have presented "state-of-the-art" theoretical 
expressions for the triple differential B — > X„ l^u decay rate and for the B — > Xg'j 
photon spectrum, which incorporate all known contributions and smoothly inter- 
polate between the "shape-function region" of large hadronic energy and small 
invariant mass, and the "OPE region" in which all hadronic kincmatical variables 
scale with Mb- The differential rates are given in a form which has no explicit 
reference to the mass of the b quark, avoiding the associated uncertainties. Depen- 
dence on TOft enters indirectly through the properties of the leading shape function, 
which can be determined by fitting the B — > Xg^ photon spectrum. This elim- 
inates the dominant theoretical uncertainties from predictions for B X,^l~u 
decay distributions, allowing for a precise determination of |Kt5|. In the shape- 
function region, we have factorized short- distance and long-distance contributions 
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at next-to-leading order in renormalization-group improved perturbation theory. 
The higher-order power corrections include effects from subleading shape func- 
tions where they are known. When integrated over sufficiently large portions in 
phase space, our results would reduce to standard OPE expressions up to yet un- 
known 0{al) terms. We have presented predictions for partial B — > Xy,l~i> decay 
rates with various experimental cuts. The elaborate error analysis that we have 
performed contains all significant theoretical uncertainties, including weak annihi- 
lation effects. We have suggested that the latter can be eliminated by imposing a 
cut on high leptonic invariant mass. 

In the second application we have derived a shape-function independent rela- 

— — * 

tion between the partial B ^ I'P decay rate with a cut on P+ = Ex — \Px\ < ^ 
and a weighted integral over the normalized B Xg'j photon-energy spectrum. 
We have calculated the leading-power contribution to the weight function at next- 
to-next-to-leading order in renormalization-group improved perturbation theory, 
including exact two-loop matching corrections at the scale Hi ~ a/^vV^cd- The 
overall normalization of the weight function is obtained up to yet unknown cor- 
rections of order Q;f(m(,). We have included power corrections from phase-space 
factors exactly, while the remaining subleading contributions have been included 
at first order in Aqcd /fnb. At this level unavoidable hadronic uncertainties enter, 
which we have estimated in a conservative way. The combined theoretical accu- 
racy in the extraction of \Vub\ is at the level of 5% if a value of A near the charm 
threshold can be achieved experimentally. 

The work presented here was already implemented by the " B-factories" : Babar 
and Belle. As of June 2006, Babar has measured \Vub\ using a combined sh and 
El cut [106], combined Mx and cut [107], and Ei cut [108]. Babar has also 
extracted the Mx spectrum for B — > Xg^ [5]. As of June 2006, Belle has measured 
\Vub\ using an Ei cut [109], Mx cut, combined Mx and cut, and cut [110]. 

We would like to conclude this work with an outlook. The reduction of the 
experimental errors encourages theorists to improve the theoretical predictions. As 
of June 2006, what theoretical developments are feasible? In general the formahsm 
developed in this work can be extended to higher orders in l/m^ and as- 

At leading power the jet function [97] and the partonic shape function [111] 
have already been calculated to second order in a^. For B — > Xg'j a next-to-next- 
to-leading order QCD calculation is underway, where the contribution of Oy^ is 
already known [104, 112]. Once the complete result is known, one could extract 
the hard function for B ^ Xgj from it. An "easier" calculation is to find the 
leading power hard function for B Xul^v- Having these ingredients would 
allow to extend the analysis of chapters 4 and 5 to a complete two loop order. 

Beyond leading power one would like to find a way to extract the subleading 
shape functions from experiment, or at least to constrain them. This would help 
to reduce a large source of uncertainty in the extraction of \Vub\- Another possible 
improvement is to find the subleading shape function beyond tree level. The anal- 
ysis performed in [75, 76] indicates that beyond tree level more subleading shape 
functions would arise. This calculation is related to the long overdue OPE calcula- 
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tion of the 0{as) power corrections corrections to -B — > Xcl~V decays. The latter 
calculation would allow to check the results of the former, once it is performed. 

Finally, as a result of the better control of the charm background, the exper- 
imental cuts are being pushed away from the shape function region and into the 
OPE region. In this intermediate region a multi scale OPE calculation, like the 
one performed in [50], is more appropriate. This calculation is yet to be done for 
semileptonic decays. 

To summarize, despite the impressive progress achieved in our understanding 
of charmless inclusive B decays, much is left to be done. Hopefully the research 
presented here would help make this progress possible. 
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Appendix A 

The hadronic tensor 

In chapter 1 we defined the current correlator T'^" as: 

T'^'^i J d^x e^^^ r{ J^'^(O), r{x)}. (A.l) 
Reinserting the complete set of hadronic states we have 

{B\T^'''\B) = i I d^xe''^'' 9{-x'^)^{B\J^^{id)\Xq){Xq\r{x)\B) 

'J V 



+ e{x') Y,{B\r{x)\x,,,){x,,,\f^mB) 



bbq 



(A.2) 



Following [12] we write the steps function as: 



1 p—iijjx^ 

^(±a;0) = ^— / duj- . (A.3) 

(these identities can be proved by contour integration, or by integrating the integral 
representation of the delta function). Translation invariance allows us to write: 

{X\r{x)\B) = e'(^^-^^)-"(X|r(0)|5) 

{B\r{x)\X) = e*(-P^-f'^)-"(^|J"(0)|X). (A.4) 

Integrating over x and uj we find that in the rest frame of the B meson: 

^ —Mb + -c/x + q — le 

^ Mb - + (/" + !€ ' ' ^ ' 

■^bbq 

Taking the imaginary part, we find using equation (1.33) 

^ Xq 

+ l^^Y.W''mXm){Xm\J^'mB)-{2n?5\PB-Px + q) (A.6) 

Xbbq 

For B decays the second sum does not contribute, since the energy of X^b^ is always 
large than Mb- 
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Using this expression for the hadronic tensor, it is easy to prove the following 
symmetry of the hadronic tensor: 

lyM^ ^ (W''^')* . (A.7) 

Using the definition of the Hermitian conjugate, we have for each term in the first 
sum in (A. 6): 

{{B\j'"'\X,){X,\J^\B)y = {X,\r\B){B\J^^'\X,) = {B\j''^\X,){X,\r\B) (A.8) 

A similar relation holds for the second sum in (A. 6). Combining them gives us 
(A.7). 
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Appendix B 

Perturbative Expressions 

B.l Anomalous dimensions 

Here we list the known perturbative expansions of the /^-function and relevant 
anomalous dimensions. We work in the MS scheme and define 



din// 



n=0 



cusp 



n=0 



n=0 



fB.ll 



as the expansion coefficients for the /3-function, the leading-order SCET current 
anomalous dimension, and the cusp anomalous dimension. To three-loop order, 
the /9-function reads [113] 



132 = 



2857 

~5r 



C1+ C]. 



205 
~18 



fc'^-fcAnf-2Ci.nf, 



1415 

~5r 



Ciy,+ (^^C,+ '^^C?)n}, 

(B.2) 



where n/ = 4 is the number of fight flavors, = 3 and Cp = 4/3. The three-loop 
expression for the cusp anomalous dimension has recently been obtained in [98]. 
The coefficients read 



r2 = i6Cf 



Ti = 8Cf 



67 _ 7r^\ _ 5^ ■ 
18 6/^9^-^ 



/245 677r^ IItt^ 11 . , 

55 



24 



27 



209 bTT^ 7 , , ^ 
T08 + ^"3^^^^^''^ 



(B.3) 



The SCET anomalous dimension 7 is explicitly known only to one-loop order. 
However, the two-loop coefficient can be extracted by noting that 7 is related to 
the axial-gauge anomalous dimension in deep inelastic scattering [50]. The result 
is 



7o 



-5Cf 
-8Cf 

125 vr^ 



(l^-?-C3)c..( 



/154 9 77r2 _ 11 
432" ^ 18" ^ T ' ^ 



B.2 Evolution factor 



The exact expression for the evolution factor reads 

lnC/(//ft,//j) = 25'r(/Xft,/Xj) - 2ar(/^fe,/^j) In 2ay (//^, /ij) , (B.5) 

where the functions of the right-hand side are solutions to the renormalization- 
group equations 



d 



din II 

d 

dln/j, 



5r(z/,/x) 



-rcusp(as(Ai)) In - , 



din /J, 



ay(u, ii) = -7'(q;,(/x)) , (B.6) 



with boundary conditions S'(z/, /i) = etc. at ;U = z/. These equations can be 
integrated using that d/dln/i — (3{as) d/das- The solutions are 



da' 



Or 



da ^^^'^ 



a) 



(B.7) 

and similarly for a^i. 

Next, we give explicit results for the Sudakov exponent Sy and the functions 
ar and in (B.5) at next-to-leading order in renormalization-group improved 
perturbation theory. We obtain 



In 



+ ... 



(B.8) 



and similarly for a^. The ncxt-to-leading order expressions for the Sudakov expo- 
nent Sr contains the three-loop coefficients P2 and T2- With r = as{n)/as{i'), it 
reads 



Tr 



An 



+ 



47r 



'-(l--- Inr) + (^-^) (1 - r + lnr) + 

r J \ro (3oJ ^ ^ 2(30 



(^0^0 /3o 



[1 — r + rlnr) 



P! P2 Pi^i ^ T2\ {1-rf 
Po Po Po^o To 



P! 


P2 


Pi 


Po 


+ . 


} 



1 — r) In r 

(B.9) 



The next-to-lcading-logarithmic evolution factor U {fih, ^i) can be obtained by com- 
bining the above expressions according to (B.5) and expanding out terms of order 

Ots- 
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Appendix C 

Partially integrated decay rates 

With the exception of the combined cut on the lepton energy Ei and the hadronic 
quantity studied in Section 4.5.5, all other partial rates investigated in our 
analysis can be derived by first integrating the triple differential decay rate (4.19) 
over the lepton energy Ei > Eq and P < P™^^ analytically, where the quantity 
pmax (^g^j^^ principle even Eq) may depend on the value of P+. The remaining 
integration over P+ is then performed numerically. In such a situation, we need to 
evaluate the partially integrated decay rate 



Changing variables from P to y defined in (4.21), the constraint P < p™^^ 
translates into the integration domain < y < ^max, where in analogy to (4.21) 
we define 

_ pma. - p^ _ p^^^ _ p^ _ _ 2Eo 

- Mb-P+ ' Mb'P+ Mb-P+' ^ ^ 

From the phase-space relation (1.5) it follows that a cut on the lepton energy has 
no effect if yo > ymax- The result of performing the integrations in (C.l) can be 
written as 

dru{yms.^,yo) _ / r;^(ymax) ; ymax<yo, 



dP+ \ r^{yo) + r^ ; ymax>yo, 



(C.3) 



where 



2 IT/ _ 12 



X fUyy^-^-- [(3 - 2y) + 6(1 - y) 7^ + y 3^3] , 
Jo 

= ^^{MB-P^fU{^,n,^J^^) p^^dyy-^'^-y. 



X [{Qyil + 1/0) - - i/o(3 + 2yo)) 

+Qy{l -y)'J2 + yoCiy - 2|/o) H • (C.4) 

When the kinematical power corrections in (4.26) are expanded as in (4.28) and 
(4.29), the resulting integrals over y can be expressed in terms of the master 
functions In{b,z) given in eq. (86) of [63]. The resulting expressions are used to 
obtain the numbers in the various tables in Section 4.5. 

We now list the values of t/o and t/max for fhc different cuts studied in Section 4.5. 
Whenever a cut Ei > Eq on the charged-lepton energy is applied, we have 
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For an additional cut P+ < Ap, we have j/max = 1 and < P+ < min(Ap, Mb — 
2Eo). For a cut on hadronic invariant mass, Mx < Mq, we have 

_ min(M..M„VP,)-P, 

Mb-P+ ^ ^ 

and < P+ < min(Mo, — 2Eq). For a cut on leptonic invariant mass, > Qq, 
we have ^ 

and < P+ < min(Ms — go, Mb — ^Eq). Finally, for the combined Mx-q^ cut we 
take the minimum of the previous two ymax values. 
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